Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	122,856	G→A	*255* (TAG→TAA)	pdhR →	pyruvate dehydrogenase complex repressor, autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	122,856	0	G	A	94.7%	50.3 / ‑3.8	19	*255* (TAG→TAA)	pdhR	pyruvate dehydrogenase complex repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAC  >  NC_000913/122747‑122971                                                                                                              |                                                                                                                    aaGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCagaag                                                                                        <  1:97987/139‑1 (MQ=255)               cATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                          <  2:691168/139‑1 (MQ=255)                   tGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTATGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgcgc                                                                      >  1:101440/1‑139 (MQ=255)                      ccTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                        >  1:724947/1‑134 (MQ=255)                      ccTTTATCGAAGAAATTTTGCACGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                        <  2:724947/134‑1 (MQ=255)                             cGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAgc                                                            >  2:173437/1‑139 (MQ=255)                                                 gAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                        <  1:328654/139‑1 (MQ=255)                                                            gagagCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTg                             >  1:360467/1‑139 (MQ=255)                                                               agCCGCCGAGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                          >  2:353002/1‑75 (MQ=255)                                                               agCCGCCGAGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                          <  1:353002/75‑1 (MQ=255)                                                                  cgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTaaaaa                                                                                                   >  1:826996/1‑63 (MQ=255)                                                                  cgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTaaaaa                                                                                                   <  2:826996/63‑1 (MQ=255)                                                                             ttCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGt            <  2:157841/139‑1 (MQ=255)                                                                              tctcTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                          <  1:679878/76‑1 (MQ=255)                                                                              tctcTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                          >  2:679878/1‑76 (MQ=255)                                                                              tctcTGCGTCGGCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGtt           <  1:542930/139‑1 (MQ=255)                                                                                       cgtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAc  <  2:360467/139‑1 (MQ=255)                                                                                                   cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                        <  1:583393/89‑1 (MQ=255)                                                                                                   cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                        >  2:583393/1‑89 (MQ=255)                                                                                                              |                                                                                                                    AAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAC  >  NC_000913/122747‑122971

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