Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,832 A→G 6.9% 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,092,8320AG 6.9% 793.0 / 35.7 289intergenic (+653/+265)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (181/86);  new base G (14/6);  total (195/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 303 total aligned reads displayed.

TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2092683‑2092981
                                                                                                                                                     |                                                                                                                                                     
tGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGTTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGa                                                                                                                                                       <  2:29410/150‑1 (MQ=21)
    gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCaa                                                                                                                                                   >  1:5699358/1‑150 (MQ=34)
    gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTGTTGGTTACTTTGGCCACGACCaa                                                                                                                                                   >  2:4658540/1‑150 (MQ=21)
      gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGATCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGTCGTGGCCAAAg                                                                                                                                                 <  1:1421603/150‑1 (MQ=11)
      gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                                 >  1:7120391/1‑150 (MQ=34)
         tGGATGCCCGTCACGCGCGGGGCCCCTGGTTGGGGCCTGCCGCGCGGTCCTCCCGGGAGCAAGCCGGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                              <  1:3038372/150‑1 (MQ=11)
         tGGAGGCCCGTCAAGCGCTGGGCCCCTGGATGGGGCCTGCTGCGCGGTCCTCCCGGGTGCAAGCCGGGTGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                              <  1:2526567/150‑1 (MQ=11)
          ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAAc                                                                                                                                             >  2:2170711/1‑150 (MQ=34)
              gCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgc                                                                                                                                         >  1:2612120/1‑150 (MQ=34)
               cccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgcc                                                                                                                                        >  2:3783255/1‑150 (MQ=34)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  2:1840633/1‑150 (MQ=34)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  1:5499282/1‑150 (MQ=34)
                ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                       >  2:6262536/1‑150 (MQ=34)
                 cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCgg                                                                                                                                      >  1:8716507/1‑150 (MQ=34)
                  gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGa                                                                                                                                     >  2:5123772/1‑150 (MQ=34)
                   tCACGCGCTGGGCCCCTGGATGGGGCCTGCAGGGGGGTCCTCCGGGGAGCACGCCCGGTGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGTTGGTAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGaa                                                                                                                                    <  2:4241246/150‑1 (MQ=11)
                     aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAgg                                                                                                                                  >  2:7846581/1‑150 (MQ=34)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                >  1:2360568/1‑150 (MQ=34)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                >  2:4142099/1‑150 (MQ=34)
                       gcgcTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGc                                                                                                                                >  2:2103846/1‑150 (MQ=34)
                            gggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGaaa                                                                                                                           >  2:1875384/1‑150 (MQ=34)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  1:4872824/1‑150 (MQ=34)
                                cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTg                                                                                                                       >  2:8423752/1‑150 (MQ=34)
                                 ccAGGATGGGGCCTGCAGCGCGGTCCTCCGGGTAGCATGCCTGGAGTGAGGGTACCCGGGAGGTTAGGGTAGGGGCGGTTTGATGGTAGCCGACGGTTTTGTTTACTTTGGCCAGGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                      <  2:5726522/150‑1 (MQ=11)
                                   cggtgggggCCTGCCGGGGGGGCCGCCCGGGAGCATGCCGGTAGGGCGGGGACCCGGGAGGGTAGCGTAGGGGCCGGATTCTGGGAGCCGCGGGTTTTGGTTACTTGGGCCACGACCAAATTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                                    <  2:4890494/145‑1 (MQ=2)
                                   aGGTTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGTGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                                    <  2:5807852/150‑1 (MQ=21)
                                      tccgggCCTGCAGCGCTGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAgg                                                                                                                 >  1:2790438/4‑150 (MQ=21)
                                      tccgggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAgg                                                                                                                 >  1:6427359/4‑150 (MQ=21)
                                       tGGGGCCTGCTGCGGGGTCCTCCGGGGAGCACGCCGGGTGTGCGGGTACCCGGGTGCTTTGGGTTGGGGGGGTATTTTGGTAGCCTAGGGTTTTGTTTAGTTGGGCCACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgc                                                                                                                <  2:8469100/150‑1 (MQ=2)
                                       gggggCCTGCAGCGCGGGCCTCCGGGGAGCACGCCGGGTGTGTGGGGACCCGGGAGCGTCGGGTCGGGGCGGGATGATGGGAGCCGCCGGTTTTGGTTACTTTGGCCACGACCAAAGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGgc                                                                                                                <  2:6739918/149‑1 (MQ=11)
                                           gCCTGCTGCGGGGTCCGCCGGGTTGCAAGCCGGGTGTGTGGGGACCCCGGTGGGTTGCGTTGGGGCGGGTTGTTGGGAGCCGAGGGTTTTGGTTAGTTTGGCCACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            <  2:5244573/150‑1 (MQ=11)
                                           gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            >  1:4690169/1‑150 (MQ=34)
                                           gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            >  2:8110572/1‑150 (MQ=34)
                                           gCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                                                                                                            >  2:1547461/1‑150 (MQ=34)
                                              tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                                                                                         >  2:44622/1‑150 (MQ=37)
                                              tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACa                                                                                                         >  2:8577882/1‑150 (MQ=37)
                                               gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACAt                                                                                                        >  1:8118937/1‑150 (MQ=37)
                                                 aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCg                                                                                                      >  1:401286/1‑150 (MQ=37)
                                                    gcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGACGGGCGGAAAGGTGAATAGGCGCCTACATCGtga                                                                                                   >  1:5018277/1‑150 (MQ=37)
                                                     cgGTCCTCCGGGGAGCAAGCCCGGTGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaa                                                                                                  <  2:5018277/150‑1 (MQ=255)
                                                     cgGTCCTCCCGGGAGCAAGCCGGGAGTGGGGGGACCCGGGCGCGTAGGGTTGGGGCGGGTTGAGGGGAGCCGAGGGTTTTGGTTACTTGGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGtgaa                                                                                                  <  2:5376247/150‑1 (MQ=11)
                                                     cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATCGGCGCCTACATCGtgaa                                                                                                  >  2:8502822/1‑150 (MQ=37)
                                                           tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGaatgaatg                                                                                            >  1:1110761/1‑150 (MQ=37)
                                                           tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGaatgaatg                                                                                            >  2:1404055/1‑150 (MQ=37)
                                                             ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                          >  2:8392976/1‑150 (MQ=37)
                                                             ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTc                                                                                          >  2:2355254/1‑150 (MQ=37)
                                                                  aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACtt                                                                                     >  2:7649545/1‑150 (MQ=37)
                                                                    catgcccGGTGTGTGGGGACCCGGGTGCGTTGCTTTGGGGCCGGATGTTGGGAGCCGAGGGTTTTGGTTACTTGGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGa                                                                                   <  2:7008686‑M2/143‑1 (MQ=255)
                                                                     aagcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGAt                                                                                  >  2:3867072‑M2/7‑150 (MQ=255)
                                                                      agcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATg                                                                                 >  1:8569683‑M2/6‑150 (MQ=255)
                                                                       gcccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGc                                                                                >  1:3555805‑M2/5‑150 (MQ=255)
                                                                         ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGa                                                                              >  1:6623026‑M2/3‑150 (MQ=255)
                                                                         ccGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGa                                                                              >  1:3257759‑M2/3‑150 (MQ=255)
                                                                           ggtgtgtGGGGACCCGGGCGGGTTGCTTGGGGGCCGTTTGTTGGGAGCCGAGGTTTTTGTTTACTTTGGCAACGACCAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                            <  2:2498574/147‑1 (MQ=14)
                                                                           ggAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGaaa                                                                            <  2:6507957/150‑1 (MQ=39)
                                                                                         ccGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaa                                                              >  1:7866865/1‑150 (MQ=39)
                                                                                          cGGGGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                             >  1:2496685/1‑150 (MQ=255)
                                                                                          cGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaa                                                             >  2:7477726/1‑150 (MQ=40)
                                                                                                gtagcgtagGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaa                                                       >  1:925806/1‑150 (MQ=255)
                                                                                                  tgcgttgGGGCGGGATGTTGGGAGCCGACGTTTTGGGTTACTTTGGCCACGACCAATTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaa                                                     <  1:6480118/149‑1 (MQ=14)
                                                                                                     gtagGGGCCGTATTTTGGTAGCCGACGTTTTTGTTTAATTTGGCCACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcc                                                  <  2:8337218/150‑1 (MQ=14)
                                                                                                         gggCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaa                                              >  2:5965478/1‑150 (MQ=255)
                                                                                                                  tttgGGTGCCGACGTTTTTGGTTACTTTGCCCACGACAAATTAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                     <  2:3982101/148‑1 (MQ=14)
                                                                                                                  gatgGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaa                                     >  1:705780/1‑150 (MQ=255)
                                                                                                                     gggAGCCTAGTGTTTTGGTTACTTGGGCAACTACAAATTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcc                                  <  2:230381/150‑1 (MQ=17)
                                                                                                                     gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcc                                  >  2:1734503/1‑150 (MQ=255)
                                                                                                                     gggAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcc                                  >  2:6378839/1‑150 (MQ=255)
                                                                                                                         gCCGAGGGTTTGGTTTACTTTGGCAACGACAAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATCGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                              <  2:773768/150‑1 (MQ=17)
                                                                                                                         gCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                              <  2:8210960/150‑1 (MQ=255)
                                                                                                                         gCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                              >  1:194942/1‑150 (MQ=255)
                                                                                                                          ccgtgggTTTTGGTTACTTTGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaa                             <  2:3521389/145‑1 (MQ=255)
                                                                                                                           ctaggtttttGTTTATTTTGGCAACCACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                            <  1:6900034/144‑1 (MQ=17)
                                                                                                                           cGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAAcccaaaagcccaaaagcccaaaagcccaaaagcccaaaag                            >  1:632848/1‑150 (MQ=255)
                                                                                                                            gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                           >  1:5330703/1‑150 (MQ=255)
                                                                                                                            gACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACccaaaagcccaaaagcccaaaagcccaaaagcccaaaagc                           >  2:4451690/1‑150 (MQ=255)
                                                                                                                               gtttttGTTTTCTTTGCCCCCTACAAATTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAaaagcccaaaagcccaaaagcccaaaagcccaaaagccca                        <  1:8138097/148‑1 (MQ=17)
                                                                                                                               ggTTTTCGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAaaagcccaaaagcccaaaagcccaaaagcccaaaagccca                        >  2:3328102/1‑150 (MQ=255)
                                                                                                                                 ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGACGGGCGGAAAGGTGACTAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      >  2:7127001/1‑150 (MQ=255)
                                                                                                                                 ttttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      >  2:2340510/1‑150 (MQ=255)
                                                                                                                                 ttttGGTTACTTTGGCCACGAAAAATTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAagcccaaaagcccaaaagcccaaaagcccaaaagcccaaa                      <  1:1608281/150‑1 (MQ=255)
                                                                                                                                    tGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCt                   <  2:8072038/150‑1 (MQ=255)
                                                                                                                                      gTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 >  2:3091362/1‑150 (MQ=255)
                                                                                                                                      gTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 >  1:7388361/1‑150 (MQ=255)
                                                                                                                                      gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 >  1:136752/1‑150 (MQ=255)
                                                                                                                                      gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 >  2:2597587/1‑150 (MQ=255)
                                                                                                                                      gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCttt                 >  1:7534073/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  2:3503569/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                >  1:64476/1‑150 (MQ=255)
                                                                                                                                       ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAACAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  2:123106/150‑1 (MQ=255)
                                                                                                                                       ttACTTTGGCCACCACAAAATTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  1:771869/150‑1 (MQ=255)
                                                                                                                                       ttACTTTGCCCACGACAAAATAACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTa                <  1:2171123/150‑1 (MQ=18)
                                                                                                                                          cTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAgcg             >  2:7324243/1‑150 (MQ=255)
                                                                                                                                             tGGCCACTACAAAATAACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCg          <  1:3721294/150‑1 (MQ=18)
                                                                                                                                             tGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCg          <  2:896703/150‑1 (MQ=255)
                                                                                                                                               gCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGt        >  2:6054662/1‑150 (MQ=255)
                                                                                                                                                 cACGACCGAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTc      >  2:6801990/1‑150 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTc      <  2:5842873/150‑1 (MQ=255)
                                                                                                                                                 cACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAAGGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTc      <  2:8716507/150‑1 (MQ=255)
                                                                                                                                                     aCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAt  <  2:2845417/150‑1 (MQ=255)
                                                                                                                                                     aCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAt  >  1:4185356/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTTGATGCGAAAGCGCACAACCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAAGCCCAAAGCTTTAGCGCCGGTTCTGAT  >  NC_002947/2092683‑2092981

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: