Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,599,059 A→C 5.9% intergenic (‑266/+35) PP_1399 ← / ← kgtP hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,599,0590AC 5.9% 893.5 / 30.4 265intergenic (‑266/+35)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (124/125);  new base C (8/8);  total (132/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCG  >  NC_002947/1598911‑1599208
                                                                                                                                                    |                                                                                                                                                     
gggCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCaa                                                                                                                                                      <  1:5624145/150‑1 (MQ=255)
gggCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCaa                                                                                                                                                      >  2:4766713/1‑150 (MQ=255)
  gcgcGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATg                                                                                                                                                    <  2:5535040/150‑1 (MQ=255)
    gcgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATggg                                                                                                                                                  <  1:1779508/150‑1 (MQ=255)
     cgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgggg                                                                                                                                                 <  1:2997266/150‑1 (MQ=255)
         tGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAGCGcc                                                                                                                                             <  2:2377283/150‑1 (MQ=255)
          gTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGTATAGAAAAGGGGCGCCCCATTGCATGGAGCGccc                                                                                                                                            <  2:4588199/150‑1 (MQ=255)
              tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                        <  2:4339814/150‑1 (MQ=255)
              tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                        <  1:3487916/150‑1 (MQ=255)
              tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                        <  1:3822981/150‑1 (MQ=255)
              tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                        <  2:4938517/150‑1 (MQ=255)
              tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAGCGCCCCttt                                                                                                                                        <  1:4645763/150‑1 (MQ=255)
               gCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCtttt                                                                                                                                       <  2:2264177/150‑1 (MQ=255)
                cTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                      <  2:3601381/150‑1 (MQ=255)
                cTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                      <  1:2592257/150‑1 (MQ=255)
                 ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAg                                                                                                                                     >  2:1074700/1‑150 (MQ=255)
                 ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAg                                                                                                                                     >  2:4469750/1‑150 (MQ=255)
                 ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAg                                                                                                                                     <  1:2084679/150‑1 (MQ=255)
                  tCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGt                                                                                                                                    >  2:726369/1‑150 (MQ=255)
                   cAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTa                                                                                                                                   <  2:2655338/150‑1 (MQ=255)
                           aaGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTCATAGGGCGTTAAGAAGGAAAAGTCTGACTAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccagg                                                                                                                           >  2:1052764/1‑150 (MQ=255)
                             ggctggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccaggcc                                                                                                                         >  1:4796467/1‑150 (MQ=255)
                                tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGc                                                                                                                      >  1:1032719/1‑150 (MQ=255)
                                 ggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCg                                                                                                                     >  1:3648128/1‑150 (MQ=255)
                                      gggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                                <  2:1208081/150‑1 (MQ=255)
                                      gggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                                <  2:1207966/150‑1 (MQ=255)
                                         aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAg                                                                                                             >  2:1595322/1‑150 (MQ=255)
                                         aGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAg                                                                                                             <  1:4114434/150‑1 (MQ=255)
                                          gCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGt                                                                                                            <  2:2156824/150‑1 (MQ=255)
                                          gCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAAGAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGt                                                                                                            <  2:1042312/150‑1 (MQ=255)
                                           cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                           <  2:4577996/150‑1 (MQ=255)
                                           cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                           >  1:819160/1‑150 (MQ=255)
                                           cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                           <  1:4331037/150‑1 (MQ=255)
                                           cGTTCGCGCGGGCGCCAACGTTCTGCCAGACATTGCATCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                           <  1:1113258/150‑1 (MQ=255)
                                                gcgcgGGCGCCAACGTTCTGCCAGCCATTGCAGCTGGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTCCTGCCAGGCCAGCGTTACTCAGTGGTAAt                                                                                                      >  2:3682903/1‑150 (MQ=255)
                                                    gggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgc                                                                                                  <  2:2464594/150‑1 (MQ=255)
                                                    gggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgc                                                                                                  >  1:202447/1‑150 (MQ=255)
                                                    gggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgc                                                                                                  <  2:1450406/150‑1 (MQ=255)
                                                     ggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgcg                                                                                                 <  2:484161/150‑1 (MQ=255)
                                                      gcgcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                                <  1:4183340/150‑1 (MQ=255)
                                                        gcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                              <  1:1243068/150‑1 (MQ=255)
                                                        gcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                              >  1:243620/1‑150 (MQ=255)
                                                           aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgct                                                                                           >  1:1515069/1‑150 (MQ=255)
                                                             cGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgctgc                                                                                         <  1:1802275/150‑1 (MQ=255)
                                                                tCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTg                                                                                      <  2:1709146/150‑1 (MQ=255)
                                                                  tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                    <  2:4460763/150‑1 (MQ=255)
                                                                  tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                    <  1:4125990/150‑1 (MQ=255)
                                                                  tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                    <  2:3559142/150‑1 (MQ=255)
                                                                   gccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtg                                                                                   >  1:9650/1‑150 (MQ=255)
                                                                    ccagccaTTGCAGCTTGATAAGGGGGAAGTAATAGGGCGTTAAGAAGGAAAAGGCTGAATAGAAAAGGGCCGCTCCCTGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCGAGgtgt                                                                                  >  1:2376214/1‑150 (MQ=255)
                                                                       gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCtt                                                                               >  2:1646324/1‑150 (MQ=255)
                                                                          atagCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTGTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCGTATTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                            <  1:457325/147‑1 (MQ=255)
                                                                          aTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                            >  2:2053387/1‑150 (MQ=255)
                                                                           ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCAt                                                                           >  1:408904/1‑150 (MQ=255)
                                                                             gCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATgg                                                                         >  1:2620202/1‑150 (MQ=255)
                                                                               aGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCt                                                                                <  2:4689230/141‑1 (MQ=255)
                                                                               aGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCt                                                                                >  1:4689230/1‑141 (MQ=255)
                                                                                 cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                     <  1:4477271/150‑1 (MQ=255)
                                                                                 cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                     <  1:5123549/150‑1 (MQ=255)
                                                                                 cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTAGGTGTCTTTCATGGTcgc                                                                     >  1:1757567/1‑150 (MQ=255)
                                                                                   tGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGt                                                                   <  1:720393/150‑1 (MQ=255)
                                                                                     aTAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTaa                                                                 <  1:2708253/150‑1 (MQ=255)
                                                                                          gggggaTTAATAGGGCGTTAAGAAGGAACAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTATACcag                                                            >  1:4345134/6‑150 (MQ=255)
                                                                                          gTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagc                                                           >  1:1420987/1‑150 (MQ=255)
                                                                                           tGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagc                                                           >  2:5149406/1‑150 (MQ=255)
                                                                                            ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                          <  2:1976646/150‑1 (MQ=255)
                                                                                            ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                          >  2:1124522/1‑150 (MQ=255)
                                                                                            ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                          >  2:3952315/1‑150 (MQ=255)
                                                                                             gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagcag                                                         >  2:3746306/1‑150 (MQ=255)
                                                                                              aaTTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCagcagc                                                        <  2:1680998/150‑1 (MQ=255)
                                                                                                 tAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                                >  1:2804282/1‑107 (MQ=255)
                                                                                                  aaTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGc                                                    <  2:2948988/150‑1 (MQ=255)
                                                                                                   aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCa                                                   >  1:2917035/1‑150 (MQ=255)
                                                                                                                aagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCg                                      <  1:1962757/150‑1 (MQ=255)
                                                                                                                   gAAAAGTCTGAATAGAAAAGGGGCACTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTg                                   >  1:568317/1‑150 (MQ=255)
                                                                                                                     aaaGTCTGAATAGAAAAGGGGCGCTCCAGGCACTGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGAc                                 <  1:5234628/150‑1 (MQ=255)
                                                                                                                         tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                             <  1:5514162/150‑1 (MQ=255)
                                                                                                                         tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                             >  1:3594289/1‑150 (MQ=255)
                                                                                                                           tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAg                                                  <  1:5012309/127‑1 (MQ=255)
                                                                                                                             aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGt                         <  2:5210508/150‑1 (MQ=255)
                                                                                                                                  aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                    >  2:2764327/1‑150 (MQ=255)
                                                                                                                                  aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                    >  2:1347313/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                   >  2:3464366/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                   >  2:3426647/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGCGCCCCATTGCATGGAGCACCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                   >  2:3165165/1‑150 (MQ=255)
                                                                                                                                    aaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGcc                  >  2:5011648/1‑150 (MQ=255)
                                                                                                                                     aGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCg                 <  1:2464792/150‑1 (MQ=255)
                                                                                                                                      ggggCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGc                >  1:3395123/1‑150 (MQ=255)
                                                                                                                                       gggCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct               <  2:5400925/150‑1 (MQ=255)
                                                                                                                                       gggCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct               >  1:1871209/1‑150 (MQ=255)
                                                                                                                                        ggCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctc              >  1:1026558/1‑150 (MQ=255)
                                                                                                                                         gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct             >  2:4089421/1‑150 (MQ=255)
                                                                                                                                         gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct             >  2:2563433/1‑150 (MQ=255)
                                                                                                                                         gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct             >  1:2233115/1‑150 (MQ=255)
                                                                                                                                             tCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAt         >  2:5078441/1‑150 (MQ=255)
                                                                                                                                              ccATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATc        <  2:3150617/150‑1 (MQ=255)
                                                                                                                                               cATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATcc       <  2:2323472/150‑1 (MQ=255)
                                                                                                                                                aTGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATccc      >  1:908321/1‑150 (MQ=255)
                                                                                                                                                   cAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGgc   >  2:4267904/1‑150 (MQ=255)
                                                                                                                                                    aaTGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGgcg  >  1:135700/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCG  >  NC_002947/1598911‑1599208

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