Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,094,691	G→T	intergenic (‑419/+32)	rhaD ← / ← rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,094,691	0	G	T	100.0%	34.1 / NA	13	intergenic (‑419/+32)	rhaD/rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (8/5);  total (8/5)

GCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCC  >  NC_000913/4094624‑4094820                                                                    |                                                                                                                                  gCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCa                                                           <  2:779709/140‑1 (MQ=17)              gATGCGACTCTGTCGCTCCTTAACAGGCCTACATCTATTGCCGACCTGAAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCa                                              <  2:405185/140‑1 (MQ=18)                tGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTCGTGCCGCATCCGGCAGTGTTTACCCGGGGCGACTCTAAAATTCGTTCTCATAAGCCCGCACGCTCTCCAGCCAtt                                            >  2:236184/1‑140 (MQ=18)                    aCGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCt                                        >  1:788809/1‑140 (MQ=21)                     cGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTa                                       >  1:41980/1‑140 (MQ=21)                          tCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctgc                                  >  2:13138/1‑140 (MQ=21)                                        gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGAc                    <  2:788809/140‑1 (MQ=37)                                             cAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTgg               >  2:67500/1‑140 (MQ=37)                                              caactttGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGc              <  1:437589/135‑1 (MQ=38)                                               aCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCa             <  2:510685/140‑1 (MQ=37)                                                  attGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa          >  1:600276/2‑140 (MQ=38)                                                        gACCTGTAGGCCTGATAAGACACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTc    >  2:481342/1‑140 (MQ=38)                                                          cctgTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTccc  >  2:751351/4‑140 (MQ=38)                                                                    |                                                                                                                                  GCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCC  >  NC_000913/4094624‑4094820

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