Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,089,800 C→T 5.3% intergenic (+648/+270) 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,089,8000CT 5.3% 432.1 / 10.0 150intergenic (+648/+270)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base C (95/46);  new base T (3/5);  total (98/52)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 155 total aligned reads displayed.

GGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTT  >  NC_002947/2089654‑2089949
                                                                                                                                                  |                                                                                                                                                     
ggTGACGGGCTTGGATGCCCGTCAAGCGGTGGGCCCCCGGATGGGGCCTGCAGCGGGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCAc                                                                                                                                                    <  2:618121/150‑1 (MQ=255)
ggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCAc                                                                                                                                                    <  2:3592929/150‑1 (MQ=34)
  tGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                                                                                                  <  2:1232915/150‑1 (MQ=34)
      gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaa                                                                                                                                              >  2:5058583/1‑150 (MQ=34)
       ggCATGGATGCCCGTCATGCGGTGGGCCCCTGGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGCGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             <  1:4813729/150‑1 (MQ=11)
       ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                             <  2:4212037/150‑1 (MQ=34)
        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGCGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                            >  2:3990666/1‑150 (MQ=21)
        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                            >  1:1577983/1‑150 (MQ=34)
        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTATGGTTACTTTGGCCACGACCAAAg                                                                                                                                            >  1:1577841/1‑150 (MQ=21)
           tggtggCCCGTCACGCGGTGGGCCCCTGTTTGGGGCCGGCCGCGCGGTCCTCCGGGGAGCAAGCCGGGCGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  1:601052/145‑1 (MQ=11)
           tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  2:1339486/150‑1 (MQ=255)
           tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                         <  2:3507558/150‑1 (MQ=34)
             gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAcc                                                                                                                                       >  2:3684928/1‑150 (MQ=34)
              aTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGATCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAccc                                                                                                                                      >  1:4687112/1‑150 (MQ=21)
                 cccGACAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCACGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTGGGTTACTTTGGTCGTGGCCAAAGTAAccaaaa                                                                                                                                   <  1:739271/150‑5 (MQ=11)
                  ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                  >  2:2308106/1‑150 (MQ=34)
                  ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                  >  2:2308085/1‑150 (MQ=34)
                  ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTGTGGCCACGACCAAAGTAACccgccg                                                                                                                                  >  2:1657902/1‑150 (MQ=21)
                           tgcaggCCCCATCCTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                         >  1:3149589/5‑150 (MQ=255)
                           gcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                         >  1:2117613/1‑150 (MQ=34)
                           gcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                         >  2:2491672/1‑150 (MQ=34)
                            cTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGa                                                                                                                        >  2:2153101/1‑150 (MQ=34)
                              gggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaaa                                                                                                                      <  2:4072091/150‑1 (MQ=34)
                              aggCCCCATCCTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaaa                                                                                                                      >  2:3615707/2‑150 (MQ=255)
                               ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAg                                                                                                                     >  1:3945257/1‑150 (MQ=34)
                                 ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGagatcg                                                                                                                   >  1:1566310/1‑145 (MQ=21)
                                  cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                  >  1:2936860/1‑150 (MQ=34)
                                  cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                  >  2:4977356/1‑150 (MQ=34)
                                   cctggtgggggCCTGCCGGGCGGTCCTCCGGGTAGCACGCCCGGTGTGTGGGGACCCGGGTGCGTTGGGTTGGGGCCGTATGTTGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 <  1:4083199/143‑1 (MQ=11)
                                   cccggtgggggCCTGCAGCGGGGCCCCCCGGGGAGCATGCGCGGTGTGTGGGGACCCCGGTGCGTAGCGTCGGGGCGGGATGTTGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 <  1:4083179/143‑1 (MQ=11)
                                   ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  1:3410245/1‑150 (MQ=34)
                                     atcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                               >  2:194249/5‑150 (MQ=21)
                                       gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         >  2:5294606/1‑122 (MQ=255)
                                       gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  1:5294606/122‑1 (MQ=255)
                                           gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                         >  2:5191988/1‑150 (MQ=34)
                                             gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                       >  1:1357342/1‑150 (MQ=34)
                                                    ggggggTCCTCCGGGTAGCACGCCGGGTGTGGGGGGCCCCGGGCGCGTTGGGTTGGGGCGGGATGTTGGGTGCCGAGGTTTTTGTTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGt                                                                                                <  2:3410245/146‑1 (MQ=2)
                                                           cctccGGGTAGCACGCCGGTTGTGCGGGCACCCCGGAGCGTCGCGTAGGGGCCGGATGTTGGGAGCCGCCGGTTTTGGTTACTTGGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatgaa                                                                                         <  2:5169852/150‑1 (MQ=11)
                                                            ctccGGGGAGCAAGCCGGGTGTGAGGGGACCCGGGTGGTTTGCGTTGGGGCCGTATGTTGGGAGCCGCCGTTTTTGGTTACTTTGGCCACGACAAATGAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        <  2:3794250/150‑1 (MQ=11)
                                                            ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                        >  2:2448162/1‑150 (MQ=37)
                                                             tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGaatgaatg                                                                                       >  2:3586843/1‑150 (MQ=37)
                                                             tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGaatgaatg                                                                                       >  1:4578030/1‑150 (MQ=37)
                                                               ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGtg                                                                                     >  2:4929504/1‑149 (MQ=37)
                                                               ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:104535/1‑150 (MQ=37)
                                                               cagggAGCAAGCCCGGAGTGAGGGGACACCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                     >  2:3360384/3‑150 (MQ=37)
                                                                cGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  2:2949926/97‑1 (MQ=11)
                                                                cGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         >  1:2949926/1‑97 (MQ=11)
                                                                       aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                             >  2:3803079/1‑150 (MQ=37)
                                                                           ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGa                                                                         >  1:4658969/1‑150 (MQ=39)
                                                                             ggtgtgtGGGGACCCGGGTGGGTTGGTTTGGGGCGGGATGTTGGGAGCCGCGGTTTTTGGTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                       <  2:4578030/147‑1 (MQ=14)
                                                                                 ggtggggACCCCGGTGGTTTGGGTTGGGGCGTTATTTTGGGAGCCTAGGTTTTTGTTTACTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGGAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAgcgc                                                                   <  2:603574/147‑1 (MQ=11)
                                                                                      ngACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAc                                                              >  1:2000941/2‑150 (MQ=39)
                                                                                       gACCCCGGAGCGTTGCGTTGGGGCCGGTTGATGGGAGCCGACGGTTTTGGTTACTTGGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcc                                                             <  2:2166928/150‑1 (MQ=255)
                                                                                          cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                          >  2:301469/1‑150 (MQ=39)
                                                                                          cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                          >  1:1101137/1‑150 (MQ=39)
                                                                                                 cgtagcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagccca                                                   <  1:4977356/150‑1 (MQ=255)
                                                                                                          ggggCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaa                                          >  1:5037398/1‑150 (MQ=255)
                                                                                                            ggcggtatgttgGGAGCCGAGGTTTTTGTTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                        <  1:5441834/144‑1 (MQ=14)
                                                                                                             gCCGGTTTATGGTAGCCGTCGGTTTTGGTTACTTTGGCCACGAAAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                       <  2:2856014/150‑1 (MQ=14)
                                                                                                             gCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                       >  2:1396611/1‑150 (MQ=255)
                                                                                                             gCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                       >  1:3836332/1‑150 (MQ=255)
                                                                                                                ggATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccc                                    >  1:3102790/1‑150 (MQ=255)
                                                                                                                ggATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccc                                    >  1:3745722/1‑150 (MQ=255)
                                                                                                                 gatgatGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccca                                   >  1:3190804/1‑150 (MQ=255)
                                                                                                                 gatgatGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccca                                   >  1:5090896/1‑150 (MQ=255)
                                                                                                                    gttgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                <  2:5090896/148‑1 (MQ=255)
                                                                                                                    gatgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                <  2:2117613/150‑1 (MQ=255)
                                                                                                                       gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcc                             >  2:4016462/1‑150 (MQ=255)
                                                                                                                        ggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagccc                            >  2:786877/1‑150 (MQ=255)
                                                                                                                            ccGACGGTTTTGTTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                        <  2:5124272/150‑1 (MQ=255)
                                                                                                                            ccGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                        >  1:2468772/1‑150 (MQ=255)
                                                                                                                             cGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                       >  2:3109956/1‑150 (MQ=255)
                                                                                                                              gACGGTTTTGGTTACTTTGGCCACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGTATGTCTACTTGATGCGAAAGCGCACAAGccaaaagcccaaaagcccaatagcccaaaagcccaaaagc                      >  1:2440573/1‑150 (MQ=255)
                                                                                                                              gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                      >  1:2215073/1‑150 (MQ=255)
                                                                                                                               aggtttttGGTTACTTGGGCAACGACCAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCcaaaagcccaaaagcccaaaagcccaaaagcccaaaagcc                     <  2:2827973/146‑1 (MQ=255)
                                                                                                                                  gTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAaagcccaaaagcccaaaagcccaaaagcccaaaagcccaa                  >  2:4426531/1‑150 (MQ=255)
                                                                                                                                   ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                 >  2:3180708/1‑150 (MQ=255)
                                                                                                                                   ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                 >  2:2168963/1‑150 (MQ=255)
                                                                                                                                        gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt            >  2:3089850/1‑150 (MQ=255)
                                                                                                                                         ttAGTTTGGCCACGACAAATGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           <  1:5058583/150‑1 (MQ=255)
                                                                                                                                         ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           >  2:5310790/1‑150 (MQ=255)
                                                                                                                                         ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           >  1:1693004/1‑150 (MQ=255)
                                                                                                                                         ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           >  1:1693392/1‑150 (MQ=255)
                                                                                                                                         ttACTTGGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           <  2:4315469/150‑1 (MQ=255)
                                                                                                                                         taaattgggCCACCACAAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAAGGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa           <  1:912114/143‑1 (MQ=18)
                                                                                                                                          tACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAg          <  2:2607750/150‑1 (MQ=255)
                                                                                                                                           aCTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgc         >  2:118417/1‑150 (MQ=255)
                                                                                                                                           aCTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgc         <  1:2707227/150‑1 (MQ=255)
                                                                                                                                            cTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcg        >  1:1037792/1‑150 (MQ=255)
                                                                                                                                            cTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcg        >  2:1629397/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc       >  1:763157/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc       >  2:66934/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc       >  2:3175375/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc       >  1:4622679/1‑150 (MQ=255)
                                                                                                                                              ttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGcc      <  1:3586843/150‑1 (MQ=255)
                                                                                                                                               tGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagc                                      <  1:3359211/117‑1 (MQ=255)
                                                                                                                                                 gCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGt   >  2:561470/1‑150 (MQ=255)
                                                                                                                                                  ccACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt  <  2:4687029/150‑1 (MQ=255)
                                                                                                                                                  ccACGACCAACTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt  <  1:2657480/150‑1 (MQ=255)
                                                                                                                                                  ccACGACCAAAGCAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt  <  2:4687112/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTT  >  NC_002947/2089654‑2089949

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: