Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,883	0	C	T	48.4%	17.5 / 23.4	31	intergenic (+52/‑51)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/15);  new base T (15/0);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 5.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGAC  >  NC_000913/2386750‑2387029                                                                                                                                      |                                                                                                                                                      tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa                                                                                                                                     <  2:256062/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa                                                                                                                                     <  1:130958/151‑1 (MQ=255) tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGa                                                                                                                                     <  1:360121/151‑1 (MQ=255)     aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAATGCACATTTTTTTTTTCCTCA‑TCtttttttaaaaggg                                                                                                                                 >  1:205061/1‑143 (MQ=255)     aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCctcc‑tc                                                                                                                                               >  2:80650/1‑137 (MQ=255)     aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑Tc                                                                                                                                               <  1:80648/137‑1 (MQ=255)         cagacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGGAAAAAGGAAATGCTCAATTTATTTTTTCTCATTCTGtttttaaagcggaa                                                                                                                               >  1:227579/1‑140 (MQ=255)                                        ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAATGAACATTTTATTTTTTCTCA‑TCTGTTTTTAAAAGC‑GAAACATG‑AAAATAAAAGGAGGGAAAAAAAATAGGa                                                                                            >  1:447668/1‑151 (MQ=255)                                        ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAGCGAAAATCAAATTTTTTTTTTTCCCA‑TTTGTTTTTAAAAGA‑GAAAAATG‑AAATTTAACGGAGGAAAAAAAAATgggg                                                                                            >  1:484926/1‑150 (MQ=255)                                        ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGGAAAATTTTTTTTTCCTCC‑TCTGTTTTTAGATGT‑GAAAAATG‑AAAATTAACCGGGGAAAAATAAATgggg                                                                                            >  1:157446/1‑150 (MQ=255)                                        ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGGACATTTTTTTTTTCCTCA‑TCTGGTTTAAGATGC‑GGAACATG‑AAAATCAAAGGGAGAAAAAAAAATgggg                                                                                            >  2:272017/1‑150 (MQ=255)                                        ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGATTCATAAAAAAAGACATACCAATTTTTTTTTTTCTCA‑TCTGTTTTTAGATTG‑GAAACATG‑AAAATCTACGGAGGGGAAATAAATaggg                                                                                            >  1:127686/1‑150 (MQ=255)                                                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGAACATTTTATTTTTCCTCA‑TCTGTTTTTAGATGC‑GAAAAATG‑AAATTTAAAAGAGGAAAAAAAAAATGGAATGATTGGCTa                                                                                 >  1:199273/1‑151 (MQ=255)                                                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATAAACAATTTATTTTTCCTCC‑TCTGTTTTTAAAAGA‑GGAAAATG‑AAAATTAAATGAGGGAAAATAAATGGGGATTATTgggtt                                                                                 >  1:324893/1‑148 (MQ=255)                                                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCA‑TCTGTTTTTAGATGG‑GAAAAATG‑AAAATTAAAGGGGGGAAAAAAAATGGGGATGATTGGGTa                                                                                 >  2:447352/1‑151 (MQ=255)                                                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGGAATGCACAATTTATTTTTCCTCA‑TCTGTTTTTAAAAGCGGAAAAATA‑AAATTAACGGGGGGAAAAAAAAATGGG‑ATGATTggttt                                                                                 >  2:203900/1‑148 (MQ=255)                                                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAAAAAAAAAGAAATGCACATTTTATTTTTTCTCC‑TATGTTTTAAGAAGA‑GAAAAATGAAAATTTAAAGGGGGAAAAA‑AAATTGGAAAGAATGGGTa                                                                                 >  1:550217/1‑151 (MQ=255)                                                    gtgtataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCA‑TCTTTTTTTAGATGC‑GAAAAATG‑AAAATCAACGGAGGAAAAATAAATGGGAATGATTGGCTAt                                                                                >  2:366477/4‑151 (MQ=255)                                                                 ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCCCATTTTATTTTCCCTCA‑TCTGTGTTAAGATGC‑GAAACATT‑GAATTCAACGGAGGAAAAATATATGGTAATGATTGGCTATTTTGCAGAAATTg                                                                   <  1:214278/151‑1 (MQ=255)                                                                         caGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTTTTTTTTCTCC‑TCTGTTTTTAGATGC‑GAAAAATG‑AAAATCTAAGGAGGGAAAAAAAATTGGGATGATTTGGTATTTTTCAGAAAATTACTCAGaa                                                           >  1:249617/1‑151 (MQ=255)                                                                                                                  cGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                  <  1:115260/151‑1 (MQ=255)                                                                                                                           caTTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat         <  2:494894/151‑1 (MQ=255)                                                                                                                           caTTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat         <  2:110826/151‑1 (MQ=255)                                                                                                                             ttttattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc       <  1:286965/151‑1 (MQ=255)                                                                                                                               ttattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACAt     >  2:508322/1‑151 (MQ=255)                                                                                                                                tattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg    <  2:504768/151‑1 (MQ=255)                                                                                                                                tattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg    <  1:114169/151‑1 (MQ=255)                                                                                                                                tattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg    <  2:497685/151‑1 (MQ=255)                                                                                                                                tattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg    <  2:438881/151‑1 (MQ=255)                                                                                                                                tattttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg    <  2:227583/151‑1 (MQ=255)                                                                                                                                  ttttCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGac  <  1:194049/151‑1 (MQ=255)                                                                                                                                      |                                                                                                                                                      TGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGAC  >  NC_000913/2386750‑2387029

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