Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,760 A→T 6.1% intergenic (+581/+337) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,092,7600AT 6.1% 760.8 / 5.5 277intergenic (+581/+337)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (186/64);  new base T (2/20);  total (191/86)
Fisher's exact test for biased strand distribution p-value = 1.44e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGC  >  NC_002947/2092611‑2092909
                                                                                                                                                      |                                                                                                                                                     
cttaGCGCGCGGTAGTTCAGGCGACGCATATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGGATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                                       <  1:8464052/147‑1 (MQ=12)
cTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                                       <  2:5830079/150‑1 (MQ=32)
cTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                                       <  2:724458/150‑1 (MQ=17)
aacgctggctcctatcatccggcccctgcgctacgctccggggttccctcgcgcaGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                                       >  2:3776693‑M1/56‑150 (MQ=255)
  aaGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGt                                                                                                                                                     >  1:7802691/1‑150 (MQ=32)
    tcggctcccctcatccggcccctacgctacgctccggggtcccctcactccGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGa                                                                                                                                                   >  1:6324256‑M1/52‑150 (MQ=255)
    gcgcgcGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGa                                                                                                                                                   >  2:771869/1‑150 (MQ=32)
     cggctcccatcatccggcccctacgctacgctccggggtcccctcactcaGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAg                                                                                                                                                  >  1:5193349‑M1/51‑150 (MQ=255)
      gcgcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgg                                                                                                                                                 >  1:1132434/1‑150 (MQ=32)
      gcgcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgg                                                                                                                                                 <  1:1024947/150‑1 (MQ=32)
      gcgcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgg                                                                                                                                                 >  1:8527943/1‑150 (MQ=32)
      gcgcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgg                                                                                                                                                 <  1:2083000/150‑1 (MQ=32)
        gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAgggg                                                                                                                                               >  1:6023907/1‑150 (MQ=32)
           gTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAcc                                                                                                                                            >  2:2409401/1‑150 (MQ=32)
            tAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAccc                                                                                                                                           >  2:415782/1‑150 (MQ=32)
             aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGTTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAcccc                                                                                                                                          <  1:4918013/150‑1 (MQ=17)
                tCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGa                                                                                                                                       <  2:5231198/150‑1 (MQ=32)
                tCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCACGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGa                                                                                                                                       <  2:659698/150‑1 (MQ=17)
                     cGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgta                                                                                                                                  <  2:184721/150‑1 (MQ=32)
                       aCGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCTGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtagc                                                                                                                                <  2:6122377/150‑1 (MQ=14)
                       aCGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtagc                                                                                                                                <  2:4399299/150‑1 (MQ=31)
                                       cAGGAGGCCGAGCGCTGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGgatgat                                                                                                                <  2:1071986/150‑1 (MQ=11)
                                        cggcggCCGCGCGCTGGGCTTGCGGAGGGCGGTGAGGGGCATGGATGCCCGTCACGCGCGGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGatgatg                                                                                                               <  1:2593499/146‑1 (MQ=11)
                                           tggCCGAGCGCTGGGCTTGCGGAGGGAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGTAACAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa                                                                                                            <  2:1184114/149‑1 (MQ=255)
                                                gAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGa                                                                                                       >  1:8337218/1‑150 (MQ=32)
                                                 tcgctccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcgtgacgggcatccatgcccgtcacctcccgtcgcaagcccGGCGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGagcccg                                                                                                       >  1:8409231‑M2/100‑148 (MQ=255)
                                                  gcgcTGGGCTTGCGGCGGGAGGGGAGGGGCATGGATGCCCGACACGCGCTGGGCCCCAGGCTGGGGCCTGCAGCGCGGTCCTCCGGGGAGCAAGCCCGGAGTGCGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACg                                                                                                     <  1:3364233/150‑1 (MQ=11)
                                                   cgctGGGCTTGCGGCGGGAGGTGAGGGGCATGGAGGCCGGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGCCCTCCGGGGAGCAAGCCGGGTGTGCGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACgg                                                                                                    >  1:5635360‑M1/5‑150 (MQ=255)
                                                    gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGt                                                                                                   >  1:772769/1‑150 (MQ=32)
                                                     tccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccTGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGAcggc                                                                                                   >  1:7989265‑M2/96‑149 (MQ=255)
                                                     tccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGt                                                                                                   >  1:812133‑M2/96‑150 (MQ=255)
                                                      aGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGttt                                                                                                 <  2:2717909/150‑1 (MQ=34)
                                                      aGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGttt                                                                                                 >  1:7292176/1‑150 (MQ=34)
                                                       cgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGttt                                                                                                 >  2:5414456‑M2/94‑150 (MQ=255)
                                                       cgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGctt                                                                                                 >  2:5062010‑M2/94‑150 (MQ=255)
                                                       GGGCTTGCGGTGGGAGGTGACGGGCATGGATGCCCGTCATGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                                                >  2:3967314‑M1/1‑150 (MQ=255)
                                                        gggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                                                >  2:3177855‑M2/93‑150 (MQ=255)
                                                         ggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagcccTGCGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGgttttt                                                                                               >  1:2296772‑M2/92‑149 (MQ=255)
                                                             GCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTAc                                                                                          >  2:3844566‑M1/1‑150 (MQ=255)
                                                             GCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTAc                                                                                          <  1:180928‑M1/150‑1 (MQ=255)
                                                              CGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTGTTGGTTACt                                                                                         <  2:4886396‑M1/150‑1 (MQ=255)
                                                                 tgggaggTGACGGGCTTGGATGCCCGTCAAGCGCTGGGCCCCTGGTTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCTTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTg                                                                                      <  1:6539779/149‑1 (MQ=255)
                                                                 cgggtggGGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGGGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTg                                                                                      <  1:5315675/149‑1 (MQ=255)
                                                                      ggTGAGGGGCATGGATGCCCGTCAAGCGCTGGGCCCCTGGATGGGGCCTGCAGCGCGGTCCTCCCGGTAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCAc                                                                                 <  1:7843656/150‑1 (MQ=255)
                                                                       ttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACg                                                                                >  2:7093092/2‑150 (MQ=21)
                                                                       gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACg                                                                                >  1:1316770/1‑150 (MQ=34)
                                                                        tGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGTTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                               <  2:29410/150‑1 (MQ=21)
                                                                        tGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                               >  2:3721294/1‑150 (MQ=21)
                                                                         gACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGAc                                                                              >  2:1497644/1‑150 (MQ=21)
                                                                            gggCATGGTTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaa                                                                           <  2:3947131/150‑1 (MQ=21)
                                                                            gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTGTTGGTTACTTTGGCCACGACCaa                                                                           >  2:4658540/1‑150 (MQ=21)
                                                                            gggCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaa                                                                           >  1:4350218/1‑150 (MQ=21)
                                                                              gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAg                                                                         <  2:6902482/150‑1 (MQ=255)
                                                                              gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                         >  1:1465539/1‑150 (MQ=34)
                                                                              gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                         >  2:5268407/1‑150 (MQ=34)
                                                                                aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTa                                                                       >  1:4029218/1‑150 (MQ=34)
                                                                                 tGGATGCCCGTCACGCGCGGGGCCCCTGGTTGGGGCCTGCCGCGCGGTCCTCCCGGGAGCAAGCCGGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                      <  1:3038372/150‑1 (MQ=11)
                                                                                 tGGAGGCCCGTCAAGCGCTGGGCCCCTGGATGGGGCCTGCTGCGCGGTCCTCCCGGGTGCAAGCCGGGTGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                      <  1:2526567/150‑1 (MQ=11)
                                                                                  ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                     >  2:2171123/1‑150 (MQ=34)
                                                                                  ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                     >  2:625674/1‑150 (MQ=34)
                                                                                        ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                               >  1:5499971/1‑150 (MQ=34)
                                                                                         cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCgg                                                              >  1:8716507/1‑150 (MQ=34)
                                                                                          gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGa                                                             >  2:5123772/1‑150 (MQ=34)
                                                                                             aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAgg                                                          >  2:7846581/1‑150 (MQ=34)
                                                                                               gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                        >  2:2103846/1‑150 (MQ=34)
                                                                                               gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                        >  2:4142100/1‑150 (MQ=34)
                                                                                                cgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCg                                                       >  2:4315682/1‑150 (MQ=34)
                                                                                                  tgaggCCCCATGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                      >  2:6744098/4‑150 (MQ=255)
                                                                                                   tGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaa                                                    >  1:6409957/1‑150 (MQ=34)
                                                                                                    gggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaaa                                                   >  2:1875384/1‑150 (MQ=34)
                                                                                                        cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                               >  1:4872567/1‑150 (MQ=34)
                                                                                                        cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                               >  2:3508041/1‑150 (MQ=34)
                                                                                                         ccAGGTTGGGGCCTGCAGCGCGGTCCTCCCGGTAGCAAGCCCGGAGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                              <  2:8813945/150‑1 (MQ=255)
                                                                                                         ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                              >  2:6470550/1‑150 (MQ=34)
                                                                                                         ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                              >  1:5842873/1‑150 (MQ=34)
                                                                                                          cAGGCTGGGGCCTGCAGCGGGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCTACGGTTTTGGTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGaa                                             <  2:4947476/150‑1 (MQ=11)
                                                                                                          cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa                                             >  2:8369213/1‑150 (MQ=34)
                                                                                                           ctgtctGGGGCCTGCCGGGGGTTCCTCCGGGTTGAACGCGGGGTGTGTGGGGACCCGGGTGGGTCGCGTTGGGGCGGTTTTTTGGGCGCGGCGGTTTTTGGTTACTTGGGCCCCGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGaa                                             <  2:1672363/145‑1 (MQ=2)
                                                                                                           aGGTTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGTGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                            <  2:5807852/150‑1 (MQ=21)
                                                                                                               tccgggCCTGCAGCGCTGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAgg                                         >  1:2790438/4‑150 (MQ=21)
                                                                                                                tGGGGCCTGCTGCGGGGTCCTCCGGGGAGCACGCCGGGTGTGCGGGTACCCGGGTGCTTTGGGTTGGGGGGGTATTTTGGTAGCCTAGGGTTTTGTTTAGTTGGGCCACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgc                                        <  2:8469100/150‑1 (MQ=2)
                                                                                                                   ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                               >  1:7917807/1‑140 (MQ=34)
                                                                                                                    gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                    >  2:1547461/1‑150 (MQ=34)
                                                                                                                    gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                    >  2:8110572/1‑150 (MQ=34)
                                                                                                                    gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                    >  1:4690169/1‑150 (MQ=34)
                                                                                                                      ctgctgCGGGGTCCTCCGGTTACCACGCCCGTTGTGCGGGTCCCCGGGCGCGTAGGGTTGGGGCCGTATTATGGGTGCCGCCGGTTTTTGTTACTTTGGCAACGACAAATTAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTAc                                  <  2:8199100/150‑1 (MQ=2)
                                                                                                                       tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                 >  2:44622/1‑150 (MQ=37)
                                                                                                                       tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                 >  2:8577882/1‑150 (MQ=37)
                                                                                                                        gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt                                >  1:8119341/1‑150 (MQ=37)
                                                                                                                          aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCg                              >  1:6952180/1‑150 (MQ=37)
                                                                                                                               ggTCCGCCGGGGAGCACGCCGGTTGTGTGGGGACCGCGGTGGGTCGGGTTGGGGCCGTTTGTTGGGAGCGGAGGTTTTTGGTTACTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaat                         <  2:4096730/150‑1 (MQ=2)
                                                                                                                                   ctccCGGTAGCACGCCGGGTGTGTGGGGACCCCGTAGCGTAGCGTCGGGGCGGTTTGATGGGAGCCTACGTTTTTGTTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                     <  2:8405191/150‑1 (MQ=11)
                                                                                                                                      ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGACTGAATGTc                  >  2:8077697/1‑150 (MQ=37)
                                                                                                                                      ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                  >  2:2355254/1‑150 (MQ=37)
                                                                                                                                              aagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt          >  2:7267956‑M2/7‑150 (MQ=255)
                                                                                                                                                gcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc        >  2:6277823‑M2/5‑150 (MQ=255)
                                                                                                                                                gcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc        >  1:3555805‑M2/5‑150 (MQ=255)
                                                                                                                                                 cccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTAGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCg       >  1:7116621‑M2/4‑150 (MQ=255)
                                                                                                                                                  ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGa      >  1:6623026‑M2/3‑150 (MQ=255)
                                                                                                                                                    ggAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa    <  2:6507957/150‑1 (MQ=39)
                                                                                                                                                      |                                                                                                                                                     
CTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGC  >  NC_002947/2092611‑2092909

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: