Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5390CT47.2% 50.4 / 26.1 53V324V (GTG→GTAlacIlactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base T (25/0);  total (38/15)
Fisher's exact test for biased strand distribution p-value = 6.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAT  >  NC_000913/366393‑366687
                                                                                                                                                    |                                                                                                                                                       
aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                                     <  1:896561/151‑1 (MQ=255)
 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                                      >  1:296935/1‑149 (MQ=255)
             aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                                       <  2:341218/151‑1 (MQ=255)
               tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGGATTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCaa                                                                                                                                       >  2:859052/1‑149 (MQ=255)
                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGGGGCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACa                                                                                                                                     >  2:104448/1‑146 (MQ=255)
                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAAc                                                                                                                                      >  2:759837/1‑145 (MQ=255)
                        ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGAtt                                                                                                                              >  2:353014/1‑149 (MQ=255)
                            gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTGGGGGGCCGGGGGGGTTTTTTTTTTTCCAGGGAGAAGGGGGACAGCTGATTGtcc                                                                                                                        >  1:748805/1‑149 (MQ=255)
                             ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTGTTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAACTGATTGcccct                                                                                                                      >  1:567202/1‑148 (MQ=255)
                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                                                <  2:885320/149‑1 (MQ=255)
                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGGGGTTTTGGCGCCCGGGGGGGTTTTCTTTTTCCCCGTGGGACGGGCAACAGCTGATTGCCCCTCACcccc                                                                                                               >  1:195724/1‑149 (MQ=255)
                                       tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTGTTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGGAGGGGAACAGCTGGTTGCCCTTCACcccc                                                                                                               >  1:764185/1‑147 (MQ=255)
                                        gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGTATTGGGGGCCCGGGGGGGTTTTTTTTTCACCCAGGAGACGGGGCACAGCCGGTTGGCCCTTCccccccgg                                                                                                             >  1:486255/1‑147 (MQ=255)
                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTGTTGGGGGGCAGGGGGGGTTTTTTTTTTCCCCGGGAGACGGGCAAAACCCGATTGCCCCTCACCCCCCGGccc                                                                                                         >  2:781833/1‑150 (MQ=255)
                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCCCCCTGccc                                                                                                         >  2:224325/1‑150 (MQ=255)
                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCCCCTGGccc                                                                                                         >  2:221220/1‑150 (MQ=255)
                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                                        <  2:106325/151‑1 (MQ=255)
                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGATTGGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGGACAGCTTGTTGCCCCTCCCCCCCCGGccc                                                                                                         >  1:915757/1‑148 (MQ=255)
                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTGGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGGCGGGCAACAACTGATTGCCCTTCACCGCCCGGCCCCGagagaggt                                                                                               >  2:878026/1‑146 (MQ=255)
                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGCCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGCACAGCTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTg                                                                                              >  2:713345/1‑149 (MQ=255)
                                                       gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAACAACTGATTGCCCTTCACCCCCCGGCCCCGagagaggt                                                                                               >  2:244512/1‑145 (MQ=255)
                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGATTTCCCTTCACCCCCTGGCCCTGAGAGAGGTgcagca                                                                                         >  2:1116175/1‑150 (MQ=255)
                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAAAAGCGGATTGCCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCaa                                                                                        >  2:123952/1‑150 (MQ=255)
                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCaa                                                                                        >  2:258877/1‑149 (MQ=255)
                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTCTATTTTGGCGCCGGGGGGGTTTTTTTTTTACCCGTGAGACGTGCACCAGCACCTTGACCCTCCACCCCTGTTCCTCagagaggtttatcaatccg                                                                                    >  1:209794/1‑137 (MQ=255)
                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGAAGGGCAACAGCTGATTGGCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCggg                                                                                   >  1:1025718/1‑150 (MQ=255)
                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                                                  <  1:218461/151‑1 (MQ=255)
                                                                  cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGCTTTGGGGGCCGGGGGGGTTTTTTTTTTACCCCTGGGACGGGGCACAGCCTAATGTCCTTCACCCCCCCGCCCCGTGAGAGtgggcgccagccggcc                                                                                 >  1:280624/1‑135 (MQ=255)
                                                                  cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                                <  1:136403/151‑1 (MQ=255)
                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCCTGAGACGGGGAACAGCTGATTTCCCCTCCCCCCCCGGCCCTGAGAGAGGTGCAGCCAACGGTCCAc                                                                               >  1:563343/1‑147 (MQ=255)
                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                                     <  1:598537/151‑1 (MQ=255)
                                                                              aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGGTTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGGGACAGGCAAAAACCGATTGCCCTTCCCCGCCAGGCCCTCGGAGAGGTGCAGCAAGAGGGCCACGCCGGgtt                                                                       >  2:760125/1‑148 (MQ=255)
                                                                              aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAAAAGCTGATTTGCCTTCACCGCCTGGGCCTGAGAGAGTTGCAGCAAACGGGCCACGCCGGGTTTcc                                                                    >  2:740777/1‑151 (MQ=255)
                                                                                tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTCCCCGTGAGACGGGCAACAACCGATTGCCCCTTACCGCCCGGCCCTGAGAGAGGTGGAGCAAACGGGCCACGCCGGGTTTcccc                                                                  >  2:1015672/1‑151 (MQ=255)
                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                                <  2:395138/151‑1 (MQ=255)
                                                                                    gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                              <  1:416480/151‑1 (MQ=255)
                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                             <  1:1152805/151‑1 (MQ=255)
                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                                         <  2:995191/151‑1 (MQ=255)
                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                                         <  1:1119601/151‑1 (MQ=255)
                                                                                           ccAACGCGCGGGGAGAGGCGGTTTGGCTCATAGGGGCCCGGGGGGTTTTTCTTTTCATCAGTGAGACCGGGCAAAGCGGACTGGCCTTTTCCGCCTTGCCCCGCGGGAGGTGCAGCAAGCCGGTCACCCTGGGGTGCCCCAGCAGGcaca                                                        >  2:611104/1‑147 (MQ=255)
                                                                                                              ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                                    >  2:560599/1‑151 (MQ=255)
                                                                                                                 tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                                   >  1:1105291/1‑150 (MQ=255)
                                                                                                                        aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                           >  2:523826/1‑151 (MQ=255)
                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                           >  2:577047/1‑145 (MQ=255)
                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata                      >  1:155657/1‑150 (MQ=255)
                                                                                                                                  cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa                   >  1:803836/1‑148 (MQ=255)
                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa                   >  1:524799/1‑147 (MQ=255)
                                                                                                                                     ggtggtTTTTCTTTTCACCATTGAGACGGGCAACAGCTGATTGCCCTGCACCGCTTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg                 <  2:751533/149‑1 (MQ=255)
                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg               <  1:581916/151‑1 (MQ=255)
                                                                                                                                      gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATAACATGa                >  1:333607/1‑149 (MQ=255)
                                                                                                                                      gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg                 >  1:1117874/1‑148 (MQ=255)
                                                                                                                                        ggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt             <  1:759800/150‑1 (MQ=255)
                                                                                                                                         gtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGt           >  2:307292/1‑151 (MQ=255)
                                                                                                                                                    cACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGCTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATCACATTAGCTGTCTTCGGTAt  >  1:946323/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
AATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAT  >  NC_000913/366393‑366687

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: