Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,582,992	T→C	S81S (TCA→TCG)	ligA ←	DNA ligase, NAD(+)‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,582,992	0	T	C	100.0%	202.6 / NA	55	S81S (TCA→TCG)	ligA	DNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/39);  total (16/39)

TTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTT  >  minE/1582924‑1583056                                                                     |                                                                 tttttCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGa                                                                <  1:2911037/71‑1 (MQ=255)  ttttCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGa                                                                <  1:613552/70‑1 (MQ=255)     tCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAg                                                             <  1:2747720/70‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:3098928/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:2905770/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:3819178/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:3798278/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:1856667/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:1947393/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:1973869/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:4176862/71‑1 (MQ=255)        gACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCAtt                                                         <  1:2358757/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:3213542/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:3908865/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:3937309/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:434101/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:4640108/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:1838700/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:4840591/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:1449053/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:3278366/71‑1 (MQ=255)          cGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:855040/71‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTgg                                                       <  1:768032/70‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGt                                                      <  1:3731969/71‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGt                                                      <  1:3435736/71‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGt                                                      <  1:2850268/71‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGt                                                      <  1:4670592/71‑1 (MQ=255)           ggTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGt                                                      <  1:13024/71‑1 (MQ=255)                      cGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCg                                           >  1:896827/1‑71 (MQ=255)                      cGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCg                                           >  1:947659/1‑71 (MQ=255)                      cGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCg                                           >  1:3264954/1‑71 (MQ=255)                      cGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCg                                           >  1:2147337/1‑71 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:1173936/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:1214498/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:2431662/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:636602/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:1364799/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:2003723/53‑1 (MQ=255)                       gTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGCGACAGc                                                            <  1:4647899/53‑1 (MQ=255)                                              tcttcATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAgc                     <  1:1825727/69‑1 (MQ=255)                                              tcttcATCAAAAACGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAgc                     <  1:4693749/69‑1 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:1051883/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:3509805/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:1221211/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:659027/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:1696824/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:4777286/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:450151/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:427167/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:4263459/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:4094017/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc           >  1:3045838/1‑71 (MQ=255)                                                      aaaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcg            >  1:321032/1‑70 (MQ=255)                                                       aaaaCGTTATCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGcc          <  1:3346540/71‑1 (MQ=255)                                                                tCCAGCGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGtt  <  1:3274935/70‑1 (MQ=255)                                                                     |                                                                 TTTTTCAGACGGTCCTGCACACGTTTGTTGAAAGCAAGAAAGCTTTCTTCATCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTT  >  minE/1582924‑1583056

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