Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2130AC26.4% 263.2 / 36.7 215intergenic (‑356/+139)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base A (118/24);  new base C (51/0);  total (179/36)
Fisher's exact test for biased strand distribution p-value = 7.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915065‑915363
                                                                                                                                                     |                                                                                                                                                           
cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGAtt                                                                                                                                                           >  8:1897680/1‑151 (MQ=14)
cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                           <  1:1281850/151‑1 (MQ=31)
 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTg                                                                                                                                                          >  6:1393593/1‑151 (MQ=14)
 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                          <  7:1535916/151‑1 (MQ=31)
   ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑‑c                                                                                                                                                      <  3:451176/151‑1 (MQ=31)
    gCAGCCCAGATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑‑cc                                                                                                                                                     <  1:820848/151‑1 (MQ=14)
    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑‑cc                                                                                                                                                     <  4:289682/151‑1 (MQ=31)
    gCAGCCCACAAGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑‑cc                                                                                                                                                     <  1:885415/151‑1 (MQ=14)
     cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGGGCTTGg‑‑ccc                                                                                                                                                    >  6:1512996/1‑150 (MQ=9)
      aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑‑ccgc                                                                                                                                                   <  4:1379788/151‑1 (MQ=31)
          cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGccccc                                                                                                                                               <  2:893476/151‑1 (MQ=31)
          aacatgcattagcgCTGAGGGCGCGTTTTAAAGGCATGTTACCTAGAGCCCGACACGGTAGCCGTAAGGGATAAGTCGAGGGCAACGCCCACCCGCGAATTTCTTCGCCAGAGCACGCAGGTGCACGGGGCGGCGGCGATTGG‑‑CCGccccc                                                                                                                                               >  8:1548819/14‑151 (MQ=2)
            caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑g                                                                                                                                           <  8:22819/151‑1 (MQ=31)
            caTGGAGGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑g                                                                                                                                           <  1:725675/151‑1 (MQ=14)
                 aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGGGCTTGG‑‑CCGCCCCCT‑‑gcgcgc                                                                                                                                      >  8:879133/1‑151 (MQ=14)
                   gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTc                                                                                                                                    >  2:894900/1‑151 (MQ=31)
                   gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGCTTGG‑‑GCGCCCCCG‑‑gcgcgcgc                                                                                                                                    >  6:588393/1‑149 (MQ=255)
                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTg                                                                                                                               >  1:12304/1‑151 (MQ=31)
                        cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTgc                                                                                                                              >  8:464413/1‑151 (MQ=14)
                        cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTgc                                                                                                                              >  3:1138187/1‑151 (MQ=31)
                         tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTgcg                                                                                                                             >  6:1128921/1‑151 (MQ=31)
                             ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGcc                                                                                                                           >  5:1078087/1‑149 (MQ=31)
                             ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGAGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAg                                                                                                                         <  5:776809/151‑1 (MQ=14)
                              gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGt                                                                                                                        >  4:724827/1‑151 (MQ=14)
                                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTgg                                                                                                                      >  2:1797267/1‑151 (MQ=32)
                                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCGTGCGCCAGTgg                                                                                                                      >  8:489926/1‑151 (MQ=17)
                                  gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGGGGGGCGGTG‑‑GCGCCCCCC‑‑GCGCCCGCCGTgcggccgggctt                                                                                                                    >  7:1159853/1‑142 (MQ=11)
                                    tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCCCCCCCG‑‑GCGCGCGCCGTGCGCCCGTGGCCat                                                                                                                  >  7:890565/1‑151 (MQ=11)
                                      tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTATGGGCTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCGTTGCGCCAGTGGCAatat                                                                                                                >  6:1154147/1‑151 (MQ=255)
                                      tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAatat                                                                                                                >  2:1539263/1‑151 (MQ=37)
                                        cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                              >  3:941779/1‑151 (MQ=38)
                                          gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGGTGG‑‑CCGCCCCCC‑‑GCGCGCGCCGTGCGcccgggggcaaaagttg                                                                                                            >  6:1770656/1‑136 (MQ=11)
                                              tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                                                         >  8:1827666/1‑150 (MQ=255)
                                               gTTACCTCGCGCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCGTGG‑‑CCGCCCCCG‑‑GCGCGCGCCGTGCGCCAGGGGCAATATGTTGCtt                                                                                                         >  5:1825007/1‑149 (MQ=14)
                                                  aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCCGGGGCCGCCGCGATTGG‑‑CCGCCCCCG‑‑GCGCGCTGCATGCGCCAGGGGCCATAAGTTGGTTAGCTc                                                                                                    >  4:904910/1‑151 (MQ=14)
                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCCGCGGCGCTTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCCATATGTTGGTTAGCTCAt                                                                                                  >  4:1569537/1‑151 (MQ=255)
                                                         gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTTG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAgggcaagaa                                                                                               >  1:200663/1‑141 (MQ=14)
                                                         gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCGCCTTGCGCCAGTGGCAATATGTTGGGTAGCGCCTGaa                                                                                               >  1:654127/1‑149 (MQ=255)
                                                           ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCCGGGGGCGCGGGGCTTTG‑‑CCGCCCCCC‑‑GCGCCCCCCGTTGGGCCGGGGGAAATTTTTTGCttgcggaatcaaggc                                                                                            >  7:1674631/1‑138 (MQ=2)
                                                           ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaaacg                                                                                            >  3:1663429/1‑148 (MQ=255)
                                                            cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCT‑‑GCGCGCTCCATGCGCCAGTGGCAATATGTTGGTTAGGTCATGAAAggtg                                                                                          >  5:285086/1‑149 (MQ=255)
                                                                 cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTG‑‑CAGCCCCCT‑‑GCGCGCACCTTGAGCCAGAGGCAAAATGTTGCATACCTCATGAAAGCAGCGca                                                                                      >  2:967165/1‑150 (MQ=17)
                                                                 cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCAAGAAAGGAGCGca                                                                                      >  3:407409/1‑150 (MQ=255)
                                                                     aGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCCCCAGTGGCAAAATGTTGCTTAGCTGATGAAAGGAGCGcaacaa                                                                                  >  2:1419083/1‑150 (MQ=255)
                                                                       ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCT‑‑GCGCGCGCCGTGCGCCACTGGCAATATGTTGCTTAGCGCAAGCAAGGAGCGcaaca                                                                                   >  6:202762/1‑147 (MQ=17)
                                                                       ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCA‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaaca                                                                                   >  3:954420/1‑147 (MQ=255)
                                                                       ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCGCTTGG‑‑CCCCCCCCT‑‑GCGCCCTCCATGCCCCAGTGGCCATATGTTGGGTAGCTGATGAAAGGAGCGccaca                                                                                   >  8:1604746/1‑147 (MQ=17)
                                                                        cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAAGGGGCCGCGGCGAATGG‑‑CCGCCCCCT‑‑GCCCGCTCCTTTCGCCAATGGCAATATTTTTCATATCTCATGAAAGGAGCGCAACAAgatg                                                                              >  5:556043/1‑151 (MQ=17)
                                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGCATGG‑‑CCCCCCCCA‑‑GCCCGCGCCGTGCGCCCGGGGCAATAATTTGCGTAGCTCAAGAAACGAGCGCCACAAga                                                                                >  8:1726884/1‑148 (MQ=14)
                                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCCGCCGCGATTTG‑‑CCCCCCCCA‑‑GCGCGCTCCATTCGCCAATGGCAATATTTTTCATAAATAAATAAAAGAGCGcaacaaaaat                                                                              >  8:648372/1‑146 (MQ=9)
                                                                            aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGatgatg                                                                           >  4:1025965/1‑150 (MQ=255)
                                                                            aGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCCACAAgatga                                                                             >  7:1115706/1‑148 (MQ=255)
                                                                             gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGG‑‑GCGCCCCCT‑‑GCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAtgatga                                                                          >  7:441463/1‑150 (MQ=255)
                                                                              ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCCGGGGCGGGGGCGCGTGG‑‑CCGCCCCCC‑‑GCGCGCGCCCTGCGCCCGGGGCCATATTTTGGGTTGGACATGAAATGTGCGccacaacctgttcaa                                                                         >  4:1216007/1‑141 (MQ=2)
                                                                                      ggtgggCACGGCGCAGGGGCGTTTTTGTTCGCCAGAGCCGGGGGGGCTGGGGGCGGGGGCTTTTGG‑‑CCGCCCCGT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                <  3:1395908/147‑1 (MQ=21)
                                                                                      cGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCA‑‑GCGCGCTCCTTGCGCCAGTGGCAATATTTTTCTTATCTCATGAAAGGAGCGCAACAAGATGATGAATAAACatat                                                                >  3:1812184/1‑151 (MQ=18)
                                                                                            cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCT‑‑GCGCGCGCCTTGCGGCAGTGGCAATATTTTGCGTAGCTCATGAAAGGAGCGCCACAAGATGCTGAATCAACATATaacaa                                                           >  1:1876822/1‑150 (MQ=21)
                                                                                            cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCGCTTGG‑‑CCCCCCCCG‑‑GCGCGCGCCGTGCGCCACTGGCAATATGTTGCTTAGCGGATTAAAGGAGCGCCACAAGATGATGAATCAACATATaacaac                                                          >  4:804364/1‑151 (MQ=21)
                                                                                             aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGG‑‑CCGCCCCCG‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCGTAGCGCATGCAAGGGGCGGCACAAGATGATGAATCACCATATaacaac                                                          >  7:1273231/1‑150 (MQ=18)
                                                                                              ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGCGCGCAACAAGATGATGAATCAACATATAaccaca                                                         >  1:137626/1‑150 (MQ=255)
                                                                                               cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCCGCGGCGCTTGG‑‑CCCCCCCCC‑‑GCGCGCGCCGTGCGCCCGGGGGAAAATTTTGCTTAGCTGATGAAATGAGCGCAACAAGAAGATGAATCAACAAATAAAAACATc                                                       >  1:1647166/1‑151 (MQ=18)
                                                                                               cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCCGGGGCGGCGGGGCTTGG‑‑GCGCCCCCG‑‑GCGCGGGCCGTGCGCCCCGGGCAATATTTTTGTTTGGGCATGAAAGGGGCGGAACAAAAATATGTATGAACATATAaaaacaa                                                        >  4:266737/1‑149 (MQ=9)
                                                                                               cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                        <  6:1746169/151‑1 (MQ=255)
                                                                                                gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGGTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAAAAACATc                                                       >  5:372545/1‑150 (MQ=255)
                                                                                                   gggggtttaTTTTGATGGCCGGGGCGGGGGGGTGTGGGGGGCGTCTGCTTTGCG‑‑CCCCCCTTG‑‑GCGCGCTGCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                                   <  4:988008/143‑1 (MQ=255)
                                                                                                   cAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTa                                                    >  2:1152882/1‑150 (MQ=255)
                                                                                                      cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTTG‑‑GCGCCCCCA‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGCGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                 >  7:987591/1‑150 (MQ=255)
                                                                                                      cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                                 >  3:203835/1‑150 (MQ=255)
                                                                                                       ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCGCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAACTTaaaaa                                                >  6:773839/1‑150 (MQ=255)
                                                                                                       ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCCCCCCCT‑‑GCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTAATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                >  2:182276/1‑150 (MQ=255)
                                                                                                         cgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑‑CCGCCCCCG‑‑GCGCGCGCCTTGCGCCAGTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAAATATAAAAACATCTTaaaaaaac                                             >  1:1056342/1‑150 (MQ=25)
                                                                                                         cgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGAATGG‑‑CCCCCCCCA‑‑GCCCGCGCCGTGCGCCAGTTGCAATATTTTTGTTAGCGGAAGAAAAGGGCGCCACAATATAATGAATAAAAAAATAACAACCAATTaaaaaaa                                              >  3:293769/1‑150 (MQ=9)
                                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaacg                                            >  3:1450936/1‑149 (MQ=255)
                                                                                                           aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                           >  5:349825/1‑151 (MQ=255)
                                                                                                                   cGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAGATCTTAAAAAAAGGCCTGACAtt                                   >  2:1689171/1‑151 (MQ=255)
                                                                                                                     gggctgCCGGGGGGGGCGGGGGGGGGGGGCGTGTGGGCCCCCCCCCC‑‑CCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                   <  4:355091/146‑1 (MQ=255)
                                                                                                                     ccAGAGCCCGGGGGTGCAGGGGGCGGCGGGGATTGG‑‑CCGCCCCCT‑‑GCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGCAAGGAGCGCAACCAGATGATGAATCAACATATAACAAAATCTTAAAAAAAGGCCGGACATTAc                                 >  1:1815177/1‑151 (MQ=37)
                                                                                                                     ccAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGG‑‑CCCCCCCCA‑‑GCGCGCTACTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAATGCGCGCAACAAGATGATGAATAAACATATAACAACAACATAAAAAAACGCCTGACATTa                                  >  2:451458/1‑150 (MQ=255)
                                                                                                                         agCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATATTAAAAAAAGGCCTGACATTACGcc                              >  6:1034208/1‑150 (MQ=255)
                                                                                                                           cccGGGGGTGCAGGGGGCGGCGGCGGTTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATTTTTCTTAGCTCATGAAAGGAGGGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGCCTGCAATTACGCCAgg                           >  8:1228779/1‑151 (MQ=37)
                                                                                                                             cGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGGGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGCATCAACATATTACAACCTCTTTCAAAAAGGCGGGCCCTTACGGCAGGcc                         >  5:1497008/1‑151 (MQ=21)
                                                                                                                              gggggTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCt                        >  8:1164419/1‑151 (MQ=255)
                                                                                                                                  ggggtgaggggggGCGGGGATTGG‑‑GGCCCCCCGGGGCGCGCTCCTTGCGGCAGTTGCAATATGTTGCTTATCTCATTATAGGAGCTCAACAAGATGATTAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                       <  6:42734/144‑1 (MQ=255)
                                                                                                                                  gTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                    >  5:435272/1‑151 (MQ=255)
                                                                                                                                   tGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                   >  1:1409397/1‑151 (MQ=255)
                                                                                                                                   tGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                   >  6:1637784/1‑151 (MQ=255)
                                                                                                                                     cAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGAAGAAACAACATATAACAACATATTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGcat                 >  5:343293/1‑149 (MQ=255)
                                                                                                                                      aGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                >  8:783512/1‑151 (MQ=255)
                                                                                                                                             cggcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta           >  5:1676254/1‑150 (MQ=255)
                                                                                                                                             cggcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta           >  5:716194/1‑150 (MQ=255)
                                                                                                                                             cggcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta           >  2:1213454/1‑150 (MQ=255)
                                                                                                                                             cggcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCGGCCATTACGCCAGGCCTTCTGCGtgaatta           >  3:1249514/1‑150 (MQ=255)
                                                                                                                                              ggcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACGAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat         >  1:380592/1‑151 (MQ=255)
                                                                                                                                               gcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc        >  1:272247/1‑151 (MQ=255)
                                                                                                                                               gcggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc        >  7:133957/1‑151 (MQ=255)
                                                                                                                                                cggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc        >  7:271559/1‑150 (MQ=255)
                                                                                                                                                cggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc        >  5:1251479/1‑150 (MQ=255)
                                                                                                                                                 ggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCATGCCTTCTGCGTTAATTAATCAc      >  2:1000292/1‑151 (MQ=255)
                                                                                                                                                 ggcgATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc      >  3:1756374/1‑151 (MQ=255)
                                                                                                                                                 ggcgAGTGC‑‑CCGCCCCCG‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc      <  3:679029/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                           
CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑‑CCGCCCCCT‑‑GCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTT  >  NC_000913/915065‑915363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: