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

TAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122727‑122989                                                                                                                                  |                                                                                                                                      tAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTc                                                                                                                              >  1:181966/1‑139 (MQ=255)        gaagaaGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTaaa                                                                                                                       >  2:280647/1‑139 (MQ=255)                   cgAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                                            <  1:1216673/139‑1 (MQ=255)                     aaGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          <  1:5900/139‑1 (MQ=255)                                cgcCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          >  2:73559/1‑128 (MQ=255)                                cgcCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          <  1:73559/128‑1 (MQ=255)                                 gcCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATc                                                                                              <  2:1109898/139‑1 (MQ=255)                                   cATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaa                                                                                            <  2:181966/139‑1 (MQ=255)                                       tGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgcgc                                                                                        <  2:474636/139‑1 (MQ=255)                                             ttATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATaa                                                                                  >  2:843851/1‑139 (MQ=255)                                                   aagaaATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAgcgc                                                                            >  2:846825/1‑139 (MQ=255)                                                          tttGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACt                                                                     <  2:838244/139‑1 (MQ=255)                                                              cTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGt                                                                                                                          <  1:1062904/82‑1 (MQ=255)                                                              cTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGt                                                                                                                          >  2:1062904/1‑82 (MQ=255)                                                               tCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACg                                                                <  1:350232/139‑1 (MQ=255)                                                                     gAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                                          <  1:846825/139‑1 (MQ=255)                                                                          cGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTc                                                     <  1:843851/139‑1 (MQ=255)                                                                          cGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTc                                                     <  2:882193/139‑1 (MQ=255)                                                                                   agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          >  1:618790/1‑77 (MQ=255)                                                                                   agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          <  2:618790/77‑1 (MQ=255)                                                                                   agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  1:215928/1‑139 (MQ=255)                                                                                   agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  2:330748/1‑139 (MQ=255)                                                                                                 ttCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGc                                    <  1:8449/133‑1 (MQ=255)                                                                                                 ttCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGc                                    >  2:8449/1‑133 (MQ=255)                                                                                                           cgtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAc                    <  2:462130/139‑1 (MQ=255)                                                                                                           cgtcTGGAGCAACGAAAGAATTAATGATTTTTCTAGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGAc                    <  1:280647/139‑1 (MQ=255)                                                                                                            gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          >  1:814682/1‑52 (MQ=255)                                                                                                            gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCagaag                                                                                                          <  2:814682/52‑1 (MQ=255)                                                                                                             tcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTGAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTc                            >  2:931474/1‑129 (MQ=255)                                                                                                             tcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTc                            <  1:931474/129‑1 (MQ=255)                                                                                                             tcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAg                  >  2:635252/1‑139 (MQ=255)                                                                                                                      aCGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGAtgttgt                                                                                                   <  2:1433108/49‑1 (MQ=255)                                                                                                                      aCGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGAtgttgt                                                                                                   >  1:1433108/1‑49 (MQ=255)                                                                                                                       cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAAAGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGaaa        >  1:666430/1‑139 (MQ=255)                                                                                                                            gAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaggc                                                                                          <  2:110734/52‑4 (MQ=255)                                                                                                                            gAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCaggc                                                                                          >  1:110734/1‑49 (MQ=255)                                                                                                                            gAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCa   >  2:386260/1‑139 (MQ=255)                                                                                                                             aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  2:468956/139‑1 (MQ=255)                                                                                                                             aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  2:672631/139‑1 (MQ=255)                                                                                                                             aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  1:383998/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      TAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122727‑122989

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