Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG40.0% 130.7 / 29.2 86T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/42);  new base G (34/0);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 1.11e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAT  >  NC_000913/366378‑366671
                                                                                                                                                        |                                                                                                                                                   
aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑Gggtgg                                                                                                                                                   <  1:683903/151‑1 (MQ=255)
 aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑Gggtggt                                                                                                                                                  <  1:185652/151‑1 (MQ=255)
   ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGttt                                                                                                                                                <  1:294682/151‑1 (MQ=255)
    cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTTGGCGCCC‑GGGGGGgttt                                                                                                                                               >  1:794458/1‑151 (MQ=255)
     tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTG‑TTTGGCGCCC‑GGGGGGgtttt                                                                                                                                              >  1:669649/1‑151 (MQ=255)
       gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTttttc‑t                                                                                                                                           <  2:691257/151‑1 (MQ=255)
            ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCa                                                                                                                                      <  2:81631/151‑1 (MQ=255)
            ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCa                                                                                                                                      <  2:512756/151‑1 (MQ=255)
            ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCa                                                                                                                                      <  1:1382954/151‑1 (MQ=255)
             cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGGCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCa                                                                                                                                      <  2:396246/150‑1 (MQ=255)
               aaTGAGTGAGCTAACTCACATTACTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGA‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑GTTTCTCCa                                                                                                                                   <  2:128095/151‑1 (MQ=255)
               aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCAc                                                                                                                                     <  1:459517/149‑1 (MQ=255)
                 tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTT‑TGGGTA‑TTGGGGGAAG‑GGGGGGTgtttt‑ttttagta                                                                                                                                    >  1:845318/1‑145 (MQ=255)
                 tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGTGTA‑TTTGGGGCCA‑GGGGGGGTTTTT‑TTTTTAcc                                                                                                                                    >  1:709119/1‑148 (MQ=255)
                 tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCCTA‑TTGGGCGGCA‑GGGGGGGttttc‑ttttttc                                                                                                                                     >  2:443597/1‑144 (MQ=255)
                   agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTga                                                                                                                              <  2:576085/151‑1 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGTATTTGGGGGCCG‑GGGGGGTTTTT‑‑TTTTTACCA‑GTg                                                                                                                               >  1:1209219/1‑149 (MQ=255)
                      gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTgaga                                                                                                                            <  1:941041/150‑1 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACgg                                                                                                                         <  1:743599/151‑1 (MQ=255)
                         cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACg                                                                                                                          <  2:22929/149‑1 (MQ=255)
                          tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGTGTA‑TTGGGGGGCA‑GGGGGGGTTTTT‑TTTTTACCA‑GTGAGACggg                                                                                                                        >  1:178233/1‑150 (MQ=255)
                           aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGTAGAGGC‑GGTT‑TGCGTA‑TTGGGCGGCC‑GAGTGGGTTTTT‑TTTTTAccc                                                                                                                                   >  1:54152/1‑138 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAg                                                                                                                  <  2:1287039/151‑1 (MQ=255)
                                 caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTTAGACGGGCAACAg                                                                                                                  <  2:1092861/149‑1 (MQ=255)
                                  aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAGCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTg                                                                                                               <  2:1415365/151‑1 (MQ=255)
                                  aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGT‑TGCGTA‑TTTGGCGCCC‑GGGGGGGTTTTC‑TTTTTACCA‑GTGAGACGGGCAACAGCTg                                                                                                               >  1:805931/1‑151 (MQ=255)
                                      aTTGCGTTGCGCTCACTGGCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTTCATTAATGAATCGGCCGCCGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCGA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTg                                                                                                           <  2:49072/151‑1 (MQ=255)
                                       ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGTA‑TTTGGCGCCC‑GGGGGGGTTTTT‑TTTTCACCC‑ATGAGACGGGCAAAAGCTGAtt                                                                                                            >  1:889504/1‑149 (MQ=255)
                                       ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGGA‑TTTGGCGCCC‑GGGGGGGTTTTT‑TTTTTACCA‑GGGAGACGGGGAACAGCCGAtt                                                                                                            >  1:560058/1‑149 (MQ=255)
                                             tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCG‑TGGGTG‑TTTGGCGCCA‑GGGGGGGTTTTT‑TTTTTACCA‑ATGAGACGCGGAAcca                                                                                                                  >  2:1127/1‑135 (MQ=255)
                                                gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTTGTA‑TTTGGCGCCC‑GGGGGGGTTTTT‑TTTTTACCA‑GTGAGAGGGGGAACAACTGATTTCCCTTTAcc                                                                                                  >  1:687671/1‑150 (MQ=255)
                                                 cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGAT‑TTAGCT‑CTTGGCGCCA‑GGGGGGGTTTTCTTTTTCACC‑‑GTGAGACAGGGAACAACTAAGTGCCCCTTAcccc                                                                                                >  1:595703/1‑149 (MQ=255)
                                                   cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                               <  1:1004553/150‑1 (MQ=255)
                                                   cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                              <  2:350327/151‑1 (MQ=255)
                                                   cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                              <  2:980177/151‑1 (MQ=255)
                                                     cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTCTCGGTT‑TGGGGCCCCT‑GGGGGGGTTTT‑‑TTTTTAACA‑CGGTGACGGGCAACAACCGAATGcccactaccccccgg                                                                                            >  2:1470088/1‑139 (MQ=255)
                                                     cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGTA‑TTGGGGGCCC‑GGGGGGGTTTTC‑TTTTTACCA‑GTGAGAGGGGCAACAGACGGTTGCCCTTCACCGCCTgg                                                                                            >  1:116844/1‑151 (MQ=255)
                                                       gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTGGTA‑TTTGGGGCCC‑GGGGGGGTTTTC‑TTTTTACCA‑GTGAGACGGGGAACAGCTGATTTCCCTTCACCGCCTGGcc                                                                                          >  1:1167617/1‑151 (MQ=255)
                                                        cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGGT‑TTGGGGGCCC‑GGGGGGGTTTTT‑TTTTTCCCA‑AGGGGAAGGGCCACAGCCGGTTGGCCtttccccccctgccc                                                                                         >  2:635572/1‑138 (MQ=255)
                                                          cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGTGTA‑TTGGGGGCCA‑GGGGGGGTTTTT‑TTTTTTCCA‑CGGAGGCGGGCAAAAGCCGGTTGCCCTTCCCCCCCCGGGCCt                                                                                        >  2:1363468/1‑150 (MQ=255)
                                                          cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTg                                                                                       <  2:499679/151‑1 (MQ=255)
                                                             tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGTA‑TTTGGGGCCC‑GGGGGGGTTTTT‑TTTTTACCA‑GTGAGAAGGGGAACAGCTGATTTCCCTTCACCGCCCGGGCCTg                                                                                       >  1:238171/1‑148 (MQ=255)
                                                              ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                     <  1:1405622/149‑1 (MQ=255)
                                                               tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagag                                                                                   <  2:513954/150‑1 (MQ=255)
                                                                 cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGTA‑TTTGGGGCCA‑GGGGGGGTTTTT‑TTTTCACCA‑GGGAGACGGGGAACAGCTGATTTCCCTTCACCGCCCGGCCCCGagagag                                                                                 >  1:1371371/1‑150 (MQ=255)
                                                                 cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGTGCGCCA‑GGGGGGGTTTTC‑TTTTCACCA‑GGGGGACGGGGAACAGCCGATTGCCCTTCCCCCCCTCGCCCTGcgagag                                                                                 >  2:1517124/1‑150 (MQ=255)
                                                                         aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTTGGCGCCC‑GGGGGGGTTTTT‑TTTTCACCA‑GTGAGAAGGGCAACAACTTATTGCCCCTCACCGCCTGGCCCTGGGAGAGGTgcag                                                                           >  1:1108768/1‑148 (MQ=255)
                                                                          aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGT‑TTCGGA‑TTTTGCGCCC‑GGGGGGGTTTTT‑TTTTCTCCC‑GTGAGACGGGCAACAGCTTATTGTCCTTCCCCCCCCGGGCCCGAGAGGGGTgcag                                                                           >  1:658838/1‑147 (MQ=255)
                                                                            ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                     <  1:448389/151‑1 (MQ=255)
                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCC‑GGGGGGGTTTTT‑TTTTCACCA‑GTGAGAAGGGCAACAGCTGATTGCCCTTCCCCCCCCGGGCCCGGGAGAGGTGCAGCAAGccg                                                                    >  1:1051070/1‑149 (MQ=255)
                                                                                tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGT‑TGCGGC‑TTGGGCGCCA‑GGGGGGTTTTTC‑TTTT‑ACCACGGGAGACGGGGACCAGCCGGTTTCCCTTCCCCCCCCTGCCCCGGGAGAGGTGCCACAAGGAGTcc                                                                 >  2:68290/1‑151 (MQ=255)
                                                                                tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCT‑TGCGTA‑TTGGGCGGCC‑GGGGGGGTTTTC‑TTTTCAACA‑GAGAGGCGGGCGAAAGCCGCGTGCCCTTCACCGCCTGGCCCTtagag                                                                                   >  2:706966/1‑133 (MQ=255)
                                                                                    gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTTGGCGCCA‑GGGGGGGTTTTT‑TTTTTACCA‑GGGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGc                                                             >  1:1097235/1‑151 (MQ=255)
                                                                                      taGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGGA‑TTGGGCGGCA‑GAGAGGGTTTTC‑TTTTCACCC‑GAGATACGGGCACCCGATGGCTGGCCCCCCCCCCCCGGCCCTGGGAGAGCTgcacgcagcgggcccc                                                               >  2:931571/2‑134 (MQ=255)
                                                                                       aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTTGGCGCCA‑GGGGGGGTTTTT‑TTTTTACCA‑GTGGGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGCTgg                                                          >  1:1042021/1‑151 (MQ=255)
                                                                                         cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTCGGA‑TTTGGCGCCA‑GGGGGGGTTTTC‑TTTTCTCCA‑GTGAGACGGGCAACAGCTGAGTGCCCCTCACCCCCTGGCCCCGTGAGAGGTGCAGCCAGCGGGCCACGCTgg                                                          >  2:281923/1‑149 (MQ=255)
                                                                                         cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTTGGCGCCC‑GGGGGGTTTTTT‑TTTTTACCA‑GTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCTgggt                                                        >  1:461211/1‑149 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGCCGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                      <  2:672906/151‑1 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                       <  2:974620/150‑1 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                        <  1:1285369/149‑1 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                      <  1:421317/151‑1 (MQ=255)
                                                                                              ttAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTT‑TTTGGCGCCC‑GGGGGGGTTTTT‑TTTTCACCA‑GTGAGACGGGGAACAGCTGATTGTCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCCGGGTTGccc                                                   >  1:394429/1‑151 (MQ=255)
                                                                                                 aTGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TTGGCG‑ATGTGCGCCC‑GGGGGGGTTTTC‑TTTTTACCA‑GTGAGAAGGGGCACAGGTGGGTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGGAGCAAACGGGCCCCGCAGGTTTGcccc                                                  >  2:516641/1‑149 (MQ=255)
                                                                                                 aTGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                <  1:168812/151‑1 (MQ=255)
                                                                                                   gAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                              <  2:769426/151‑1 (MQ=255)
                                                                                                       cGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                           <  2:1466564/150‑1 (MQ=255)
                                                                                                       cGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                          <  2:935553/151‑1 (MQ=255)
                                                                                                           cAACGCGCGGGGAGAGGC‑GGTT‑TTTGGA‑TTGGGCGCCAGGGGGGGTTTTTT‑TTTTCCCCC‑G‑GGGACGGGCAACACCCGATTGCCCCTTACCCCCCCGCCCTGGGAGAGTTGCCGCAAGAGGGCCCCGCTGGGTTGTCCCCGCAGGCGgaa                                       >  2:1550817/1‑150 (MQ=255)
                                                                                                            taCGCGCGGGGAGAGGC‑GGTT‑TGACGA‑TTGGGCGCCC‑CGGGGGGTTTTT‑TTTTCACCA‑CAGTGACGGGCAACCCCCGCTTTCCCCCCAACCCCCGGCCCCGCGGCAGTTGCCGCAACCGGTCCACCCTGGGTTTCCCccgc                                               >  2:1041942/2‑141 (MQ=255)
                                                                                                             aCGCGCGGGGAGATGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                                    <  2:483437/151‑1 (MQ=255)
                                                                                                                   ggggAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACATCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                              <  1:837353/151‑1 (MQ=255)
                                                                                                                   ggggAGAGGC‑GGTT‑TGCGTA‑TTGGCCGCCG‑GGGTGGTTTTTC‑TTTTCCCCA‑GTGAGACGGGCAACAGTGGTTTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGAAAGCGGTCCACGCTGGTTTGCCCAAGCAGGCGAAAATCCTGtt                               <  2:1310951/150‑1 (MQ=255)
                                                                                                                    gggAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                             <  2:1071567/151‑1 (MQ=255)
                                                                                                                         agGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                        <  1:277726/151‑1 (MQ=255)
                                                                                                                           gC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                      <  2:1049199/151‑1 (MQ=255)
                                                                                                                               gTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑GTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                   <  1:100375/151‑1 (MQ=255)
                                                                                                                                tt‑tGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                   >  2:569090/1‑150 (MQ=255)
                                                                                                                                    gCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                <  1:782042/150‑1 (MQ=255)
                                                                                                                                      gTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg             >  2:43798/1‑151 (MQ=255)
                                                                                                                                          ttGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGCAAATCCTGTTTGATGGTGGTTACCGGCggg           >  2:150459/1‑150 (MQ=255)
                                                                                                                                            gggCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata        <  2:110967/151‑1 (MQ=255)
                                                                                                                                              gcgcCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat       >  1:273792/1‑150 (MQ=255)
                                                                                                                                               cgcCA‑GGGTGGTTTTTG‑TTTTCACCA‑GGGAGAAGGGCAACAGGTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGGGGTTAACGGCggggtgt       >  2:169787/1‑145 (MQ=255)
                                                                                                                                               cgcCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa     <  1:379222/151‑1 (MQ=255)
                                                                                                                                               cgcCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGCTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAACCTGTTTGATGGTGGTTAACGGCGGGgta        >  1:737950/1‑148 (MQ=255)
                                                                                                                                                gcCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa     >  2:1348288/1‑150 (MQ=255)
                                                                                                                                                 ccA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa     >  1:1474241/1‑149 (MQ=255)
                                                                                                                                                  cA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt  <  2:794075/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTT‑TGCGTA‑TTGGGCGCCA‑GGGTGGTTTTTC‑TTTTCACCA‑GTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAT  >  NC_000913/366378‑366671

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