Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	377,432	0	C	G	50.0%	1.4 / 28.3	16	intergenic (+121/+103)	yaiL/frmB	DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base G (8/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATT  >  NC_000913/377188‑377675                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    tGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACTCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCt                                                                                                                                                                                                                                               >  1:246382/1‑251 (MQ=255)                                   ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACg                                                                                                                                                                                                             >  1:23646/1‑251 (MQ=255)                                                                   gCAGGATGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTTCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCg                                                                                                                                                                             <  1:137955/251‑1 (MQ=255)                                                                                                                                                                                              aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg                                                                                                                                                                                                    >  1:167035/1‑105 (MQ=255)                                                                                                                                                                                              aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg                                                                                                                                                                                                    <  2:167035/105‑1 (MQ=255)                                                                                                                                                                                                       gggtacaCCGGATAAGGCGTCCCCGCCCGC‑TCCGGCAGGCGGGCCC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa                                         <  2:147800/245‑1 (MQ=255)                                                                                                                                                                                                                     atttgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTATTGATAtt                             >  1:534026/24‑251 (MQ=18)                                                                                                                                                                                                                        tgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCAt                          >  1:513249/21‑251 (MQ=18)                                                                                                                                                                                                                              cccccccacccGGCCGTCGTGCGC‑CATGCCTGATGCGACGTTTTTCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCATTCGTGCTATTATCAACGCATATTCAGTTTATTGTCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTTATAACGGATTAACCTCTTGATATTCATCTCCTgg                   <  2:246382/242‑1 (MQ=255)                                                                                                                                                                                                                                taaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                  >  1:144981/13‑251 (MQ=21)                                                                                                                                                                                                                                  aacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGACGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGACAATATTCTCGCCAATAAAACTAGAGACAAAATAATAGCTGAGATCATAACCCTCGTAATAACGGCTTAACGACTTGGTATTCCTCTCCttggcg                >  1:81996/11‑246 (MQ=17)                                                                                                                                                                                                                                    catagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                  >  1:36749/9‑247 (MQ=25)                                                                                                                                                                                                                                          taaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCCCGGGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc        >  1:352085/3‑251 (MQ=25)                                                                                                                                                                                                                                            aGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGGAGTGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGAAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCCTAACCCTCGTGATAAAGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTcc       >  2:187822/1‑249 (MQ=25)                                                                                                                                                                                                                                             gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTTTGATCCTAACCCTCGTGATAACGGATTAACGTCTTGATATTAATCTCCTGGCAGATCTTTTCa       >  2:359808/1‑249 (MQ=25)                                                                                                                                                                                                                                               cGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGAtt  <  2:23646/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    TGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATT  >  NC_000913/377188‑377675

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