Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,090,581 C→T 5.8% intergenic (+298/‑297) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,090,5810CT 5.8% 585.9 / 18.6 208intergenic (+298/‑297)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (141/55);  new base T (7/5);  total (148/60)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 215 total aligned reads displayed.

GGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2090435‑2090730
                                                                                                                                                  |                                                                                                                                                     
ggTGACGGGCTTGGAGGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCAc                                                                                                                                                    <  2:2851353/150‑1 (MQ=21)
ggTGACGGGCATGGAGGCCCGTCAAGCGCTGGGCCCCCGGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCAc                                                                                                                                                    <  1:2716167/150‑1 (MQ=255)
 ttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                   >  2:2085266/2‑150 (MQ=21)
 ttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                   >  1:4210565/2‑150 (MQ=21)
 ttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                   >  2:1266135/2‑150 (MQ=21)
 gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACg                                                                                                                                                   >  2:3015721/1‑150 (MQ=34)
   gACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAc                                                                                                                                                 >  1:5117807/1‑150 (MQ=34)
   gACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAc                                                                                                                                                 >  2:6097126/1‑150 (MQ=21)
    aCGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGAGGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                >  2:5642457/1‑150 (MQ=21)
    aCGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAcc                                                                                                                                                >  2:4301772/1‑150 (MQ=21)
      gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGAccca                                                                                                                                              >  1:6757447/1‑148 (MQ=21)
        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTCGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                            >  1:1921885/1‑150 (MQ=21)
        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAg                                                                                                                                            >  2:6382325/1‑150 (MQ=34)
          ctggtggCCCGTCACGCGCTGGGCCCCCGGCTGGGGCCTGCCGCGCGGTCCGCCGGGCAGCACGCCGGGAGTGTGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTa                                                                                                                                          <  1:3718154/144‑1 (MQ=11)
          aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTa                                                                                                                                          >  1:690099/1‑150 (MQ=34)
           tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTaa                                                                                                                                         <  2:3754784/150‑1 (MQ=255)
           tGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTaa                                                                                                                                         <  2:5033572/150‑1 (MQ=34)
                gCCCGTCACGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCGCCCGGGAGAAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTTGTTACTTTGGTCGTGGCCAAAGTAAccaaa                                                                                                                                    <  1:4944516/150‑4 (MQ=11)
                 cccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACccgcc                                                                                                                                   >  1:3343963/1‑150 (MQ=34)
                 cccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACccgcc                                                                                                                                   >  1:3574590/1‑150 (MQ=34)
                  ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACccgccg                                                                                                                                  >  1:595196/1‑150 (MQ=34)
                   cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCgg                                                                                                                                 >  1:5413774/1‑150 (MQ=34)
                     tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGaa                                                                                                                               >  2:1425889/1‑150 (MQ=34)
                       aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAgg                                                                                                                             >  1:3059625/1‑150 (MQ=34)
                            gggggCCCCAGGATGGGGCCTGCCGCGGGGTCCGCCGGGGTGCACGCCCGGTTTGTGGGGACCCCGGAGCGTAGCGTTGGGGCGGGATGAGGGGCGCCGACGGCTTTGGTTACTTTGGCCACGACCAAATTAACCCGCCGGAAGGGCGGa                                                                                                                        <  2:1091855/148‑1 (MQ=11)
                                gCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAgg                                                                                                                    >  1:5678919/1‑150 (MQ=34)
                                  cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAACGGTg                                                                                                                  >  1:6100375/1‑150 (MQ=21)
                                  cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCAGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGACAggag                                                                                                                  >  1:3816227/1‑148 (MQ=19)
                                   ccatcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  1:456467/7‑150 (MQ=21)
                                   ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  1:2389940/1‑150 (MQ=34)
                                   ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  2:4811786/1‑150 (MQ=34)
                                   ccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                                                                                                 >  2:6482941/1‑150 (MQ=34)
                                    cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTCGGGGCCGTATGTGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGaa                                                                                                                <  2:5127918/150‑1 (MQ=255)
                                    cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGACAGGTGaa                                                                                                                >  2:3277662/1‑150 (MQ=21)
                                     atcctGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                               >  2:332993/5‑150 (MQ=21)
                                     aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGTCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAAt                                                                                                               >  2:774445/1‑150 (MQ=21)
                                           gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGgcgc                                                                                                         >  1:1125451/1‑150 (MQ=34)
                                            ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGcc                                                                                                        >  1:3012225/1‑150 (MQ=34)
                                            ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGcc                                                                                                        >  1:1310380/1‑150 (MQ=34)
                                                tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGGAAGGGGAATAGGCGCCTGCa                                                                                                    >  1:2899696/1‑150 (MQ=255)
                                                 gCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCAt                                                                                                   >  1:1611717/1‑150 (MQ=37)
                                                    gcgcGGTCCTCCCGGGAGCACGCCGGGCGTGCGGGGACCCCGGAGCGTCGCGTAGGGGCCGGATGCGGGGAGCCGCCGGCTTTGTTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt                                                                                                <  2:489770/150‑1 (MQ=11)
                                                       cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                             >  2:3085691/1‑150 (MQ=37)
                                                       cgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                             >  1:1630092/1‑150 (MQ=37)
                                                       cgGTCCGCCCGGGAGCATGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGCCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGtgaa                                                                                             <  2:6265625/150‑1 (MQ=255)
                                                          tcctccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGACAGGTGACTAGGCGCCTGCATCGtgaatga                                                                                          >  2:3774138/1‑150 (MQ=255)
                                                            ctccCGGGAGCAAGCCGGGTGTGAGGGGACCCCGGAGCGTAGGGTAGGGGCGGGATGTGGGGCGCCGACGGCTTTGTTTACTTTGGCCACGACCAATTTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        <  2:6698285/150‑1 (MQ=11)
                                                            cccccGGGAGCACGCCGGGCGTGGGGGGACCCGGGCGCGTCGCGTCGGGGCGGGTTGTGGGGAGCCGAGGGCTTTGGTTACTTGGGCCACGACCAATGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat                                                                                        <  2:2499718/148‑1 (MQ=11)
                                                             tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGGAAGGTGAATAGGCGCCGGCATCGTGaatgaatg                                                                                       >  1:798638/1‑150 (MQ=255)
                                                             tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                       >  2:2969263/1‑150 (MQ=37)
                                                             tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg                                                                                       >  2:3611448/1‑150 (MQ=37)
                                                               ccGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTc                                                                                     >  2:442209/1‑150 (MQ=37)
                                                                       aaGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTATACTTGAt                                                                             >  1:5192558/1‑150 (MQ=37)
                                                                        aGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtg                                                                            >  2:3462299/1‑150 (MQ=37)
                                                                         gCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtat                                                                           >  1:6696226/1‑148 (MQ=37)
                                                                             gttgtgtGGGGACCCGGGTGGGTTGGGTTGGGGGCGGTTGTGGGGCGCCGCGGGTTTTGGTTATTTGGGCCACCAAAAAACAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGaaa                                                                       <  1:3870443/147‑1 (MQ=2)
                                                                             gtcgtgtGGGGACCCGGGAGCGTTGGGTTGGGGCCGTATGTGGGGAGCGGACGGCTTTTGTTACTTGGGCCACCACCAATGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGaaa                                                                       <  2:2033125/147‑1 (MQ=14)
                                                                               cggggggggCCCCGGGCGCGTTGCGGGGGGGCCGTTTGAGGGGCGCCGCCGGCTTTGTTTACTTTGGCCACGCCCAACGAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgc                                                                     <  2:4914306/145‑1 (MQ=14)
                                                                                    ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                                >  1:4730507/1‑150 (MQ=39)
                                                                                    ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGcaca                                                                >  1:1017023/1‑150 (MQ=39)
                                                                                    ggggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCAACGTGAATGAAAGTCTAATTTATGTGAAAGCGcaca                                                                >  1:1775778/1‑150 (MQ=255)
                                                                                      ggACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGGAAGGGGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTTAAAGCGCACAAc                                                              >  1:1367829/1‑150 (MQ=255)
                                                                                          cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCaa                                                          >  2:1594924/1‑150 (MQ=39)
                                                                                          cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTATACTTGATTTGAAAGCGCACAACCCaa                                                          >  1:576438/1‑150 (MQ=255)
                                                                                           ccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAAAGAATGTCTACTTTATGTTAAAGCGCACAACCCaaa                                                         >  1:5318022/1‑150 (MQ=255)
                                                                                                   tagcgtagGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaa                                                 >  2:3470322/1‑150 (MQ=255)
                                                                                                      gtttggggCCGTATGTGGGGAGCCGACGTCTTTGTTTATTTTGGCCACGACAAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagc                                              <  2:3518096/146‑1 (MQ=14)
                                                                                                      ggtgggggCCGTTTGGGGGGAGCCGTGGGTTTTGTTTTTTTGGCCCACCACAAACGAACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagc                                              <  2:2260109/146‑1 (MQ=11)
                                                                                                            ggCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaagctcaaaa                                        >  1:6190318/1‑150 (MQ=255)
                                                                                                                 tttttggggCGCCGTCGGGTTTGGTTAGTTTGGCCTCGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGATAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCttt                                   <  2:4680988/145‑1 (MQ=14)
                                                                                                                      ggggAGCCGCGGGTTTTTTTTACTTGGGCCACGACAAATTAAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              <  2:1200546/150‑1 (MQ=14)
                                                                                                                      ggggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              <  2:4168119/150‑1 (MQ=255)
                                                                                                                      ggggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGGGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              >  1:40426/1‑150 (MQ=255)
                                                                                                                      ggggAGCCGACGGCTTTGGTTACTTGGGCCACGACCAACGTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc                              <  2:2624099/150‑1 (MQ=255)
                                                                                                                         gAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCGAAAGCGTTAGCGCCgg                           >  2:3814932/1‑150 (MQ=255)
                                                                                                                         gAGCCGACGACTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCgg                           <  2:4205998/150‑1 (MQ=255)
                                                                                                                              gACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTg                      >  2:6591775/1‑150 (MQ=255)
                                                                                                                                   cttggtttACTTTGCCCACGACCAAATTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGCCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg                 <  2:2431000/144‑1 (MQ=17)
                                                                                                                                    tttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCggg                >  1:5053964/1‑150 (MQ=255)
                                                                                                                                     ttGTTTACTTGGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc               <  1:6167088/150‑1 (MQ=255)
                                                                                                                                     ttGGTTAGTTTGCCCCCCACAAACTTACCCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc               <  1:185419/150‑1 (MQ=18)
                                                                                                                                     ttGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc               <  2:4210565/150‑1 (MQ=255)
                                                                                                                                     tgggTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc               <  2:690099/148‑1 (MQ=255)
                                                                                                                                       ggTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt             >  2:1780121/1‑150 (MQ=255)
                                                                                                                                       ggTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGt             >  2:4563686/1‑150 (MQ=255)
                                                                                                                                        gTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTa            >  1:1155530/1‑150 (MQ=255)
                                                                                                                                        gTTACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTa            >  1:3963727/1‑150 (MQ=255)
                                                                                                                                        gTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAGGTGAAAGCGCACAACCCAAAGGCCCAACAGCTCAAAAGCGTTAGCGCCCGTGCTGATCGGGCgga            >  2:6857360/1‑148 (MQ=17)
                                                                                                                                         ttACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTAc           >  2:5730510/1‑150 (MQ=255)
                                                                                                                                          tACTTTGGTCGTGGCCAAAGTAACCCGCCGGAAGGGCGGTAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTTATGTTAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACa          >  1:904969/1‑150 (MQ=18)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGACAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTTATGTGAAAGCGCACAACACAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTc       >  1:484177/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTATCGCCGGTTCTGATCGGGCGTACAGTc       >  2:2935130/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTc       >  2:6749811/1‑150 (MQ=255)
                                                                                                                                             tttGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTc       >  1:2427474/1‑150 (MQ=255)
                                                                                                                                                ggCCACGACCAACGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGt    <  2:1247547/150‑1 (MQ=255)
                                                                                                                                                ggCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGt    <  2:510325/150‑1 (MQ=255)
                                                                                                                                                  ccACGCCCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGttt  <  2:5497306/150‑1 (MQ=255)
                                                                                                                                                  ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGCATAGGCGCCGGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACCGTCCGttt  >  1:6615653/1‑150 (MQ=255)
                                                                                                                                                  ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGttt  >  1:424498/1‑150 (MQ=255)
                                                                                                                                                  ccACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGCCCGttt  >  1:6606953/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTT  >  NC_002947/2090435‑2090730

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: