Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,062,561	C→T	L248F (CTT→TTT)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,062,561	0	C	T	100.0%	110.1 / NA	35	L248F (CTT→TTT)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/22);  total (13/22)

CCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATG  >  minE/2062499‑2062623                                                               |                                                               cccccAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGa                                                         <  1:133027/70‑1 (MQ=255)         tACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCgcagca                                                <  1:2062908/71‑1 (MQ=255)         tACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCgcagca                                                <  1:1633240/71‑1 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGtttt                                                              <  1:3592131/56‑1 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGtttt                                                              <  1:2832540/56‑1 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGtttt                                                              <  1:4897374/56‑1 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGcagcag                                               >  1:29588/1‑71 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGcagcag                                               >  1:2686336/1‑71 (MQ=255)          aCCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGcagcag                                               >  1:4360134/1‑71 (MQ=255)             tGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGa                                             >  1:2018372/1‑70 (MQ=255)             tGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGa                                             >  1:2473364/1‑70 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:1602060/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:3511340/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:3211108/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:1205356/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:3965371/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:2134005/71‑1 (MQ=255)                         aaaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAAccc                                <  1:2049649/71‑1 (MQ=255)                          aaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCGGCAGGATAACGCTGAACCCa                               <  1:550521/71‑1 (MQ=255)                          aaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCa                               <  1:2001022/71‑1 (MQ=255)                          aaCAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCa                               <  1:1389667/71‑1 (MQ=255)                                  ttttCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGcagcag                                               <  1:2929285/47‑1 (MQ=255)                                  ttttCCATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGcagcag                                               <  1:1626903/47‑1 (MQ=255)                                       cATAAAGGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTtgt                  <  1:3657712/71‑1 (MQ=255)                                       cATAAAAGCCGCCAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTtgt                  <  1:921103/71‑1 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGCAGGCta      <  1:6068/71‑1 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCta      >  1:2805576/1‑71 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCta      <  1:21659/71‑1 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCta      <  1:1374883/71‑1 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCta      <  1:87209/71‑1 (MQ=255)                                                   cAGCTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCta      <  1:1315361/71‑1 (MQ=255)                                                      cTTTACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTatat   >  1:3561069/1‑71 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTTCTTGTGGTCGAAGGCTATATg  >  1:2139327/1‑71 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTatat   >  1:4531267/1‑70 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATg  >  1:3913414/1‑71 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATg  >  1:4564069/1‑71 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATg  >  1:4808135/1‑71 (MQ=255)                                                       tttACGGTTTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATg  >  1:1569714/1‑71 (MQ=255)                                                               |                                                               CCCCCAAATACCTGAACTCGCCGGAAACAGACATTTTCCATAAAGGCCGCCAGCTTTACGGTCTTTATGAAGCGCAGCAGGATAACGCTGAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATG  >  minE/2062499‑2062623

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