Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG43.3% 170.5 / 41.5 104T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (7/52);  new base G (45/0);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 4.31e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366377‑366673
                                                                                                                                                       |                                                                                                                                                    
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGggtg                                                                                                                                                     <  1:980375/151‑1 (MQ=255)
aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGggt                                                                                                                                                      <  2:1754929/150‑1 (MQ=255)
 aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGggtgg                                                                                                                                                    <  2:45510/151‑1 (MQ=255)
   gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGtt                                                                                                                                                  <  2:572218/151‑1 (MQ=255)
    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGttt                                                                                                                                                 <  2:23751/151‑1 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGATTTGGGG‑GCCGGGGGGGtttt                                                                                                                                                >  1:1119441/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGgtttt                                                                                                                                               >  1:945358/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGgtttt                                                                                                                                               >  2:1884119/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCCGGGGGGgtttt                                                                                                                                               >  2:860043/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGAGTA‑TTTGGC‑GGCAGGGGGGttttt                                                                                                                                               >  1:336725/1‑151 (MQ=255)
      tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGG‑TTTGGC‑GGCCGGGGGGgtttt                                                                                                                                               >  1:434048/1‑151 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGttttt                                                                                                                                               <  2:609695/150‑1 (MQ=255)
       ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGGG‑TTTGGG‑GGCAGGGGGGgttttt                                                                                                                                              >  1:489391/1‑150 (MQ=255)
        gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGCGTA‑TTTGGG‑GCCCGGGGGGgtttt                                                                                                                                               >  1:785482/1‑149 (MQ=255)
          gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTttttcttt                                                                                                                                           <  1:1426544/151‑1 (MQ=255)
            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGG‑GCCCGGGGGGGttttttttt                                                                                                                                          >  2:40914/1‑150 (MQ=255)
             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGGttttttttt                                                                                                                                          >  1:1305183/1‑149 (MQ=255)
              cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGT‑TTTGGG‑GCCCGGGGGGGgtttttttt                                                                                                                                          >  1:39247/1‑148 (MQ=255)
               tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGG‑GCCCGGGGGGTTTTTTTTTTTAcc                                                                                                                                      >  1:1030114/1‑151 (MQ=255)
                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCAGGGGGGGtttttttttttcc                                                                                                                                      >  2:1591628/1‑146 (MQ=255)
                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                <  2:715450/151‑1 (MQ=255)
                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                               <  2:960219/150‑1 (MQ=255)
                          cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGTTTTTTTTTTCACCAGTGAGACggg                                                                                                                           >  2:232652/1‑151 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTCCCAGTGAGACggg                                                                                                                           >  2:1233325/1‑150 (MQ=255)
                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTCACCAGGGAGACggg                                                                                                                           >  2:641544/1‑150 (MQ=255)
                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                                                                                                           <  2:667774/149‑1 (MQ=255)
                               tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGGGGC‑TTGGGC‑GCCAGGGGGGttttttttttttcccgatggacg                                                                                                                             >  1:1789491/1‑131 (MQ=255)
                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                    <  1:1222771/151‑1 (MQ=255)
                                  caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                                   <  1:1717351/151‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                               <  1:537073/151‑1 (MQ=255)
                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                                <  2:912081/150‑1 (MQ=255)
                                       aTTGCGTTGCGCTCACTACCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                              <  1:155205/151‑1 (MQ=255)
                                        ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGGCCGGGGGGGTTTTTTTTTTC‑CCCGTGAGACGGGCAACAACTGGtt                                                                                                               >  1:798247/1‑149 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGC‑TTGGGC‑GGCAGGGGGGGTTTTTTTTTTACCCGGGAGAAGGGCAACAGCTGAtt                                                                                                               >  1:825575/1‑148 (MQ=255)
                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGcc                                                                                                            >  2:548019/1‑151 (MQ=255)
                                          gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAACCGATTTccc                                                                                                           >  2:1889871/1‑151 (MQ=255)
                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAAAAGATGATTTccc                                                                                                           >  2:452190/1‑150 (MQ=255)
                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTTCCCAAGGAGACGGGCAAAAGCTGATTGCCCtt                                                                                                         >  2:1859634/1‑149 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                                                                                                       <  1:1644625/150‑1 (MQ=255)
                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                                                                                                       <  2:1875981/150‑1 (MQ=255)
                                                 gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                                                                                                    <  1:748279/151‑1 (MQ=255)
                                                  cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGTGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGTTTTTTTTTTTACCAGGGAGTCGGGGAACAGGTGATTGCCCTTTACCGc                                                                                                   >  2:686368/1‑151 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGA‑TTTGGG‑GCCAGNNNNNTNTTTTNTNTTGNCCGGGAGACGGGCAACAGGTGGGTGGCCTTCCCCGcc                                                                                                  >  1:1322683/1‑150 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTT‑TTTGGC‑GGCCGGGGGGGTTTTTTTTTCACCCGGGCGACGGGCAAAAGCTGCTTTTCCTTCAcc                                                                                                     >  2:642721/1‑147 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GGCCGGGGGGGTTTTTTTTTCACCCGTGAGACGGGCGACAGCTGATTGCCCTTCACcccc                                                                                                  >  1:323926/1‑150 (MQ=255)
                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                                  <  2:963540/150‑1 (MQ=255)
                                                      cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                <  1:251142/150‑1 (MQ=255)
                                                       tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGG‑TTTGGT‑GCCCGGGGGGGTTTTCTTTTTCCCAGGGGGGAGGGCAACAGCTGATTTCCCTTTCCCGCCCgg                                                                                               >  2:1586328/1‑150 (MQ=255)
                                                           cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑GCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGGAACACCTGAGTGTCCTTTACCCCCTGGgcc                                                                                            >  1:780726/1‑149 (MQ=255)
                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGGCCTTCACCGCCCGGccc                                                                                            >  1:162123/1‑148 (MQ=255)
                                                                  cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCTGagagag                                                                                    >  1:1553849/1‑150 (MQ=255)
                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                  <  1:1507966/151‑1 (MQ=255)
                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                               <  2:1832668/151‑1 (MQ=255)
                                                                        ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                             <  2:285115/151‑1 (MQ=255)
                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                           <  2:1613608/151‑1 (MQ=255)
                                                                          aaaCCTGTCGTGCCAGCTACATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                           <  2:1524725/151‑1 (MQ=255)
                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                          <  2:1305390/151‑1 (MQ=255)
                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                          <  1:1265147/151‑1 (MQ=255)
                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGG‑GCCCGGGGGGGTTTTCTTTTTACCCCTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTGCAGCaa                                                                           >  1:1747152/1‑150 (MQ=255)
                                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                         <  2:1221111/151‑1 (MQ=255)
                                                                              cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGTA‑TTTGGC‑GGCCGGGGGGGTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCCTCCCCCCCCGGCCCTCAGGGAGGTGCCGCAAGCgg                                                                       >  1:530080/1‑151 (MQ=255)
                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                       <  1:635945/149‑1 (MQ=255)
                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                      <  2:1623509/150‑1 (MQ=255)
                                                                                  cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                   <  1:821616/151‑1 (MQ=255)
                                                                                   gcgctaGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCAc                                                                  <  1:586586/146‑1 (MQ=255)
                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGt                                                           <  2:448373/151‑1 (MQ=255)
                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGGTTTTTTTTTCACCAACGAGACGGGCCACAGCTGACTGCCCTTCACCCCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACG‑Ctcggtt                                                         >  1:1515316/1‑144 (MQ=255)
                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGtt                                                          <  2:527757/148‑1 (MQ=255)
                                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCcag                                                  <  2:258853/151‑1 (MQ=255)
                                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGGGGGTTTTCTTTTTACCAGGGAGACGGGCAACAGCCGATTTCCCCTCACCCCCCGGCCCTGAGAGAGTTGCAGCAAACGGTCCACG‑CTGGGTTGcccc                                                    >  2:1613534/1‑149 (MQ=255)
                                                                                                    gAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGC‑GCCCGGGGGGTTTTTCTTTTTACCAGGGAGAAGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGCCGGGTTTGccc                                                     >  1:715943/1‑147 (MQ=255)
                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAgg                                              <  2:659454/151‑1 (MQ=255)
                                                                                                       tCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑GCCCGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCCGCCCTGAGAGAGTTTCAGCAAGCGGGCCACG‑CTGGTTTGCCCCAGCAgg                                              >  1:1224038/1‑150 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGc                                             <  2:1538314/150‑1 (MQ=255)
                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCg                                            <  1:52437/151‑1 (MQ=255)
                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCAGGGGGGGTTTTCTTTTTACCCGTGAGAAGGGGGACAGGTGATTGCCCTTCAACGCCCGGCCCCGGGAGAGGTGCAGCAAGCGGGCCACG‑CTGGTTTGCCCCCGCCGGCg                                            >  1:524184/1‑150 (MQ=255)
                                                                                                           ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGaaa                                         <  2:94127/151‑1 (MQ=255)
                                                                                                            cAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCAGGGGGGGTTTTTTTTTCAACCGAGAGACGGGGAACAGCCGATTGGCCCTCCCCCCCCGGCCCTGGGAGAGTTGCAGCAAACGGGGCACG‑CTGGGTTGCCCCCGCAGGCGa                                           >  1:1939164/1‑148 (MQ=255)
                                                                                                            cAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGG‑GGCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCCGATTTCCCCTCACCCCCCGGCCCTGGGAGGGGTGCAGCAAACGGGCCCCG‑CTGGGTTGCCCCCGCCAGCGaaa                                         >  1:1629594/1‑150 (MQ=255)
                                                                                                            cAACGCGCGGGGAGAGGCGGGTTTCGGA‑TTTGGC‑GCCCGGGGGGTTTTTCTTTTCACCCGGGAGACGGGCAACAGCTGGTTGCCCCTCACCGCCTGGCCCCGTGAGAGGTGCAGCAAGAGGTCCACG‑CTGGTTTTCCCCCGCAgggggaa                                         >  1:643347/1‑145 (MQ=255)
                                                                                                             aaCGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAAt                                       <  1:1866523/151‑1 (MQ=255)
                                                                                                             aaCGCGCGGGGAGAGGCGGTTTGCGCA‑TTGGGC‑GCCAGGGGGGGTTTTTTTTTTCCCAGTGGGAAGGGCAAAAGCTGATTGCCCCTCACCGCCCGGCCCTGAGGGAGGTGGAGCAAACGGGCCCCC‑CCGGGTTGCCCCACCAGGGGgaaa                                        >  1:1394248/1‑150 (MQ=255)
                                                                                                               cgcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATcc                                     <  2:1627007/151‑1 (MQ=255)
                                                                                                                gcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCt                                    <  2:61726/151‑1 (MQ=255)
                                                                                                                 cgcgGGGAGAGGCGGTTTGCGTA‑TTTGGG‑GCCCGGGGGGGTTTTTTTTTTCCCAGGGAGACGGGGAACAGCTGATTGGCCTTCCCCGCCCGGCCCCTTGAGAGTTGTCGCCAGCGGTCCCCC‑CCGGGTTTCCCCCGCAGGCGGAAAttc                                     >  1:525888/1‑147 (MQ=255)
                                                                                                                  gcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCt                                    <  2:1712662/149‑1 (MQ=255)
                                                                                                                         gagGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           >  2:316309/1‑151 (MQ=255)
                                                                                                                           ggCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                         >  1:1285232/1‑151 (MQ=255)
                                                                                                                              ggTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                      <  1:476241/151‑1 (MQ=255)
                                                                                                                               gTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                      >  2:933732/1‑150 (MQ=255)
                                                                                                                                 ttGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                    <  2:780036/150‑1 (MQ=255)
                                                                                                                                   gCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                  <  1:1175850/150‑1 (MQ=255)
                                                                                                                                          tGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa            <  2:785634/150‑1 (MQ=255)
                                                                                                                                           gggC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg             >  1:1095400/1‑148 (MQ=255)
                                                                                                                                            ggC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggggt           >  1:1601662/1‑147 (MQ=255)
                                                                                                                                            ggC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat         <  1:452118/151‑1 (MQ=255)
                                                                                                                                                gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata        >  1:1573186/1‑149 (MQ=255)
                                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg   <  2:352345/151‑1 (MQ=255)
                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa  <  1:641423/151‑1 (MQ=255)
                                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa     >  1:484488/1‑145 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑GCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACG‑CTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366377‑366673

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