Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 8.3% 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,2260TG 8.3% 771.5 / 4.1 359intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (227/52);  new base G (13/19);  total (269/90)
Fisher's exact test for biased strand distribution p-value = 2.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 374 total aligned reads displayed.

CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGT  >  NC_000913/915077‑915375
                                                                                                                                                       |                                                                                                                                                     
caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                                       <  2:5085919/150‑1 (MQ=32)
   ggctgtgGGCTGAGGGCGCGTTTTACAGGGATGGTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                                    <  2:9524376/147‑1 (MQ=17)
   ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                                    >  2:7584007/1‑150 (MQ=32)
    gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc                                                                                                                                                   <  2:4483365/150‑1 (MQ=32)
      tgtgGGGTGTGGGCGGGTTTTACCGGGAGGTTACCGCGCGCCGGACCGGGTAGCCGTATGGGATACGTCGCGGGCACCGCGCCGGGGCGATTTTGTTCGCCATAGCCGGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcgc                                                                                                                                                 <  2:9078002/150‑1 (MQ=11)
        tgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTc                                                                                                                                               >  1:3592553/1‑150 (MQ=32)
                 ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCCTGCGCCa                                                                                                                                      >  1:9327379/1‑150 (MQ=12)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  1:1142079/1‑150 (MQ=32)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  2:7052999/1‑150 (MQ=32)
                   cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                                    >  2:9206314/1‑150 (MQ=32)
                   cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                                    >  2:10709697/1‑150 (MQ=17)
                      gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGc                                                                                                                                 >  1:1363266/1‑150 (MQ=18)
                       ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                >  1:8214018/1‑150 (MQ=34)
                       ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                >  2:5744666/1‑150 (MQ=34)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                          >  1:8806477/1‑150 (MQ=39)
                              gggTTGTTACCGCGCGCCCGCCCCGGGACCCGTTAGGGTTTAGTCGTGGGCCCCGCGGCGCGGCGTTTTTTGTCGGCATCGGCCGGGGGGGGCGGGGGGGGGCGGCGCTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                          <  2:4006508/150‑1 (MQ=255)
                              gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGtt                                                                                                                         >  2:2468698/1‑150 (MQ=38)
                               ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                        >  2:1496242/1‑150 (MQ=40)
                               ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                        >  1:2586115/1‑150 (MQ=40)
                                 aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                                      >  1:7755182/1‑150 (MQ=255)
                                    ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                   >  1:9815116/1‑150 (MQ=255)
                                    ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                   >  2:3152941/1‑150 (MQ=255)
                                      aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                                 >  2:3363225/1‑150 (MQ=255)
                                       ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGCTTGGGCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATAAGTTGCATAtctc                                                                                                                >  1:2108681/1‑150 (MQ=14)
                                        cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCCATATGTTGCGTAGCTCa                                                                                                               >  1:2059492/1‑150 (MQ=255)
                                        cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCCATATGTTGCTTAGCTCa                                                                                                               >  1:7390444/1‑150 (MQ=255)
                                        cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCCGGGGCAATATGTTGGTTAGCTCa                                                                                                               >  1:4014975/1‑150 (MQ=14)
                                        cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCAAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGGGGCAATATGTTGCTTAGCTCa                                                                                                               >  1:6026429/1‑150 (MQ=255)
                                         tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                              >  1:10614042/1‑150 (MQ=255)
                                         gcgcgCCCGCCGGGGTAGCGTTATGGTATATGTCGCGGCCACCGCGCCGGGGGTTTTTTGTTCGCCAGCGCCGGGGGGTGCGGGGGGGGGCGGGTTTGGGCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                              <  2:4513366/149‑1 (MQ=255)
                                             gcccgcccGGGTCGCCGTACGGTATAAGTGGCGGGCACCGCGCCGCGGCTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaaa                                                                                                          <  2:787357/150‑1 (MQ=17)
                                             gcccgcccCGGTCGCCGTAAGGGATACGTCGTGGCCACCGCGCCGCGGCGATTTTTTGCGCCAGCGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaca                                                                                                          <  2:7356569/150‑3 (MQ=17)
                                                        tAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGAGCGccacc                                                                                               >  2:978664/1‑149 (MQ=17)
                                                        tAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCCCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaaca                                                                                               >  2:2575626/1‑150 (MQ=255)
                                                           ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATTGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCAATATGTTTCATAGCTCATGAAAGGAGCGCAACAAga                                                                                            >  2:9333641/1‑150 (MQ=17)
                                                            cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                           >  2:10090485/1‑150 (MQ=255)
                                                            cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                           >  2:2959124/1‑150 (MQ=255)
                                                              tAAGGGATAAGTCGAGGGCACCTCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTTGCTCCGGACAGGAGCGCCACAAgatga                                                                                         >  2:7169002/1‑150 (MQ=17)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                                         >  2:489536/1‑150 (MQ=255)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                                         >  2:6423332/1‑150 (MQ=255)
                                                                  ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCCGCGCCCACCATGCGCCCCTGGCAAAAATTTTCGTATCTCAAGAAAGGGGCGCAACAAGAAGATGaa                                                                                     >  1:9529974/1‑150 (MQ=14)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCCCGCGCCATGCGCCAGGGGCAATATTTTGCATAGCTCATGAAAGGAGCCCAACAACAAGATGAAt                                                                                    >  1:4385333/1‑150 (MQ=17)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                                    >  2:477107/1‑150 (MQ=255)
                                                                    aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATc                                                                                   >  2:11053953/1‑150 (MQ=255)
                                                                      atgtggtgggCACCGCGCTGGGGCTTTTTTGTTCGCCTGTGCGCGGGGGTTCTGGGGGCGGGGGGGATTGGCGGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                                 <  2:6281317/143‑1 (MQ=18)
                                                                       aGTCGTGGGCACCGCGCCGCGGCGATTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                                                                                <  2:7435192/150‑1 (MQ=18)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCCCGCCATGCGCCCCTGGCAAAATTTTTCGTATCTCATGAAAGGAGCGCAACAACATTATGAATCAACaa                                                                              >  1:9104201/1‑149 (MQ=17)
                                                                           gtgggCCCGGCGCGGGGGCGTTTTGGTGCGCCTGCGCCGGGGGGTGGGGGGGGGGGGGGGGATGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                            <  2:2626193/148‑1 (MQ=21)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGTTTTGCCGCCCCCAGCGCGCGCCGTGCGCCCGTGGCAAAATGTTGGTTATCTCATGGAAGGAGCGCACCAAGATGATGAATCAACatat                                                                            >  1:10924672/1‑150 (MQ=18)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAGCGCGCTCCGTTCGCCACTGGCAAAATGTTTGTTAGCTACTTTAATGAGCGCAACAAAATGATGGATaaaaaaaa                                                                            >  1:324669/1‑145 (MQ=14)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                            >  2:1644297/1‑150 (MQ=255)
                                                                            gggccACCGCGCCGGGGGTTTTTTGTGGGCCCGTGCGGGGGGTTGGTGGGGGGGGGGGGGTTTGGCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAtata                                                                           <  2:6925900/146‑1 (MQ=19)
                                                                               gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCCATTGGCCCCCCCCAGCGCGCACCATGCGCCAGTGGCAATAATTTTCTTAGCTCCTGAAAGGAGCGCAACAAGATGATTAATAAAAAAATaaaa                                                                        >  1:7142695/1‑148 (MQ=17)
                                                                                   cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                                    >  2:10977701/1‑150 (MQ=255)
                                                                                          cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                              <  2:10856547/150‑1 (MQ=255)
                                                                                             gttttttGTCCGCCAGCGCCGGGGGGTGCGGGGGGCGGCGCCTCTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                                          <  1:1736961/147‑1 (MQ=255)
                                                                                             ctaTTTTGTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                                          <  2:7390542/148‑1 (MQ=255)
                                                                                             cgATTTTGTTCGCCAGAGCCCGGGGGGGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACCTATACCAACATCTTaaaaaaa                                                          <  2:6465823/150‑1 (MQ=255)
                                                                                              gATTTTGTTCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         <  2:1816581/150‑1 (MQ=255)
                                                                                                    gTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                                   >  2:5518459/1‑150 (MQ=255)
                                                                                                     ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGAc                                                  >  2:7766890/1‑150 (MQ=255)
                                                                                                       cGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                                                >  2:3484795/1‑150 (MQ=255)
                                                                                                        gCCAGAGCCCGGGGGGGCAGGGGGGGGGGGGGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAGGCCTGACATTa                                              <  2:10231608/150‑1 (MQ=37)
                                                                                                         ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGACATTa                                              >  2:4310646/1‑150 (MQ=255)
                                                                                                          cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                             >  2:10556544/1‑150 (MQ=255)
                                                                                                               gggggggggTGGTGGGGGGGGGGGGGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                         <  1:545521/146‑1 (MQ=255)
                                                                                                                cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTGAAAAAAAGGCCTGACATTACGCCAg                                        >  2:8865797/1‑150 (MQ=255)
                                                                                                                  gggggggGCTGGGGGGGGGGGCTTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGGTGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc                                      <  1:10108372/149‑1 (MQ=255)
                                                                                                                   gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                                     >  2:5866058/1‑150 (MQ=255)
                                                                                                                       gTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                                 >  1:78883/1‑150 (MQ=255)
                                                                                                                        tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                                >  2:8375122/1‑150 (MQ=255)
                                                                                                                        tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAGCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                                >  2:8094740/1‑150 (MQ=255)
                                                                                                                            cgggggCGGCGGCGTTGGGCCGCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                             <  2:8239866/149‑1 (MQ=255)
                                                                                                                             ggggggggggGCGTTCGGCGCCCCCGCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTTATGTCAGGCCtttttttaa                            <  1:3461069/150‑8 (MQ=255)
                                                                                                                               gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                          >  2:7575608/1‑150 (MQ=255)
                                                                                                                               gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                          >  1:2690104/1‑150 (MQ=255)
                                                                                                                               gggCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                          >  1:10335206/1‑150 (MQ=255)
                                                                                                                                ggcgggggcGTGGGGCCCCCCCGCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                         <  2:5326225/150‑1 (MQ=255)
                                                                                                                                ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                         >  2:11134737/1‑150 (MQ=255)
                                                                                                                                  cggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta                       >  2:1680553/1‑150 (MQ=255)
                                                                                                                                    gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                     >  1:659346/1‑150 (MQ=255)
                                                                                                                                    gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                     >  2:5991718/1‑150 (MQ=255)
                                                                                                                                    gcggcgATCGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGATCATGAAAGGAGCGCAACAATATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                     <  2:232909/150‑1 (MQ=255)
                                                                                                                                      ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                   >  2:6231951/1‑150 (MQ=255)
                                                                                                                                        ggttgggCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                 <  1:6507982/145‑1 (MQ=255)
                                                                                                                                        cgttgggcggcCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACATCATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                 <  2:6653008/145‑1 (MQ=255)
                                                                                                                                         tttgggCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                <  2:8754039/146‑1 (MQ=255)
                                                                                                                                         gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                >  1:5335253/1‑150 (MQ=255)
                                                                                                                                          aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt               >  2:2170603/1‑150 (MQ=255)
                                                                                                                                          aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt               >  2:10369481/1‑150 (MQ=255)
                                                                                                                                          aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt               <  1:7410195/150‑1 (MQ=255)
                                                                                                                                           ttGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt              >  2:3567603/1‑150 (MQ=255)
                                                                                                                                              gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           >  2:278009/1‑150 (MQ=255)
                                                                                                                                              gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           >  2:8658491/1‑150 (MQ=255)
                                                                                                                                              gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           >  1:2983431/1‑150 (MQ=255)
                                                                                                                                                  ccccccGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagc       <  1:9671980/150‑1 (MQ=255)
                                                                                                                                                   ccccGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  1:2263268/150‑1 (MQ=255)
                                                                                                                                                    cccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag     >  1:3889151/1‑150 (MQ=255)
                                                                                                                                                       tGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt  >  2:9309146/1‑150 (MQ=255)
                                                                                                                                                       tGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt  >  1:620341/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGT  >  NC_000913/915077‑915375

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: