Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	46.7%	52.2 / 26.8	45	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (3/21);  new base G (21/0);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 6.70e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAA  >  NC_000913/366377‑366668                                                                                                                                                       |                                                                                                                                               aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggt                                                                                                                                                <  2:928287/150‑1 (MQ=255)  aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggtgg                                                                                                                                              <  2:966603/151‑1 (MQ=255)       tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTC‑TTGGGCGCCCGGGGGGgtttt                                                                                                                                         >  1:659808/1‑151 (MQ=255)       tGGGGTGCCTAAGGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGGGCGGGGAGAGGCGGTTTGCGAA‑TTGGGAGTCGGGGGGGgtttt                                                                                                                                         >  1:249010/1‑151 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGA‑TTGGGGGCCAGGGGGGGttttttttt                                                                                                                                    >  1:242479/1‑150 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                           <  2:519934/151‑1 (MQ=255)                        gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGGGGCCGGGGGGGttttttttttacc                                                                                                                                 >  2:101745/1‑140 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGTCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCa                                                                                                                   <  1:801846/151‑1 (MQ=255)                               cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTG‑TTTGGCGCCCGGGGGGGTTTTTTTTTCTCCAGGGAGACggggga                                                                                                                  >  1:713578/1‑147 (MQ=255)                                 cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                <  2:714449/150‑1 (MQ=255)                                    aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                                                                                            <  2:704313/151‑1 (MQ=255)                                     ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGG‑TTTGGCGGCCGGGGGGTTTTTCTTTTTCCCCGTGAGACGGGGaa                                                                                                                  >  1:812751/1‑144 (MQ=255)                                           gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                     <  2:574152/151‑1 (MQ=255)                                            cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGT‑TTTGGGGCCCGGGGGGGTTTTCTTTTTCCCCGTGAGACGGGCAAAAGCTGATTGccc                                                                                                     >  2:588887/1‑150 (MQ=255)                                                gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                                                                                                 <  1:469111/150‑1 (MQ=255)                                                gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                                                                                                 <  1:162792/150‑1 (MQ=255)                                                gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                                                                                                 <  1:766405/150‑1 (MQ=255)                                                  gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTTGTG‑TTTGGCGGCAGGGGGGGTTTTTTTTTTACCAGAGAGATGGGCAACAGCTCATTGTCCTTCAcc                                                                                               >  1:884801/1‑150 (MQ=255)                                                    tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGCACAGCCGGTTGCCCCTCACCGcc                                                                                            >  1:774912/1‑151 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                             <  1:825394/149‑1 (MQ=255)                                                      aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                          <  1:86684/151‑1 (MQ=255)                                                        tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGGGAGACGGGCAACAGCTGATTGGCCTTCACCCCCCgg                                                                                         >  1:566295/1‑150 (MQ=255)                                                         gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                                                                                        <  1:69309/150‑1 (MQ=255)                                                            cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGccc                                                                                      >  2:676898/1‑149 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                                 <  1:1019612/150‑1 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCGTGGGACGGGGCACAACTGAATGCCCTTCACCCCCCGGccc                                                                                      >  1:719218/1‑145 (MQ=255)                                                                      tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGTTTTTTTTTTTTCCAGTGGGACGGGGAACAGCCGGTTGCCCCTCACCCCCCTGCCCCGggagaggt                                                                            >  1:473203/1‑147 (MQ=255)                                                                         ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGGGAGAAGGGCAACAGCTGATTGCCCTTTCCCGCCTGGgcccggg                                                                                  >  2:699276/1‑136 (MQ=255)                                                                           aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTCTCCAGTGAGGCGGGGAACAGATGATTTGCCCTCCCCGCCGGGGCCTTggagaggttca                                                                         >  1:984403/1‑142 (MQ=255)                                                                           aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                     <  2:249021/151‑1 (MQ=255)                                                                             aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGGGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCCGGGAGAGGTGCAGCaa                                                                     >  2:199507/1‑149 (MQ=255)                                                                                 gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                <  1:743623/150‑1 (MQ=255)                                                                                  tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAACTTATTGCCCCTCACCGCCTGGGCCTGAGAGAGTTGCAGGAAGCCGTcc                                                              >  2:277802/1‑151 (MQ=255)                                                                                      gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGTTTTTTTTTTTACCAGGGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCCGAGAGAGTTGCAGCAAGCGGTCCACGc                                                          >  2:227925/1‑151 (MQ=255)                                                                                                      aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                          <  2:450579/151‑1 (MQ=255)                                                                                                        tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                         <  2:513687/150‑1 (MQ=255)                                                                                                         cGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCCGATTGCCCCTCACCCCCCGGCCCCGAGAGAGGTGCAGCAAACGGGCCACGCCGGGTTGGCCCCGCCgg                                         >  1:655537/1‑149 (MQ=255)                                                                                                         cGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                               >  2:186911/1‑143 (MQ=255)                                                                                                             cAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAACAGGTGATTGCCCTTCACCCCCTTGCCCTGGGGGAGGTGCCGCAAGCGGTCCACGCTGGGTTGGCCCAGCACGGGaaa                                    >  2:100044/1‑150 (MQ=255)                                                                                                                          gagGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                      <  1:676847/151‑1 (MQ=255)                                                                                                                               ggTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                 <  1:511275/151‑1 (MQ=255)                                                                                                                                gTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                <  2:721979/151‑1 (MQ=255)                                                                                                                                        a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg         <  1:362362/150‑1 (MQ=255)                                                                                                                                        a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg        >  1:521479/1‑151 (MQ=255)                                                                                                                                           tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa       <  1:415657/150‑1 (MQ=255)                                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCGGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata     >  1:271858/1‑148 (MQ=255)                                                                                                                                                gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa  >  1:452245/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                               AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAA  >  NC_000913/366377‑366668

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