Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 17.5% intergenic (‑369/+126) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2260TG17.5% 160.8 / 3.5 90intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (49/10);  new base G (8/8);  total (70/20)
Fisher's exact test for biased strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGG  >  NC_000913/915078‑915373
                                                                                                                                                     |                                                                                                                                                   
aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTg                                                                                                                                                    >  2:1805422/1‑150 (MQ=32)
      gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTgcgcgcg                                                                                                                                              >  1:558610/1‑149 (MQ=17)
       tgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTc                                                                                                                                             >  2:1530347/1‑150 (MQ=32)
          gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                          >  2:1844941/1‑150 (MQ=32)
          gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                          >  2:2267824/1‑150 (MQ=32)
          gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                          >  2:2267846/1‑150 (MQ=32)
          gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                          >  2:2267983/1‑150 (MQ=32)
           cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                                         >  2:202199/1‑150 (MQ=32)
           cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCGGCGCGCTCCTTg                                                                                                                                         >  2:70562/1‑150 (MQ=17)
            tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                                        >  2:1201956/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGcccg                                                                                                                                   >  1:1706035/1‑148 (MQ=17)
                   gctttttACTGGTTTGATACCACGGGCCCCCCCCGTTTTCCGTACGGGACCGGTCGTGGGCTCGGCCCACCGGCGTTTTTGTGCCCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCCGCGCGCTCCTTGCGCCAGTg                                                                                                                                 <  1:1385946/147‑1 (MQ=2)
                   gctttttACTGGGTTGTTACCTCGCGCCCGACCGGGCAGCCGTATGGGATACGTCGCGGGCACCGCGCCGCGGCTATTTTGTCCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCGAT‑TGGCCGCCCCCCGCGCGCTCCTTGCGCCAGTg                                                                                                                                 <  1:1385870/147‑1 (MQ=11)
                   gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg                                                                                                                                 >  2:1536108/1‑150 (MQ=32)
                    cgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTgg                                                                                                                                >  1:2021884/1‑150 (MQ=32)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGAT‑TGGCCGCCCCCGGCGCGCTCCTTGCGCCCGGGGCAATATGt                                                                                                                        >  1:772041/1‑150 (MQ=11)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                        >  1:1784555/1‑150 (MQ=39)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                        >  2:253299/1‑150 (MQ=255)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGCAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                        >  2:1792267/1‑150 (MQ=255)
                              ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                      >  2:2843435/1‑150 (MQ=40)
                              ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                      >  2:2259773/1‑150 (MQ=40)
                                 tGTTACCTCGCGCCGGACCGGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGGGAT‑TGGCCACCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                                   <  2:1570403/150‑1 (MQ=14)
                                 tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                                   >  2:1761779/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  1:1605929/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  1:2652903/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  2:489158/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  2:2261035/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  2:2866870/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  1:225348/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                 >  2:37097/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCCGGCGCGCGCCGTGCGCCAGGGGCAATATGTTGCGTAg                                                                                                                 >  1:1971771/1‑150 (MQ=14)
                                    tACCTCGCGCCCGACCGGTTAGCCTTATGGTATACGTCGTGGGCACCGCGCAGCGGGTTTTTTTTTCGCCAGAGCCGGGGGGGGCCGGGGGGGGGGGCGAT‑TGGCCGCCCCATGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                                                <  2:240983/150‑1 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTATCTCa                                                                                                             >  2:2423949/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                             >  2:2423971/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCGGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                             >  1:1121667/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAAAATGTTGCTTAGCTCa                                                                                                             >  1:640254/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGACAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGGTCa                                                                                                             >  2:783932/1‑150 (MQ=255)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                            >  2:266847/1‑150 (MQ=255)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                            >  2:266936/1‑150 (MQ=255)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCCATATGTTGCTTAGCTCAt                                                                                                            >  2:266977/1‑150 (MQ=255)
                                              ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                                      >  2:2190726/1‑150 (MQ=255)
                                               cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                                     >  2:1554830/1‑150 (MQ=255)
                                                 aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                                   >  2:62032/1‑150 (MQ=255)
                                                  cccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                                  >  2:1863563/1‑150 (MQ=255)
                                                          ccGTAAGGGTTACGGCGAGGGCCCCGCGCAGCGGCGTTTTTGTTCGCCTGAGCCGGGGGGTGCCGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                          <  2:1871844/150‑1 (MQ=17)
                                                            gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCCCCCCCAGCGCCCTCCATGCGCCAGTGGCAAAAATTTGGTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                                        >  1:957524/1‑150 (MQ=17)
                                                             tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                                       >  2:2065760/1‑150 (MQ=255)
                                                             tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGCT‑TGGCCGCCCCCCGCGCCCCCCCTGCCCCCCAGGCACAATTTTTCTTATCTAAAGAAAGGAGCGcaacaacaata                                                                                       >  1:1229846/1‑145 (MQ=14)
                                                             tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGCG‑TGGCCCCCCCCCGCGCGCGCCGTGCGCCCCGGGCAAAATTTTTATTATGTCAAGAAATGAGCGcaacaaaaata                                                                                       >  1:1229825/1‑145 (MQ=2)
                                                               aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TTGCCGCCCCCAGCGCGCTCCGTGCGCCAGGGGCAAAATGTTGCTTAGCTCAGGAAAGGAGCGCAACAAAatgatg                                                                                     >  1:2839619/1‑150 (MQ=18)
                                                               aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCGCATGAAAGGAGCGCAACAAGatgatg                                                                                     >  1:2839612/1‑150 (MQ=255)
                                                                 ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGCT‑TTGGCCCCCCCAGCGCGCGCCTTGCGCCCCTGGCCAAAATTTGCTTAGGTCCTGAAAGGAGCGCAACAAGatgatgta                                                                                   >  1:779746/1‑148 (MQ=17)
                                                                 ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TTGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTATCTCATGAAAGGAGCGCAACAAGATGATGaa                                                                                   >  1:2383411/1‑150 (MQ=255)
                                                                  gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGA‑TGGCCGCCCCCAGCGCGCTCCTTGCGCCCCGGGCCAAATGTTGCGTAGCTCATGAAAGGAGCGCAACAAGAAGATGaaa                                                                                  >  1:2550218/1‑149 (MQ=17)
                                                                     aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCCCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGTCAGGAGCGCAACAAGATGATGAATCaa                                                                               >  2:1271883/1‑150 (MQ=255)
                                                                       gTCGCGGGCACCGGGCTGGGGGTTTTTTGGTCGCCAGAGCGGGGGGGGGCGGGGGGGGGGGGCGTT‑GGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                             <  2:183266/150‑1 (MQ=18)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAATATGTTTGTTATGTAAAGAAAGGAGCGCAACAAGATGATTAATAAACaa                                                                            >  1:2880492/1‑149 (MQ=17)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGAG‑TGGCCGCCCCCAGCGCCCCCCGTGCGCCCCTGGCAAAATTTTGCTGTTCTCAGGTAAGGAGCGCAACAACATGGTTACTCAACaa                                                                            >  1:634730/1‑149 (MQ=9)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCAGCGCGCTCCATGCGCCAGGGGCAAAATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGTATCAACat                                                                            >  1:2840843/1‑150 (MQ=18)
                                                                          gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACCAGATGATGAATCAACatat                                                                          >  2:1875032/1‑150 (MQ=255)
                                                                            gggCACCGCGCAGCGGCGTTTTTGTTCGCCTGTGCCGGGGGGTGCTGGGGGCGGGGGGTTT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                        <  2:2404754/150‑1 (MQ=21)
                                                                                    cgcAGGGGCGTTTTTGTTCGCCTGTGCCGGGGGGTGCTGGGGGGGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCt                                                                <  2:2120633/150‑1 (MQ=21)
                                                                                    cgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGAATAAACATATAACAACATCt                                                                >  2:1760988/1‑150 (MQ=255)
                                                                                       aGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TTGCCGCCCCCTTGGCGCTCCTTTCGCCAGTGGCAATATGTTTCTTATCTCATGAAAGGAGCGCAAAAAGATTATTAATCAAAATATAAAAACATATTaa                                                             >  2:1252798/1‑150 (MQ=18)
                                                                                           gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCAGCGCGCTCCTTGCGCCAGAGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACATAAAACAACATCTTaaaaaa                                                         >  2:2836330/1‑150 (MQ=255)
                                                                                            cgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGATGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAACTTaaaaaaa                                                        >  2:511674/1‑150 (MQ=255)
                                                                                               ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                                     >  2:2441410/1‑150 (MQ=255)
                                                                                                    ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCGGCGCCCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACCTATAACAACATCTTAAAAAAAGGccgggc                                                >  2:1015938/1‑146 (MQ=255)
                                                                                                        ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACATCATCTTAAAAAAAGGCCTGACATTa                                            >  2:156830/1‑150 (MQ=255)
                                                                                                        ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                            >  2:436152/1‑150 (MQ=255)
                                                                                                         ctgcgCCCGGGGGTGCTGGGGGGGGC‑GCGTTGTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                           <  2:1453326/146‑1 (MQ=255)
                                                                                                                       gctgggggcgggggcTTT‑TGGCCGCCCCCTGCGCGCTCCTTGAGCCAGTGGCACTATGTTGCTTATCTCATGAAATGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTATAAAAAGGCCTGACATTACGCCTTGCCTTCTGc                             <  2:662785/147‑1 (MQ=255)
                                                                                                                       gctgggggCGGGGGGGTG‑TGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                             <  1:2601809/147‑1 (MQ=38)
                                                                                                                            gggCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                        >  2:1814827/1‑150 (MQ=255)
                                                                                                                                 ggggcgCT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                   <  1:422574/148‑1 (MQ=255)
                                                                                                                                    ggggtggggCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                 <  1:1547195/143‑1 (MQ=255)
                                                                                                                                      gAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc              >  1:430608/1‑150 (MQ=255)
                                                                                                                                      gAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc              >  1:517096/1‑150 (MQ=255)
                                                                                                                                       gtgggcccgccCCGTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc              <  2:2446225/144‑1 (MQ=255)
                                                                                                                                            cccccccccGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc         <  1:1015938/146‑1 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCATCATATAACAACATCTTAAAAAAAGGCCTGACCTTACGCCAGGCCTTCTGCGTGCAGTAATCACGCTTTTcc        >  1:1601091/1‑150 (MQ=39)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTGATCACGCTTTTcc        >  1:2366897/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  1:587800/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  2:778615/1‑150 (MQ=255)
                                                                                                                                              cgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca       >  2:1141350/1‑150 (MQ=255)
                                                                                                                                              cccccccGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca       <  1:2259773/148‑1 (MQ=255)
                                                                                                                                              cccccccGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca       <  1:55275/148‑1 (MQ=255)
                                                                                                                                                  ccccGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag   <  1:2467014/150‑1 (MQ=255)
                                                                                                                                                  cccGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag   <  1:2261035/150‑1 (MQ=255)
                                                                                                                                                   ccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg  <  2:453681/150‑1 (MQ=255)
                                                                                                                                                   ccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg  <  2:453682/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGG  >  NC_000913/915078‑915373

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑