Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	863,000	G→T	S171Y (TCT→TAT)	ybiY ←	putative pyruvate formate‑lyase activating enzyme YbiY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	863,000	0	G	T	94.4%	52.4 / ‑5.6	18	S171Y (TCT→TAT)	ybiY	putative pyruvate formate‑lyase activating enzyme YbiY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (11/6);  minor base C (0/1);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATGATCTCCTCAACGCTTTTCACTTCACCACACACGGTTA  >  NC_000913/862756‑863222                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               gggaagaaaaagactgtggCCAACGTGCAGCTCGTCGTCGGCAAAATCGGTAATGGCTTTTTAAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATCTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCATCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGGCGGCAACGTGGTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGACAAgagg                                                                                                                                                                                                                          <  2:56857/233‑4 (MQ=255)        aaaTGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAGTCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGGTTTTTCAGGTTATCCAGCGCTCTGGCGGCGGTACCGTCGGTCCGCTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGatata                                                                                                                                                                                                                   <  2:226724/251‑1 (MQ=255)            gAATTTCGCCCACGTGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGCTAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTc                                                                                                                                                                                                               >  1:209938/1‑251 (MQ=255)                          tGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAA‑GCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGc                                                                                                                                                                                                   >  1:194878/1‑248 (MQ=255)                                 cgtcgGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGGTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGAGTTTTAAATCGGCCAGAACCAGATCGATATAGGGCAGATAAGGGGAGCTATAATTCCAAGGCACATGCAGACAGGtt                                                                                                                                                                                          >  2:185528/1‑251 (MQ=255)                                                          ttttACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGCGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATgcctgcct                                                                                                                                                                 >  2:121036/1‑251 (MQ=255)                                                                                    aGCCCTGAATCCGCGGCCCGCGGCTAATGATTTTTTTGCCCGCTGCGCCGAGTTTTTTCAGGTTATCCAGCTCTCTGGCGGCGTTACCCTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGCTAAGGGGCGATATATTTCCACGGCACATGCAGACAGCTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGGAGcgccatc                                                                                                                                       >  2:12380/1‑251 (MQ=255)                                                                                      cccTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTGTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgc                                                                                                                                     <  1:223935/251‑1 (MQ=255)                                                                                                 cGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGACGGTCCACTGTTTAAACGGCGCGTCGGCACCGTGTTTTAAATCGGCACGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCt                                                                                                                          >  1:226405/1‑251 (MQ=255)                                                                                                             aaTGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGGCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCTCCATCGCCATTTCCGTCTTCATAACGGGCt                                                                                                              <  2:98573/251‑1 (MQ=255)                                                                                                             aaTGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTACCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCt                                                                                                              <  2:194878/251‑1 (MQ=255)                                                                                                                                        ttttttCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAccgccgc                                                                                   >  1:138287/1‑251 (MQ=255)                                                                                                                                                ggTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCg                                                                            >  1:189754/1‑250 (MQ=255)                                                                                                                                                 gTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                                                                          <  2:67545/251‑1 (MQ=255)                                                                                                                                                                                   cACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCAt                                        <  2:89031/251‑1 (MQ=255)                                                                                                                                                                                            aaCGGAGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATGATctcctc                               >  1:167914/1‑251 (MQ=255)                                                                                                                                                                                                                  ttAAATCGGCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATGATCTCCTCAACGCTTTTCACTTCACcacac         >  2:149139/1‑251 (MQ=255)                                                                                                                                                                                                                          gCAAGAAACAGATCGATATAGGGCAGATAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAACGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATGATCTCCTCAACGCTTTTCACTTCACCACACACGGTTa  >  2:199134/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               GGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATGATCTCCTCAACGCTTTTCACTTCACCACACACGGTTA  >  NC_000913/862756‑863222

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: ‑