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	100.0%	110.7 / NA	35	*255* (TAG→TAA)	pdhR	pyruvate dehydrogenase complex repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/20);  total (15/20)

CGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122745‑122989                                                                                                                |                                                                                                                                      cgAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                                            >  1:186156/1‑139 (MQ=255)                ccATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          >  1:345274/1‑126 (MQ=255)                ccATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          <  2:345274/126‑1 (MQ=255)                 cATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                                            <  1:320662/139‑1 (MQ=255)                 cATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                                            <  2:330668/139‑1 (MQ=255)                   tCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                                          <  2:568354/139‑1 (MQ=255)                      ggCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCa                                                                                       <  1:130278/139‑1 (MQ=255)                            tATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                                          <  2:740634/130‑1 (MQ=255)                            tATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                                          >  1:740634/1‑130 (MQ=255)                                        tttGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACt                                                                     <  2:356894/139‑1 (MQ=255)                                        tttGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACt                                                                     <  1:491309/139‑1 (MQ=255)                                                   gAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                                          <  2:156/139‑1 (MQ=255)                                                      gTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTg                                                       <  2:186156/139‑1 (MQ=255)                                                                gagCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTg                                                       <  2:572322/129‑1 (MQ=255)                                                                gagCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTg                                                       >  1:572322/1‑129 (MQ=255)                                                                 agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACt                                                                     >  2:105788/1‑114 (MQ=255)                                                                 agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACt                                                                     <  1:105788/114‑1 (MQ=255)                                                                 agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  1:98163/1‑139 (MQ=255)                                                                   ccgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCtt                                          >  2:660688/1‑139 (MQ=255)                                                                           aGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCga                                  <  2:311493/139‑1 (MQ=255)                                                                                        tcgtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTa                                                                    >  2:513208/1‑92 (MQ=255)                                                                                        tcgtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTa                                                                    <  1:513208/92‑1 (MQ=255)                                                                                        tcgtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGa                     <  1:44845/139‑1 (MQ=255)                                                                                         cgtcTGGAGCAACGAAAGAAATAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAc                    <  2:294459/139‑1 (MQ=255)                                                                                          gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAGTTATCCagaag                                                                                                          <  2:81778/52‑1 (MQ=255)                                                                                          gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAGTTATCCagaag                                                                                                          >  1:81778/1‑52 (MQ=255)                                                                                          gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACa                   >  2:703826/1‑139 (MQ=255)                                                                                           tcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAg                  >  2:157694/1‑139 (MQ=255)                                                                                              ggAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGtt               >  2:167864/1‑139 (MQ=255)                                                                                                     cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGaaa        >  2:694792/1‑139 (MQ=255)                                                                                                     cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGaaa        >  1:137293/1‑139 (MQ=255)                                                                                                          gAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCa   >  1:444283/1‑139 (MQ=255)                                                                                                           aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  1:737535/139‑1 (MQ=255)                                                                                                           aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  2:98163/139‑1 (MQ=255)                                                                                                           aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  1:660688/139‑1 (MQ=255)                                                                                                                |                                                                                                                                      CGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122745‑122989

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