Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,089,805 A→G 5.2% intergenic (+653/+265) 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,8050AG 5.2% 456.0 / 9.7 154intergenic (+653/+265)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (99/46);  new base G (3/5);  total (102/52)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 158 total aligned reads displayed.

TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2089656‑2089954
                                                                                                                                                     |                                                                                                                                                     
tGACGGGGATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                                                                                                       <  2:2761404/150‑1 (MQ=21)
tGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                                                                                                       <  2:1232915/150‑1 (MQ=34)
    gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaa                                                                                                                                                   >  1:3341389/1‑150 (MQ=34)
     ggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaaa                                                                                                                                                  <  2:4212037/150‑1 (MQ=34)
      gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                                 >  1:1577983/1‑150 (MQ=34)
      gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                                 >  1:4890680/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)
          ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                             >  1:3969557/1‑150 (MQ=34)
           gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAcc                                                                                                                                            >  1:2207869/1‑150 (MQ=34)
            aTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGATCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAccc                                                                                                                                           >  1:4687112/1‑150 (MQ=21)
            aTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGATCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAccc                                                                                                                                           >  1:4687029/1‑150 (MQ=21)
             tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATTATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccg                                                                                                                                          >  2:4151953/1‑150 (MQ=21)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  2:2308085/1‑150 (MQ=34)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  2:5093936/1‑150 (MQ=34)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTGTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  2:1657902/1‑150 (MQ=21)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                <  2:5422015/150‑1 (MQ=34)
                         gcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCgg                                                                                                                              >  2:2491672/1‑150 (MQ=34)
                            aggCCCCATCCTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaaa                                                                                                                           >  2:3615707/2‑150 (MQ=255)
                               ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGagatcg                                                                                                                        >  1:1566310/1‑145 (MQ=21)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  1:3423353/1‑150 (MQ=34)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  1:2936860/1‑150 (MQ=34)
                                 cctggtgggggCCTGCCGGGCGGTCCTCCGGGTAGCACGCCCGGTGTGTGGGGACCCGGGTGCGTTGGGTTGGGGCCGTATGTTGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                      <  1:4083199/143‑1 (MQ=11)
                                   atcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                                    >  2:195057/5‑150 (MQ=21)
                                     gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                              >  2:5294606/1‑122 (MQ=255)
                                           gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            >  2:2442750/1‑150 (MQ=34)
                                           gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            >  1:1357342/1‑150 (MQ=34)
                                            ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTa                                                                                                           >  2:3736839/1‑150 (MQ=35)
                                              tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                                                                                         >  2:4050922/1‑150 (MQ=37)
                                               gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt                                                                                                        >  2:3607444/1‑150 (MQ=37)
                                                  ggggggTCCTCCGGGTAGCACGCCGGGTGTGGGGGGCCCCGGGCGCGTTGGGTTGGGGCGGGATGTTGGGTGCCGAGGTTTTTGTTTACTTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGt                                                                                                     <  2:3410245/146‑1 (MQ=2)
                                                        tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatga                                                                                               >  1:697850/1‑150 (MQ=37)
                                                         cctccGGGTAGCACGCCGGTTGTGCGGGCACCCCGGAGCGTCGCGTAGGGGCCGGATGTTGGGAGCCGCCGGTTTTGGTTACTTGGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaatgaa                                                                                              <  2:5169852/150‑1 (MQ=11)
                                                          ctccGGGTAGCACGCCGGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGTATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             <  2:4415361/150‑1 (MQ=11)
                                                          ctccGGGGAGCAAGCCGGGTGTGAGGGGACCCGGGTGGTTTGCGTTGGGGCCGTATGTTGGGAGCCGCCGTTTTTGGTTACTTTGGCCACGACAAATGAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             <  2:3794250/150‑1 (MQ=11)
                                                          ctccCGGGTGCACGCCGGGCTTGTGGGGACCCCGGCGCGTCGCGTTGGGGCGGTTTTATGGGAGCCGAGGGTTTGGGTACCTTTGGCCACGACCAATGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             <  2:2577307/150‑1 (MQ=11)
                                                          ctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTgaatgaat                                                                                             >  1:2876466/1‑150 (MQ=37)
                                                             ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGtg                                                                                          >  2:4929504/1‑149 (MQ=37)
                                                             ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                          >  2:104535/1‑150 (MQ=37)
                                                              cGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                              <  2:2949926/97‑1 (MQ=11)
                                                                ggAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTAc                                                                                       <  2:4993618/150‑1 (MQ=37)
                                                                    cAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATTATGGGAGCCGACGGTTTTGGTTACTTTGGCAACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGa                                                                                   <  2:3576997/150‑1 (MQ=255)
                                                                     aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                                  >  2:5155842/1‑150 (MQ=37)
                                                                       gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc                                                                                >  2:608441/1‑150 (MQ=38)
                                                                           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)
                                                                                     gACCCCGGAGCGTCGCTTTGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcc                                                                  <  2:2166817/150‑1 (MQ=255)
                                                                                        cccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaa                                                               >  1:1101137/1‑150 (MQ=39)
                                                                                                  agcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaa                                                     >  1:5041580/1‑150 (MQ=255)
                                                                                                        ggggCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaa                                               >  1:5037398/1‑150 (MQ=255)
                                                                                                         gggcgttatgttgGGAGCGGATGTTTTTGGTTACTTTGGACAGTACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGTATAGGCGCCTACATCGTGTATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaa                                              <  2:3063058/143‑1 (MQ=14)
                                                                                                          ggcggtatgttgGGAGCCGAGGTTTTTGTTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaa                                             <  1:5441834/144‑1 (MQ=14)
                                                                                                           gCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                            >  2:1396611/1‑150 (MQ=255)
                                                                                                           gCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaag                                            >  1:3836332/1‑150 (MQ=255)
                                                                                                              ggATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagccc                                         >  1:3102790/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)
                                                                                                                  gatgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                     >  1:5166861/1‑150 (MQ=255)
                                                                                                                     gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcc                                  >  2:4016462/1‑150 (MQ=255)
                                                                                                                      ggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagccc                                 >  2:786877/1‑150 (MQ=255)
                                                                                                                          ccGTCGGTTTTGGTTACTTTGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                             <  2:4862121/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)
                                                                                                                                 ttttGGTTAGTTTGGCCACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      <  2:3969557/150‑1 (MQ=255)
                                                                                                                                 ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      >  2:3180708/1‑150 (MQ=255)
                                                                                                                                   ttGTTTACTTTGGCCACGACCAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGc                    <  1:1745052/150‑1 (MQ=255)
                                                                                                                                      gTTACTTTGGCCTCGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 <  2:2879390/150‑1 (MQ=255)
                                                                                                                                       tttcTTGGGCCACGACAAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  1:43461/147‑1 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  1:3182149/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  1:1693004/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  1:2973726/1‑150 (MQ=255)
                                                                                                                                       ttACTTGGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  2:4315469/150‑1 (MQ=255)
                                                                                                                                         aCTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgc              >  2:118417/1‑150 (MQ=255)
                                                                                                                                         aCTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgc              <  2:2207869/150‑1 (MQ=255)
                                                                                                                                          cTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcg             >  1:1037792/1‑150 (MQ=255)
                                                                                                                                          cTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcg             >  2:1629397/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  2:66934/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:826100/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:4622679/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:4641553/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  1:763157/1‑150 (MQ=255)
                                                                                                                                           tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGACCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcgc            >  2:3175265/1‑150 (MQ=255)
                                                                                                                                            ttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGcc           <  1:3586843/150‑1 (MQ=255)
                                                                                                                                             tGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagc                                           <  1:3359211/117‑1 (MQ=255)
                                                                                                                                             tGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagc                                           >  2:3359211/1‑117 (MQ=255)
                                                                                                                                                ccACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt       <  2:4687029/150‑1 (MQ=255)
                                                                                                                                                ccACGACCAACTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt       <  1:2657480/150‑1 (MQ=255)
                                                                                                                                                ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt       >  2:5238531/1‑150 (MQ=255)
                                                                                                                                                ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTATAGCGCCGGtt       >  1:3824923/1‑150 (MQ=255)
                                                                                                                                                ccACGACCAAAGCAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGtt       <  2:4687112/150‑1 (MQ=255)
                                                                                                                                                 cACGACAAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTc      <  2:1880449/150‑1 (MQ=255)
                                                                                                                                                     aCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAt  >  2:613130/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2089656‑2089954

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: