Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	466,776	C→T	G179G (GGG→GGA)	ybaE ←	uncharacterized protein YbaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	466,776	0	C	T	90.9%	27.6 / NA	11	G179G (GGG→GGA)	ybaE	uncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/3);  minor base G (0/1);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCGCCCATAAACGGCTGTAACAGAGTGCGCAGCTCACCTGGGGCCAGTTG  >  NC_000913/466595‑467021                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCgg                                                                                                                                                                                  >  1:10600/1‑251 (MQ=255)                                                    cGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATAGATGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCACTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGc                                                                                                                               >  1:64496/1‑251 (MQ=255)                                                                    aTCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTccaag                                                                                                                                                                                                                                  >  1:441716/1‑133 (MQ=255)                                                                                                        cTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACg                                                                           >  1:204792/1‑251 (MQ=255)                                                                                                                aGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGt                                                                                                          >  1:17230/1‑212 (MQ=255)                                                                                                                                  cATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCg                                                 <  1:343267/251‑1 (MQ=255)                                                                                                                                                                   gATAAAAGGCCCAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGATGAGTATTATTCTCGAAGCGGGTCAGGCCGGAACATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCTTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCGCCCATAAACGGCTGTACCAGAGTGCGCAGCTc                 >  1:322564/1‑250 (MQ=255)                                                                                                                                                                       aaaGTCCCAGCGTAAGCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCTCCAATCGGGCGCTGAGTATTATTAGCGAAGCGGGTCAGGCCGTAAAATATCTGCCCGGCGTGATTCTTCTCTGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAAt                                                                                   <  1:256656/180‑1 (MQ=255)                                                                                                                                                                            cccAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAACATATCTGACCGGAGAGATGCTGCTCGGCACGCCCGGGCCAAAAGCCTGGTTGTAGTGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCGCCCATAAACGGCTGTAACAGAGTGCGCAGCTCACCTGGGGcc       >  1:441708/1‑251 (MQ=255)                                                                                                                                                                             ccAGCGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGGTTTGCCATTGTCCGCCCATAAACGGCTGTAACAGAGTGCGCAGCTCACCTggg          <  1:429022/247‑1 (MQ=255)                                                                                                                                                                                 cGTAATCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCGCCCATAAACGGCTGTAACAGAGTGCGCAGCTCACCTGGGGCCAGTTg  <  1:30999/251‑1 (MQ=255)                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCCATTGTCCGCCCATAAACGGCTGTAACAGAGTGCGCAGCTCACCTGGGGCCAGTTG  >  NC_000913/466595‑467021

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