Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,094,703	C→T	intergenic (‑431/+20)	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,703	0	C	T	93.8%	40.3 / ‑3.6	16	intergenic (‑431/+20)	rhaD/rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (9/6);  total (10/6)
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.67e-01

GCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCC  >  NC_000913/4094624‑4094832                                                                                |                                                                                                                                  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)                                                                     gATAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACACTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGc   <  1:504823/140‑1 (MQ=255)                                                                      aTAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGcc  >  2:188230/1‑140 (MQ=255)                                                                      aTAAGGCACTTGTGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACATTGCCAGTACATTTCCCAGACCGCCTGcc  >  1:253240/1‑140 (MQ=255)                                                                                |                                                                                                                                  GCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCC  >  NC_000913/4094624‑4094832

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: ‑