Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5260TG45.8% 59.6 / 21.6 49T329P (ACC→CCC) lacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (6/20);  new base G (22/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAT  >  NC_000913/366378‑366671
                                                                                                                                                         |                                                                                                                                                 
aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTT‑G‑TA‑TTTGGGGCC‑Cgggggg                                                                                                                                                 >  2:762728/1‑151 (MQ=255)
     tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGATTGG‑G‑TG‑TTGGGGGGC‑CGGGGGGgtttt                                                                                                                                            >  1:718417/1‑151 (MQ=255)
        ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTttttctt                                                                                                                                         <  1:507969/151‑1 (MQ=255)
         gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑GT‑TTTGGGGCC‑CGGGGGGGtttttttt                                                                                                                                        >  2:435050/1‑151 (MQ=255)
           gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTttttcttttc                                                                                                                                      <  1:117786/151‑1 (MQ=255)
           gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTttttctttt                                                                                                                                       <  1:588929/150‑1 (MQ=255)
                 tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGG‑G‑TA‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTTAcc                                                                                                                                   >  2:415519/1‑148 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCG‑G‑TT‑TTGGGGCCC‑GGGGGGGTTTTTTTTTTC‑CCAGt                                                                                                                                >  1:85260/1‑148 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑TATTTGGGGGCC‑GGGGGGGttttttttttacc                                                                                                                                    >  2:758832/1‑142 (MQ=255)
                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑TA‑TTTGGCGGC‑CGGGGGGGTTTTTTTTTTTCCAGTg                                                                                                                               >  2:613869/1‑149 (MQ=255)
                     tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑GA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTAACAGt                                                                                                                                >  2:137313/1‑147 (MQ=255)
                      gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGGGCC‑AGGGGGGGTTTTTTTTTTACCCGTgaga                                                                                                                            >  1:215251/1‑150 (MQ=255)
                      gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                            <  1:530384/150‑1 (MQ=255)
                        gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGG‑GT‑TTGGGGGCC‑GGGGGGGTTTTTTTTTTC‑CCCGTGAGACgg                                                                                                                         >  1:399817/1‑151 (MQ=255)
                         cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTG‑G‑TA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCCGTGAGAAggg                                                                                                                        >  2:15259/1‑151 (MQ=255)
                          tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGG‑G‑TA‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTTTCCAGTGAGAAggg                                                                                                                        >  2:568844/1‑150 (MQ=255)
                              tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTTGGGGCC‑AGGGGGGTTTTTCTTTTTACCAGGGAGACGGGCAAc                                                                                                                    >  2:596907/1‑150 (MQ=255)
                               cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑GTTT‑TTGGGGGCC‑GGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGGAACa                                                                                                                   >  2:820924/1‑150 (MQ=255)
                                 caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                                 <  2:913703/150‑1 (MQ=255)
                                 caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                                <  1:179907/151‑1 (MQ=255)
                                    tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCACCGCGCGGGGAGAGGC‑GGCTTTC‑G‑TC‑ATTTGCGCC‑AGGGCGGTTTTTCCTTTTACCCAGGGGAAGGGCAACacatgg                                                                                                              >  1:471826/1‑145 (MQ=255)
                                     aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑GCTT‑TTGGGGGCC‑‑GGGGGGGTTTTCTTTTTCCCCGCGAGACGGGGAAACACTGAtt                                                                                                            >  2:535547/1‑151 (MQ=255)
                                          cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                       <  2:619244/151‑1 (MQ=255)
                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑TA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTTCCCGTGAAACGGGGAAAAACCGATTTCCCTTCACcccc                                                                                               >  2:678440/1‑149 (MQ=255)
                                                         ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                        <  1:35740/151‑1 (MQ=255)
                                                                      gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGT‑G‑TG‑ATTGGCGCC‑AGGGGGGGTTTTTTTTTTACCAGCGGGAAGGGGAACAGCTGGTTGCCCCTTACCGCCTGGCCCCGGGAGAGTTgc                                                                             >  2:155189/1‑149 (MQ=255)
                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGT‑G‑TA‑TTTGGCGCC‑AGGGGGGGTTTTTTTTTTACCAATGAGACGGGGAACAACTGATTGACCCTCACAGCCCGGCCCTGAGAGAGGTgcag                                                                           >  2:165343/1‑149 (MQ=255)
                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑CG‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGCAACAGCTGATTGCCCCTCCCCCCCCGGCCCTGAGAGAGGTg                                                                              >  1:683285/1‑146 (MQ=255)
                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTGTCG‑G‑GA‑TTTGGCGCCGCGGGGGGTTTTTTTTTTC‑CCCGGGAGGAGGGGAACAGCTCATTGCCCTTCACCCCCCGGCCCTGagagaggtagaggcagcgg                                                                    >  1:346943/1‑138 (MQ=255)
                                                                                tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                                  <  2:580114/150‑1 (MQ=255)
                                                                                     ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCTTTT‑G‑TA‑TTGGGGGCC‑CCGGTGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCTGGCCCTGGGGGAGGTGCAGCAAGCGGccccc                                                               >  1:874272/1‑146 (MQ=255)
                                                                                         cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTC‑G‑TA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCAGTGGGACGGGCAACAGCTGATTTCCCTTCACCCCCTGGCCCCGAGAGAGGTGGAGCAAGCGGGCCCCGCTgg                                                          >  2:776602/1‑149 (MQ=255)
                                                                                           gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTG‑G‑TA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTCCCCGTGAGACGGGGAAAAGCCGATTGCCCTTTACCCCCCGGCCCCGAGAGAGTTGGAGCAAGCGGGCCACGCTGGgtt                                                       >  2:76643/1‑150 (MQ=255)
                                                                                             aTTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGTCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                        <  2:812395/147‑1 (MQ=255)
                                                                                                       cGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                          <  2:943061/151‑1 (MQ=255)
                                                                                                          ccAACGCGCGGGGAGAGGC‑GGTTTGG‑G‑GA‑TTTGGGGCC‑AGGGGGGGTTTTCTTTTTACCAGTGAGAAGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCCCGCTGGGTTGCCCCCGCCGGCGaa                                        >  2:56664/1‑150 (MQ=255)
                                                                                                               gcgcgGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                  <  1:269695/151‑1 (MQ=255)
                                                                                                                      gagaGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                           <  2:662292/151‑1 (MQ=255)
                                                                                                                      gagaGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                           <  1:609264/151‑1 (MQ=255)
                                                                                                                           gC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                      <  2:715496/151‑1 (MQ=255)
                                                                                                                              ggTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                    >  2:910864/1‑151 (MQ=255)
                                                                                                                                  tGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                >  2:108219/1‑151 (MQ=255)
                                                                                                                                           ttGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCgggg          >  1:800396/1‑150 (MQ=255)
                                                                                                                                           ttGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg           <  1:921530/150‑1 (MQ=255)
                                                                                                                                             gggCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata        <  1:137303/151‑1 (MQ=255)
                                                                                                                                               gcgcC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat         >  2:618492/1‑148 (MQ=255)
                                                                                                                                                  cc‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa   <  1:894366/151‑1 (MQ=255)
                                                                                                                                                   c‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt  <  2:900351/151‑1 (MQ=255)
                                                                                                                                                      gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa   <  1:79663/148‑1 (MQ=255)
                                                                                                                                                      gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATAACa   >  2:487879/1‑148 (MQ=255)
                                                                                                                                                         |                                                                                                                                                 
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGC‑G‑TA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAT  >  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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑