Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,756,565	0	C	T	21.4%	46.8 / 11.4	28	intergenic (‑54/+111)	aldB/yiaY	aldehyde dehydrogenase B/L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base T (4/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

GTGTGCGCCCATTTATCTTTCACTTTGTGCGCAGCATCCAGCGCCAGATCGATGTCTCGTTTGCCCGAAGACGCCACTTCGCACAGCAGCTGCCCGGTCACCGGCGTCAGATTCTGGTAATACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATGGCGTTAATGC  >  NC_000913/3756270‑3756860                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        tgctgcgcacaaagtgaaagataaatgggcgcAGCATCCAGCGCCAGATCGATGTCTCGTTTGCCCGAAGACGCCACTTCGCACAGCAGCTGCCCGGTCACCGGCGTCAGATTCTGGTAATACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGCAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCg                                                                                                                                                                                                                                                                                                     >  2:209641/29‑300 (MQ=255)      cgcCCATTTATCTTTCACTTTGTGCGCAGCATCCAGCGCCAGATCGATGTCTCGTTTGCCCGAAGACGCCACTTCGCACAGCAGCTGCCCGGTCACCGGCGTCAGATTCTGGTAATACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGATTTCGCCAGt                                                                                                                                                                                                                                                                                                <  2:192360/300‑1 (MQ=255)                                          cgcCAGATCGATGTCTCGTTTGCCCGAAGACGCCACTTCGCACAGCAGCTGCCCGGTCACCGGCGTCAGATTCTGGTAATACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACgg                                                                                                                                                                                                                                                           <  2:80459/301‑1 (MQ=255)                                                                               tCGCACAGCAGCTGCCCGGTTACCGGCGTCAGATTCTGGTAATACTCGCCGTAGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTTAGGGGGAAACCATACTAGCCGGGCTTAATCTGTGCTGAAGGGGGAGTATTGGTCATGATCGTCTCCTTGCGGTGTTAGGTATAACAAGGGTAGACGTGACTAGCGAAATCTTTGCCAGTCGGTAACAGCTTTACACAAGCATACCGGAATTTCAGGAAAAGGGTTAATGAAATAGTTGAATAAgagtcacaca                                                                                                                                                                                                                      >  1:257194/1‑292 (MQ=255)                                                                                                  tCAACGGCGTGAGATTCTGGTAGTACTCGCCGTCGGCAGGGTCTACCCATGCTGCGACAAGAAAGTTGTCATCGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCTGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGa                                                                                                                                                                                                   <  2:257513/301‑1 (MQ=255)                                                                                                                         tACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGcgcg                                                                                                                                                                            >  2:108873/1‑301 (MQ=255)                                                                                                                                    ggCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCca                                                                                                                                                                 >  2:100822/1‑301 (MQ=255)                                                                                                                                         gggCAAACCAATCGGCGACAAAAAGGTTGTCATAGCGGGCTTTAAACTGGAAGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGCTCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAAt                                                                                                                                                             <  1:209392/300‑1 (MQ=255)                                                                                                                                                                     gTCATAGCGGGCTTTTAACTTGAGTGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGACTTTACGGGCAAGGGTTACAGAGCTAGTTGCAAAAGCGGCTAACACTTTGCCAGCCGTTGTGGACATGCTGATTACATCGCTGCGCGCTAAATCTCCACACTTTCTTCGTGATTTGCCTGGCACGGGTTAg                                                                                                                                >  1:99466/1‑301 (MQ=255)                                                                                                                                                                               gCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCcaca                                                                                                                      >  1:138254/1‑301 (MQ=255)                                                                                                                                                                                    tAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCAt                                                                                                                 >  2:226547/1‑301 (MQ=255)                                                                                                                                                                                                                    tAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGaaa                                                                                 >  2:56785/1‑301 (MQ=255)                                                                                                                                                                                                                    ccgatcttgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGaaa                                                                                 <  2:15664/294‑1 (MQ=255)                                                                                                                                                                                                                     aaTGTGGGCGGAAGGGGTATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAAt                                                                                <  1:73592/301‑1 (MQ=255)                                                                                                                                                                                                                           tgCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGcaaagatcgg                                                                          >  1:15644/1‑294 (MQ=255)                                                                                                                                                                                                                            gCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAAtcttctt                                                                          >  2:155142/1‑300 (MQ=255)                                                                                                                                                                                                                              tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATcttcttctt                                                                       >  1:180372/1‑301 (MQ=255)                                                                                                                                                                                                                              tGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATcttcttctt                                                                       <  2:47748/301‑1 (MQ=255)                                                                                                                                                                                                                                                                    ccgatctataACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACt                                 <  1:30094/295‑1 (MQ=255)                                                                                                                                                                                                                                                                     tgAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACt                                 <  1:56731/300‑1 (MQ=255)                                                                                                                                                                                                                                                                           ataACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACtagatcgg                          >  2:30129/1‑294 (MQ=255)                                                                                                                                                                                                                                                                               tcttccgatctccgACTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAg                      <  1:56155/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                         cGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATg            <  2:138391/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                          ccgACTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAgagatcggaaga           >  2:56209/3‑290 (MQ=255)                                                                                                                                                                                                                                                                                           cgACTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATgg           >  2:106994/2‑300 (MQ=255)                                                                                                                                                                                                                                                                                            gACTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATGGCg         >  2:191843/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                             aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCCTCGCCAGTTCACGGATGGCGt        >  1:41007/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                    aaaTCTTCGCCAGTCGGTACCAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATGGCGTTAATGc  <  2:160313/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        GTGTGCGCCCATTTATCTTTCACTTTGTGCGCAGCATCCAGCGCCAGATCGATGTCTCGTTTGCCCGAAGACGCCACTTCGCACAGCAGCTGCCCGGTCACCGGCGTCAGATTCTGGTAATACTCGCCGTCGGCAGGGGCTACCCATTCGCCGCCAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACCGCGAAATCTTCTTCTTTCACGTTCAGGTCGCGTAGGCCTGCCGGGATATCCACTTTCTTCGCCAGTTCACGGATGGCGTTAATGC  >  NC_000913/3756270‑3756860

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑