Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,083 T→A 5.5% intergenic (+34/‑279) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,0830TA 5.5% 1441.1 / 48.8 454intergenic (+34/‑279)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (114/311);  new base A (4/21);  total (118/336)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 465 total aligned reads displayed.

GCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCT  >  NC_002947/954934‑955231
                                                                                                                                                     |                                                                                                                                                             
  caGTGGGCCCGCAAGGTGTTCGGCTGCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGCCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACtgn                                                                                                                                                             <  1:4946648/149‑2 (MQ=255)
        gCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGGAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCggg                                                                                                                                                      <  1:4833063/150‑1 (MQ=255)
        gCCCGCAAGGTGTTCGACTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCggg                                                                                                                                                      <  2:7536651/150‑1 (MQ=255)
         cccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCgggg                                                                                                                                                     >  1:1834630/1‑150 (MQ=255)
          ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑t                                                                                                                                                 >  2:7189641/1‑150 (MQ=255)
          ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑t                                                                                                                                                 >  1:4653711/1‑150 (MQ=255)
          ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑t                                                                                                                                                 <  1:5919243/150‑1 (MQ=255)
           cGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑tt                                                                                                                                                <  1:7938842/150‑1 (MQ=255)
           cGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑tt                                                                                                                                                <  2:1606619/150‑1 (MQ=255)
            gCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑ttt                                                                                                                                               <  2:6070504/150‑1 (MQ=255)
              aaGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑ttttc                                                                                                                                             <  2:3571766/150‑2 (MQ=255)
                  tgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑ttttttttt                                                                                                                                         >  2:3821820/1‑150 (MQ=255)
                  tgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑ttttcgtga                                                                                                                                         <  2:4570832/150‑6 (MQ=255)
                  tgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑ttttcgtga                                                                                                                                         <  2:1412486/150‑6 (MQ=255)
                  tgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCACAGTGGCCGGGG‑‑‑ttttcgtga                                                                                                                                         <  2:6333152/150‑6 (MQ=255)
                   gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTg                                                                                                                                        <  2:6672182/150‑1 (MQ=255)
                   gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTg                                                                                                                                        <  2:7954120/150‑1 (MQ=255)
                   gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTg                                                                                                                                        <  2:6260005/150‑1 (MQ=255)
                   gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTg                                                                                                                                        <  2:6877416/150‑1 (MQ=255)
                   gtTCGGCTTCTCCATCCTCACCGTCACGGCCCTGAGCGTGATGAGGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTg                                                                                                                                        <  2:8182040/150‑1 (MQ=255)
                    ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTgg                                                                                                                                       <  2:4321445/150‑1 (MQ=255)
                    ttCGGCTTCTACATCCTCACCGTCACAGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTgg                                                                                                                                       >  1:4159288/1‑150 (MQ=255)
                    gtCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTgg                                                                                                                                       <  1:5360269/149‑1 (MQ=255)
                     tCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGAGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑ATATATTTTggg                                                                                                                                      <  2:5711750/150‑1 (MQ=255)
                        gCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGt                                                                                                                                   <  2:7290951/150‑1 (MQ=255)
                        gCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGt                                                                                                                                   >  2:4138247/1‑150 (MQ=255)
                        gCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGt                                                                                                                                   >  1:6891619/1‑150 (MQ=255)
                          ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑AATTTTTTTGGGGGTAt                                                                                                                                 <  2:229844/150‑1 (MQ=255)
                                         gTCACTGCCCTGAGCGTGTTGATGGGTGTGGACTGCGAGACCGCTGCGGCTGTGCTGTTGACCTACGCACGCTGAAACTGGTGCCTGTTTCACGAAATGCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TGTTCCTGTGCCggg                                                                                                             <  2:6500984/150‑2 (MQ=255)
                                              tGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCtt                                                                                                        >  2:2833754/1‑150 (MQ=255)
                                               gCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTc                                                                                                       <  1:3821820/150‑1 (MQ=255)
                                                    gAGCGTGATGAGGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCggg                                                                                                  <  2:4871445/150‑1 (MQ=255)
                                                        gTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTaaa                                                                                              <  2:7479227/150‑1 (MQ=255)
                                                         tgatgaTGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCgggggaac                                                                                             >  1:892575/1‑145 (MQ=255)
                                                         tgatgaTGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                             <  2:3282308/150‑1 (MQ=255)
                                                          gatgatGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCggg                                                                                                                                                      >  2:6417509/1‑100 (MQ=255)
                                                          gatgatGGGTGTGGACAGCCAGACCGCTGCGGAGGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGT‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                             <  2:2593233/150‑1 (MQ=255)
                                                               tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                       <  2:4026639/150‑1 (MQ=255)
                                                                gggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                      <  1:6568678/150‑4 (MQ=255)
                                                                gggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                      <  1:1023013/150‑1 (MQ=255)
                                                                 ggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCACTTCGCGGGTAAACCCGcgaca                                                                                     >  2:1225808/1‑146 (MQ=255)
                                                                  gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                       <  1:6384908/150‑1 (MQ=255)
                                                                  gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCCACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                    <  2:8140098/150‑1 (MQ=255)
                                                                    gtgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACa                                                                                  <  1:7442941/150‑1 (MQ=255)
                                                                    gtgGACAGCCAGACCGCTGCGGATGTGCTGAGGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaagagg                                                                                  <  1:4578832/150‑8 (MQ=255)
                                                                     tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACag                                                                                 <  2:2210972/150‑2 (MQ=255)
                                                                      ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                                                                                <  2:6603035/150‑1 (MQ=255)
                                                                        aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGT‑‑TTTTTTTTTGGGGTTAT‑‑‑‑‑TTTTCCTGTGCCGGCCCGTGCCCGGGTAAACCCGAGCCAACAAgg                                                                               >  1:8075734/1‑150 (MQ=255)
                                                                        aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGCAAACCCGCTCCTACAAGGt                                                                              >  2:2738402/1‑150 (MQ=255)
                                                                         cAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGtt                                                                             <  1:6458784/150‑1 (MQ=255)
                                                                          cgCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACGGGGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaggaacc                                                                            <  2:5557665/149‑7 (MQ=255)
                                                                          aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                            <  1:8238716/150‑1 (MQ=255)
                                                                          aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                            <  2:192739/150‑1 (MQ=255)
                                                                          aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                            <  1:7763697/150‑1 (MQ=255)
                                                                            ccAGACCGCTGCGGATGTGCTGATGACCTTCGCTCGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                          <  2:4413610/150‑1 (MQ=255)
                                                                            ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                          <  2:8602627/150‑1 (MQ=255)
                                                                              aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCTGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                        <  1:2376828/150‑1 (MQ=255)
                                                                              aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                        <  1:3528088/150‑1 (MQ=255)
                                                                              aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                        <  2:5191334/150‑1 (MQ=255)
                                                                               gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCtt                                                                                                        <  2:3463147/117‑1 (MQ=255)
                                                                               gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                       <  1:6697649/150‑1 (MQ=255)
                                                                                   gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                          <  2:831353/143‑1 (MQ=255)
                                                                                    cTGCGGTTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGT‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTATGTGCTGGCCTCTTCGCGGTTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                 <  1:2535038/150‑1 (MQ=25)
                                                                                    cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                  <  2:7879332/150‑1 (MQ=38)
                                                                                    cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                  <  1:2668790/150‑1 (MQ=38)
                                                                                    cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                  <  2:505560/150‑1 (MQ=38)
                                                                                    cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                  <  1:6116315/150‑1 (MQ=255)
                                                                                     tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                 <  1:7189641/150‑1 (MQ=38)
                                                                                     tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                 <  1:7580227/150‑1 (MQ=38)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTTCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                <  1:7107253/150‑1 (MQ=255)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                <  2:5992223/150‑1 (MQ=38)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                <  1:5584584/150‑1 (MQ=38)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                <  1:4018328/150‑1 (MQ=38)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                <  1:6288793/150‑1 (MQ=38)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                <  1:3339819/150‑1 (MQ=255)
                                                                                      gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                <  2:7020658/150‑1 (MQ=255)
                                                                                       cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGT‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                <  1:3153873/150‑1 (MQ=37)
                                                                                       cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                               <  2:3034314/150‑1 (MQ=255)
                                                                                       cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                               <  1:989476/150‑1 (MQ=255)
                                                                                       cGGATGTGCTGATGACCGACGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACGGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTCT‑‑‑‑‑TTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGAACCTTGTAGGAGCggg                                                               <  2:8787656/150‑1 (MQ=25)
                                                                                       cGGATGTGCTGAGGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                               <  1:5240126/150‑1 (MQ=255)
                                                                                        ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGT‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                               <  1:50544/150‑1 (MQ=37)
                                                                                        ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TATTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                              <  1:8510897/150‑1 (MQ=255)
                                                                                         gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGAGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                  <  1:4354976/150‑1 (MQ=18)
                                                                                           tgtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTa                                                           <  2:7795629/150‑1 (MQ=255)
                                                                                            gtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc                                                          <  2:4469870/150‑1 (MQ=255)
                                                                                                      ccTACGCACGCTTAAACTGCTGCCTTTTTAACGAAAACCCCGGCCACTGTGCCGGGT‑‑‑TTTTTTTTTGGGGTTTT‑‑‑‑‑TTTTCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGTTCCTGTAGGTGGGGTTTTACCCGCGAagag                                                <  2:3702881/150‑1 (MQ=37)
                                                                                                       cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCGCCTACAAGGTTGCTGGTGGGGGGGGGTTTCCCGCGAAgggg                                               >  1:2105776/1‑150 (MQ=255)
                                                                                                          cgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCGCTCCGACAAGGGTCGGGGAGGGGCGGGGTTAACCGCGGAGAgggcc                                            >  1:1306187/1‑147 (MQ=255)
                                                                                                           gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTGAACCCCCGCCGAACAGGTTCCTGTAGGGGCGGGTGTGCCCCCGGAGAgggccg                                           >  1:5630536/1‑146 (MQ=255)
                                                                                                               gcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTCCTACAGGAACCTTGTTGGAGCggg                                                               <  1:1729125/126‑1 (MQ=25)
                                                                                                               gcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                               >  2:1729125/1‑126 (MQ=25)
                                                                                                                     cctgcggCCTGTGACACGAAAACCCCGGCCCCGGGGCGGGGT‑‑‑TTTTTTTTTGGGGGTTT‑‑‑‑‑TTTACCTGGGCCGGCAACTCCGCGGGTAAACCCGCTCCTACACGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                 <  1:3633301/149‑1 (MQ=21)
                                                                                                                     aCTGCTGCCTGTGACACGAACACCCCGGCCACTGTGCCGGGG‑‑‑TATTTTTTGGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                 <  2:5390547/150‑1 (MQ=37)
                                                                                                                                tttCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                            <  2:5630411/96‑1 (MQ=255)
                                                                                                                                       aaaaCCCCGGCACAGTGGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTGCGGGGGGGGCGGGGTTAACCGCGGAGTGGCCCCCACAGAAAAATAAAAATACCTGATTTTc               >  1:8683453/1‑150 (MQ=25)
                                                                                                                                         aaCCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAAACCGCTCCTACAAGGTTCCGGTTGGAGCGGGGTTACCCGCGGAGAGGCCAGCACAGCAAAATAAAAATAACGGGTTTTGaa             >  2:4753093/1‑150 (MQ=38)
                                                                                                                                          aCCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaa            >  2:6518240/1‑150 (MQ=255)
                                                                                                                                           ccccGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTGAAACCCCGCCTAAAAAGTTACTGTAGGATCGGGTTTTACCCCGAAGTGGCCCGCACACAAAAATAAAAAAAACAGATTTTAaaaa           >  1:2893529/1‑150 (MQ=255)
                                                                                                                                             ccGGCCACTGGGCCGGGC‑‑‑TTTTATTTTGGGGGTAT‑‑‑‑‑TTTTCCGGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa         <  2:8552955/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                             
GCAAGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑‑‑TTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCT  >  NC_002947/954934‑955231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: