Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 12.4% 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,2260TG12.4% 627.4 / 16.3 340intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (167/73);  new base G (17/28);  total (225/115)
Fisher's exact test for biased strand distribution p-value = 7.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGT  >  NC_000913/915077‑915375
                                                                                                                                                     |                                                                                                                                                     
caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                                       <  2:535168/150‑1 (MQ=32)
 aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                                      >  1:6596923/1‑150 (MQ=32)
 aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                                      <  2:9622822/150‑1 (MQ=32)
   ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAgcg                                                                                                                                                    >  1:7814652/3‑150 (MQ=17)
    gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc                                                                                                                                                   <  2:2098387/150‑1 (MQ=32)
       gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGgcgcgcg                                                                                                                                                >  1:4891324/1‑149 (MQ=17)
          ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCt                                                                                                                                             >  1:5248732/1‑150 (MQ=17)
           gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCtt                                                                                                                                            >  2:10811849/1‑150 (MQ=32)
           gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCtt                                                                                                                                            >  2:92308/1‑150 (MQ=17)
            cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCAGCGCGCTCCATg                                                                                                                                           >  2:4926544/1‑150 (MQ=12)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  2:8133631/1‑150 (MQ=32)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  2:150304/1‑150 (MQ=32)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  2:7413529/1‑150 (MQ=32)
                  gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAg                                                                                                                                     >  1:9366437/1‑150 (MQ=17)
                   cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                                    >  2:5772712/1‑150 (MQ=32)
                   cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCATGCGCCAGt                                                                                                                                    >  1:6716886/1‑150 (MQ=17)
                      gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                                 >  2:6887169/1‑150 (MQ=33)
                       ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCACTGGCa                                                                                                                                >  2:3090216/1‑150 (MQ=21)
                        tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                                               >  1:3946478/1‑150 (MQ=35)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTCAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                          >  2:6168943/1‑150 (MQ=255)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGTGGCAATATGt                                                                                                                          >  1:4696344/1‑150 (MQ=255)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                          >  1:7345177/1‑150 (MQ=39)
                               ggttGTTCCCGCGCGCCCGCCCGGGTCGCCTTATGGGATACGTCGCGGGCCCCGCGCCGGGGCTTTTTGGGGCGCCCGCGCCGGGGGGGGCCGGGGGGGGCGGCGTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                        <  2:8885620/147‑1 (MQ=9)
                                    ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                   >  1:10422008/1‑150 (MQ=255)
                                    ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                   >  1:19131/1‑150 (MQ=255)
                                      ctgtggcgcCCGCCGGGTTAGCCTTATGTGATACGTGGGGGGCCCCGCGCCGCGGGTTTTTTGTGCCCCTGCGCCGGGGGGGGCGGGGGGGGGGGGGGTTTGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                                 <  2:10575096/145‑1 (MQ=255)
                                      aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCCATATGTTGCTTAGCt                                                                                                                 >  1:2548138/1‑150 (MQ=255)
                                         tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCag                                                                                                              >  1:10575096/1‑149 (MQ=255)
                                         tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCCCCCCCGGCCCGCTCCATTCGCCAGTGGCAAAATGTTTCTTAGCTCAt                                                                                                              >  1:1975711/1‑150 (MQ=14)
                                          cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                             >  2:1639100/1‑150 (MQ=255)
                                          cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGAGCCAGTGGCAATATTTTTCTTATCTAATg                                                                                                             >  1:1634590/1‑150 (MQ=14)
                                             gcCCGACCCGGTAGCCGTATGGGTTACGTCGTGGGCACGGCGCCGCGGCTTTTTTGTTGGCCAGAGCCGGGGGGTGCTGGGGGCGGGGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaaa                                                                                                          <  2:1854774/150‑1 (MQ=17)
                                               ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCCCTGGCAAAATGTTGCTTAGCGCATGAAAgg                                                                                                        >  1:1020902/1‑150 (MQ=14)
                                                 gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                                      >  2:2657044/1‑150 (MQ=255)
                                                 gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCCGTGGCAAAATGTTGCTTAGCTCATGAAAGGAg                                                                                                      >  1:3442015/1‑150 (MQ=255)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAaatta                                                                                         >  1:9493442/1‑148 (MQ=17)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCACCGCCCTCCATGCGCCACTGGCAAAAAGTTGCTTAGCTCAAGAAAGGAGCCcaacaacaaga                                                                                         >  2:2825580/1‑145 (MQ=17)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAagaaga                                                                                         >  2:3771229/1‑150 (MQ=255)
                                                                  ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGGCCCCCCCCGCGCGCGCCGTGCGCCACTGGCAAAATTTTTAGTAGGACCAGAAAGGGGCGCAACAAGATGATGaa                                                                                     >  1:4517079/1‑150 (MQ=14)
                                                                  ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATTATGaa                                                                                     >  2:1027390/1‑150 (MQ=255)
                                                                   taTAGGTCGTGGGCACGGGGCCGGGGGTTTTTTGTGGGCCTGCGCCGGGGGGTGCTGGGGGGGGGGGCGTTGGGCGGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                                    <  2:704262/149‑1 (MQ=16)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAACTTGATGTAt                                                                                    >  2:2958672/1‑150 (MQ=18)
                                                                        gtggtgggCCCCGGGCCGGGGGGTTTTGGTGCGCCTGTGCGGGGGGGTGCTGGGGGGGGCGGCGTTGGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                               <  2:9054201/145‑1 (MQ=18)
                                                                        gtggtgggCCCCGCGCTGGGGTTTTTTTGTGCGCATGAGCGGGGGGGTGGTGGGGGGGGGGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                               <  2:4326742/145‑1 (MQ=18)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCCGGGGGCGGCGGCGGGTGGCCGCCCCCGGCGCGCGCCCTGCGCCAGGTGCCATATGTTTCTTTTCTCATGACAGGAGCGCAACAAGAAGATTAATAAACaa                                                                              >  1:5392348/1‑149 (MQ=17)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                              >  2:10067517/1‑150 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                              >  2:3606542/1‑150 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACaa                                                                              >  2:5515217/1‑149 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                              >  2:3158989/1‑150 (MQ=255)
                                                                          cGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACata                                                                             >  2:8347085/1‑150 (MQ=255)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGATACCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATTTTTCTTAGCTCATGAAAGGAGCGCAACAATATGATGAATAAAAatat                                                                            >  2:5177678/1‑150 (MQ=18)
                                                                             gggCCCCGCGCTGGGGTTTTTTTGTGGGCCGGCGCGGGGGGGTGCGGGGGGGGGCGGCGTGGGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                          <  2:9392989/150‑1 (MQ=255)
                                                                              ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaac                                                                         >  2:1385515/1‑150 (MQ=255)
                                                                                         gcggcgATTTTGTTCGCCTGCGCCGGGGGGTGCTGGGGGCGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                              <  2:10010350/150‑1 (MQ=25)
                                                                                            gctaTTTTTTTCGCCTGTGCCGGGGGGTGCTGGGGGGGGCGGCGTTTGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                           <  2:4574108/147‑1 (MQ=25)
                                                                                            gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGTGGCAATATTTTGCTTATATCATGAAAGGAGCGCAAAAAAATTAATAAACAACAAAAAACAACAACATaaaaaa                                                           >  1:6149478/1‑150 (MQ=17)
                                                                                              gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         >  2:8754976/1‑150 (MQ=255)
                                                                                              gATTTTGGCCGCCCGAGCCGGGGGGGGCTGGGGGCGGCGGCTTTTGGCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCGATATGTTGCCTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         <  2:4476182/150‑1 (MQ=255)
                                                                                                 tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                                      >  2:7608000/1‑150 (MQ=255)
                                                                                                 tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCTTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                                      >  2:4720912/1‑150 (MQ=255)
                                                                                                        gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                               >  2:10944860/1‑150 (MQ=255)
                                                                                                            gggCCGGGGGGTGCTGGGGGGGGGGGGGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                                           <  2:7754645/148‑1 (MQ=37)
                                                                                                             tgcgggggggTGATGGGGGCGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                                          <  2:819926/145‑1 (MQ=38)
                                                                                                              gCCGGGGGGTGGTGGGGGGGGGGGCTTTGGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                         <  2:6661007/150‑1 (MQ=255)
                                                                                                               gggggggtTGCGGGGGGGGGGGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGAAATATGTTGCTTTGCTCATGAAAGGAGCGCAACAAGATTATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                                        <  2:7763144/147‑1 (MQ=37)
                                                                                                                ccGGTGGAGCAGGGCGAGGCGGAGATTGTACGCCCCCAGCGCGCTACATGCGCCAGTGGCAATATTTTGAATATCACATTAAAGGACCGAAAAAAGATGATGAATAAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                       >  2:7494015/1‑150 (MQ=255)
                                                                                                                 gggggtTGCTGGGGGCGCCGGGGATTGGCCGCCCACTGCGCGATCCTTTCGACAGTGGCAATATGTTGATTATCTCTTGAAATGAGCGCAACAAGATTATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACTTTACGCCAGGc                                      <  2:7867790/149‑1 (MQ=255)
                                                                                                                  gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                                     >  2:6305810/1‑150 (MQ=255)
                                                                                                                    ggggGCGGGGGGGGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                                   <  2:3442015/150‑1 (MQ=255)
                                                                                                                      gTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGAAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAATATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                                 <  2:10063682/150‑1 (MQ=255)
                                                                                                                         gtggggggggggGCTTTTGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                              <  2:10422008/148‑1 (MQ=255)
                                                                                                                           gggggCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACACATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                            <  2:9843895/150‑1 (MQ=255)
                                                                                                                           gggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                            <  2:9938922/150‑1 (MQ=255)
                                                                                                                            ggggggggggCGTTGGCCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                           <  2:3256820/150‑1 (MQ=255)
                                                                                                                            ggggcgggggcgCTGGGCCGCCCCGCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGACAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                           <  2:8722109/150‑1 (MQ=38)
                                                                                                                            cagggggcggccAATCGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                           >  2:2324143/6‑150 (MQ=255)
                                                                                                                             gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                          >  1:3570819/1‑150 (MQ=255)
                                                                                                                             agggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                          >  1:9938266/5‑150 (MQ=39)
                                                                                                                                   cggcgATGCCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc                    <  1:562307/150‑1 (MQ=255)
                                                                                                                                      atcgccgccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                 >  1:10033189/7‑150 (MQ=255)
                                                                                                                                       gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                >  2:4094354/1‑150 (MQ=255)
                                                                                                                                        aTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt               >  1:2227724/1‑150 (MQ=255)
                                                                                                                                         ttGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt              >  2:9296867/1‑150 (MQ=255)
                                                                                                                                          gggCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt             <  1:2657044/149‑1 (MQ=255)
                                                                                                                                            gcccccCCCGCCGCGCGCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           <  1:4130097/150‑1 (MQ=255)
                                                                                                                                             ccgccCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCCCGCTTTTcc          <  2:6952594/150‑1 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  2:1989546/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  2:7128075/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  2:6388463/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  1:8702472/1‑150 (MQ=255)
                                                                                                                                             ccccccccTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          <  2:4827583/147‑1 (MQ=255)
                                                                                                                                             ccccccccGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          <  1:8792049/147‑1 (MQ=255)
                                                                                                                                              cgccCCCCCCGCGCTCCTTGCGCCAGTCGCAATATGTAGCTTAGCACATGAAAGGACCGCAACAAGATCATGAAACAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca         <  1:10067517/150‑1 (MQ=255)
                                                                                                                                              cccccccGGCGCGCTCCTTGCGCCAGTGGCACGATGTTGCTGAGCTCATGAAAGGAGAGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca         <  1:10882805/148‑1 (MQ=255)
                                                                                                                                               gccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag        >  1:766190/1‑150 (MQ=255)
                                                                                                                                               gccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag        >  1:9071082/1‑150 (MQ=255)
                                                                                                                                                 ccccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  2:10150442/150‑1 (MQ=255)
                                                                                                                                                 ccccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  2:4865271/150‑1 (MQ=255)
                                                                                                                                                  cccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag     <  2:10033189/150‑1 (MQ=255)
                                                                                                                                                     tGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt  >  2:8061787/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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: