Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,340,898	C→T	R244H (CGC→CAC)	melR ←	melibiose operon transcriptional regulator, autoregulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,340,898	0	C	T	100.0%	67.1 / NA	20	R244H (CGC→CAC)	melR	melibiose operon transcriptional regulator, autoregulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/11);  total (9/11)

TGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATA  >  NC_000913/4340766‑4341032                                                                                                                                     |                                                                                                                                       tgttgCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACGTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcg                                                                                                                                  <  1:261356/139‑1 (MQ=255)   ttgCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcgcg                                                                                                                                >  1:184025/1‑139 (MQ=255)           aCATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTa                                                                                                                             <  2:456577/134‑1 (MQ=255)           aCATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTa                                                                                                                             >  1:456577/1‑134 (MQ=255)                               aaCGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCa                                                                                                    <  1:358104/139‑1 (MQ=255)                                 cGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATg                                                                                                  <  2:383226/139‑1 (MQ=255)                                   tGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcg                                                                                                                                  <  2:82544/105‑1 (MQ=255)                                   tGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcg                                                                                                                                  >  1:82544/1‑105 (MQ=255)                                               gACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAatat                                                                                    <  1:6618/139‑1 (MQ=255)                                                     ttGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCAt                                                                              <  1:140889/139‑1 (MQ=255)                                                                                     agaATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAg                                              >  1:290790/1‑139 (MQ=255)                                                                                          aCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcg                                                                                                                                  <  1:452499/50‑1 (MQ=255)                                                                                          aCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATcg                                                                                                                                  >  2:452499/1‑50 (MQ=255)                                                                                                      aTCGCTCAGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTc                             >  1:391107/1‑139 (MQ=255)                                                                                                             aGTAACGCGCGAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCt                      <  2:184025/139‑1 (MQ=255)                                                                                                                   gcgcgAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCAt                >  2:204929/1‑139 (MQ=255)                                                                                                                    cgcgAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTg                                                                            <  2:366530/78‑1 (MQ=255)                                                                                                                    cgcgAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTg                                                                            >  1:366530/1‑78 (MQ=255)                                                                                                                    cgcgAACGTGGTTGATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATa               >  1:368743/1‑139 (MQ=255)                                                                                                                                 gATGTGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATa  <  2:40134/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATGACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGTCAGCGCCTGATCATAGTTTTCGGCAATA  >  NC_000913/4340766‑4341032

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