Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 9.8% intergenic (‑358/+137) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2150TG 9.8% 332.6 / 9.7 128intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (83/30);  new base G (2/13);  total (85/43)
Fisher's exact test for biased strand distribution p-value = 1.04e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCT  >  NC_000913/915066‑915362
                                                                                                                                                      |                                                                                                                                                   
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGAtt                                                                                                                                                     <  1:112172/150‑1 (MQ=32)
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGAtt                                                                                                                                                     <  2:3635829/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGCGGCGGCGAtt                                                                                                                                                     <  2:301295/150‑1 (MQ=17)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTg                                                                                                                                                    <  2:3400769/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCACCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTg                                                                                                                                                    <  2:618867/150‑1 (MQ=17)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTgg                                                                                                                                                   <  1:577135/150‑1 (MQ=32)
  ctcAGCCCACATCCATGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTgg                                                                                                                                                   >  1:2842329/3‑150 (MQ=12)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgc                                                                                                                                                  <  1:2880550/150‑1 (MQ=32)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgcc                                                                                                                                                 <  2:1420599/150‑1 (MQ=32)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgcc                                                                                                                                                 <  2:1421900/150‑1 (MQ=32)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgcc                                                                                                                                                 <  2:696805/150‑1 (MQ=32)
     aGCCCACATGGATGTGGGCTGAGGGCGCGTGTTACAGGGATGTTACCTCGCGCCCGACCCGGTCGCCGTACGGGATCAGTCGCGGGCACCGCGCAGCGGCGCTTTTGTTCGCCAGAGCCCGGGGGTGCT‑GGGGGCGGCGGCGAGGGgccg                                                                                                                                                >  2:3153444/1‑150 (MQ=11)
     aGCCCACATCCATGTGGGCTGAGGGCGCGTTTTAAAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgccg                                                                                                                                                >  1:3274363/1‑150 (MQ=12)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgccgc                                                                                                                                               >  2:220307/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgccgcc                                                                                                                                              >  1:2775330/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGgccgcc                                                                                                                                              >  2:2403476/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGccc                                                                                                                                             <  1:2357886/150‑1 (MQ=32)
         cactTGGTTGGGGGCTGTGGGCGCTTTTTACTGGGATGTTACCTCGCGCCCGACCCGGTAGCCTTATGGGATATGTGGTGGGCACCGCGCAGCGGCGTTTTTGTTCGCCAGAGCCGGGGGGTGCT‑GGGGGGGGCGGCGATTGGCCGcccc                                                                                                                                            <  1:1123521/150‑1 (MQ=11)
           caTGGTTGTGGGCTGTGGGCGCGTTTTACTGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCGCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCT‑GGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                          <  1:2641888/150‑1 (MQ=11)
           caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                          <  2:3558267/150‑1 (MQ=32)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                >  1:799754/1‑150 (MQ=32)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                >  1:799774/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                              >  1:641044/1‑150 (MQ=32)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                              >  2:1735491/1‑150 (MQ=32)
                        tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                             >  2:840096/1‑150 (MQ=32)
                          aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgcgc                                                                                                                           >  2:3837629/1‑150 (MQ=32)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                       >  2:2780043/1‑150 (MQ=32)
                               gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg                                                                                                                      >  1:3670895/1‑150 (MQ=32)
                                cgTTTTACGGGGTTGTCCCCGCGCGCCCGCCCGGGTCGCGGTCCGGGATACGCGGTGGGCCCCGCGCCGCGGCGTTTTGGTTGGCCTGCGCCGGGGGGGGCT‑GGGGGCGGCGGGGATGGGCCGCCCCGGGCGCGCTCCTTGCGCCAGTgg                                                                                                                     <  2:3138669/150‑1 (MQ=255)
                                     tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAata                                                                                                                >  2:2552816/1‑150 (MQ=37)
                                      aCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAatat                                                                                                               >  1:2343620/1‑150 (MQ=38)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                             >  1:2724119/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAatattt                                                                                                             >  1:1226325/1‑148 (MQ=255)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCACGGGCAATATGTTg                                                                                                           >  2:1452325/1‑150 (MQ=255)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCTTGCGCCAGTGGCAATATGTTg                                                                                                           >  1:2672443/1‑150 (MQ=255)
                                             tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTTGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTTCtt                                                                                                        >  1:3254096/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGTTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCCATATGTTGCTTAg                                                                                                      >  2:381038/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                      >  2:3068609/1‑150 (MQ=255)
                                               tccccTCGCGCCCGACGGGTTCGCCTAATGGCATCTGGCGGGGCCCCCGCGCCGGGGTGTTTTTGTCCGCCAGCGCGGGGGGGGGCG‑GGGGGCGGGGGTTTTGGGCCGCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                      <  2:2081205/147‑1 (MQ=9)
                                                 aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                    >  2:1108093/1‑150 (MQ=255)
                                                   cTCGCGCCCGACCCGGTAGCCGTAAGGGAAAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCACGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTAGCCCCAGTGGCAATATGTTGAGTAGCTCa                                                                                                  >  1:2391213/1‑150 (MQ=14)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                 >  2:150717/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                 >  1:1078286/1‑150 (MQ=255)
                                                     cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                >  1:2416206/1‑150 (MQ=255)
                                                         cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAg                                                                                            >  2:1611536/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCTCGCAGCGGCGATTATGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCACCGAGCGCCAGTGGCAATATTTTGCTTAGCTCAAGAAAgg                                                                                           >  1:2424617/1‑150 (MQ=17)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATTAAAGGAg                                                                                         >  1:316589/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                         >  2:3379458/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                         >  2:2169825/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCTCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                        <  2:4211700/150‑1 (MQ=255)
                                                               ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCATGCGCCCCTGGCAAAATGTTTCTTAGCTCAAGAAAGGAgcgc                                                                                      >  1:1309282/1‑150 (MQ=17)
                                                                cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑TGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCAATCGCCAGGGGCAAAATGTTGCTTAGCTCAAGAAAGGAGCGca                                                                                     >  1:343176/1‑150 (MQ=17)
                                                                   tAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCACTGGCAATATGTTGCTTAGCTCATGAAAGGTGCGcaaca                                                                                  >  1:1781318/1‑150 (MQ=255)
                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCCGCGCGCTCCGTTCGCCAGTGGCAAAATGTTGCTTAGCGCAAGAAACGAGCGCCACAAga                                                                               >  1:1559591/1‑150 (MQ=17)
                                                                          atgggATACGTCGTGGGCACGGGGCTGGGGGGTTTTTGTTGGCCTGAGCCGGGGGGTGCT‑GGGGGGGGGGGCGATTGGCGGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                           <  2:276559/148‑1 (MQ=18)
                                                                          acgggATACGGCGGGGGCCCCGCGCCGCGGCGATTTTTTTCGCCTGTGCCGGGGGGGGCC‑GGGGGGGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                           <  2:276543/148‑1 (MQ=17)
                                                                           aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCCCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATTTTTCTTATCTCATTAAAGGAGCGCAAAAAGattatg                                                                          >  1:4254295/1‑150 (MQ=17)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCCGCGCGCGCCGTGCGCCCGGGGCCAAAAGTTGCTTATCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                       >  1:845076/1‑150 (MQ=18)
                                                                                 aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                    >  2:1789616/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCC‑GGGGGCGGCGGCGCTTGGCCCCCCCCAGCGCGCGCCGTGCGCCAGGGGCCAAAATTTTCTTATCACATGAAAGGAGCGCAACAACATGATGAATCaaaaa                                                                 >  1:907222/1‑147 (MQ=17)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                 >  2:598729/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                 >  2:4130137/1‑150 (MQ=255)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                               >  2:2402825/1‑150 (MQ=255)
                                                                                        gggCACCGCGCTGGGGGGTTTTTGTTCGCCAGAGCCGGGGGGTGCT‑GGGGGCGGGGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                             <  2:2775321/150‑1 (MQ=21)
                                                                                        gggCACCGCGCCGCGGCGTTTTTGTTCGCCAGAGCCGGGGGGTGCT‑GGGGGGGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                             <  2:2775330/150‑1 (MQ=21)
                                                                                            aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                         >  2:876800/1‑150 (MQ=255)
                                                                                            aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                         >  2:435442/1‑150 (MQ=255)
                                                                                                      ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                               >  2:711641/1‑150 (MQ=255)
                                                                                                        cttttttGTTCGCCTGGGCCGGGGGGGGCT‑GGGGGGGGGGGGTTTGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                             <  2:497395/147‑1 (MQ=255)
                                                                                                           ttttGTTCGCCTGAGCCCGGGGGTGCA‑GGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                          <  2:359420/150‑1 (MQ=255)
                                                                                                             ttGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                                        >  2:3218566/1‑150 (MQ=255)
                                                                                                             ttGTCCGCCAGAGCCGGGGGGTGCT‑GGGGGCGGCGGCGATGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                                        <  2:644936/150‑1 (MQ=37)
                                                                                                                   gCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAAAGTTGCTTAGCTCATGAAAGGAGAGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCAGAAAtt                                  >  2:2528503/1‑150 (MQ=25)
                                                                                                                   gCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCACCTTGCGCCCGTGGCAATATGTTGCTTAGCGCATGAAAGGACCGCAACAAGATGAAGAAGAAACAAATAACAACATCATAAAAAAAGGCCTGAccat                                  >  2:1488796/1‑147 (MQ=21)
                                                                                                                   ccctgcgCCGGGGGGGGGT‑GGGGGGGGGGGCGCGGGCCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                  <  1:612142/144‑1 (MQ=255)
                                                                                                                      tgcgcgggggggTGTT‑GGGGGGGGGGGCTTTGGCCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                               <  2:158078/143‑1 (MQ=37)
                                                                                                                      agagCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATACCAACATCTTAAAAAAAGGCCTGACATTACg                               >  2:3851673/1‑150 (MQ=255)
                                                                                                                         gcgggggggTGCG‑GGGGGGGGGGGAGTTGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                            <  2:2568693/146‑1 (MQ=38)
                                                                                                                         gCCCGGGGGTGCT‑GGGGGGGGCGGCTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                            <  2:3176302/150‑1 (MQ=255)
                                                                                                                          cgggggggTGGG‑GGGGGGGGGGGCGTGGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                           <  1:3145072/147‑1 (MQ=38)
                                                                                                                            gggggggGCT‑GGGGGCGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                         <  2:1073881/149‑1 (MQ=255)
                                                                                                                             gggggTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                        >  2:1100130/1‑150 (MQ=255)
                                                                                                                             gggggTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                        >  2:2987885/1‑150 (MQ=255)
                                                                                                                                  tGCT‑GGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                   <  2:1078286/150‑1 (MQ=255)
                                                                                                                                  tGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                   >  2:785188/1‑150 (MQ=255)
                                                                                                                                    cA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                 >  2:1797933/1‑150 (MQ=255)
                                                                                                                                    cA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                 >  1:3732115/1‑150 (MQ=255)
                                                                                                                                    cA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                 >  2:1797949/1‑150 (MQ=255)
                                                                                                                                    cA‑GGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                 >  1:3732080/1‑150 (MQ=255)
                                                                                                                                      gggggggggggGCGATTGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                <  1:4191576/149‑1 (MQ=255)
                                                                                                                                        ggggCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta              <  2:2418812/150‑1 (MQ=255)
                                                                                                                                         gggCGGCGGCGTTGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa             <  1:1751599/150‑1 (MQ=255)
                                                                                                                                          ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat            >  1:128126/1‑150 (MQ=255)
                                                                                                                                           ggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt           >  2:2663431/3‑150 (MQ=255)
                                                                                                                                           gcggcggcgTTGGCCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt           <  1:3161239/150‑1 (MQ=255)
                                                                                                                                                ggtgtttGCCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      <  2:3978588/145‑1 (MQ=255)
                                                                                                                                                ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      >  1:631074/1‑150 (MQ=255)
                                                                                                                                                ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa      >  1:2350586/1‑150 (MQ=255)
                                                                                                                                                   gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc   >  2:1211344/1‑150 (MQ=255)
                                                                                                                                                    aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt  >  2:732406/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCA‑GGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCT  >  NC_000913/915066‑915362

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: ‑