Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,109,964	G→A	C185C (TGC→TGT)	mdoC ←	membrane protein required for modification of periplasmic glucan

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,109,964	0	G	A	100.0%	77.3 / NA	26	C185C (TGC→TGT)	mdoC	membrane protein required for modification of periplasmic glucan
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/9);  total (17/9)

TGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAAT  >  W3110S.gb/1109909‑1110031                                                        |                                                                    tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:1640990/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:1654804/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:77810/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:3626028/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:2372225/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:3024814/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:244715/1‑71 (MQ=255) tGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      >  1:2833860/1‑71 (MQ=255)      ggATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCTAGACATAAAAAAATCACCGATAg                                                 <  1:787581/71‑1 (MQ=255)      ggATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAg                                                 <  1:1053191/71‑1 (MQ=255)      ggATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAg                                                 <  1:2861392/71‑1 (MQ=255)           cacaATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTa                                            >  1:1096719/1‑71 (MQ=255)                  aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAAt                                                          >  1:2819978/1‑50 (MQ=255)                  aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAAt                                                          >  1:2646622/1‑50 (MQ=255)                  aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAAt                                                          >  1:469483/1‑50 (MQ=255)                  aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATc                                       >  1:58667/1‑69 (MQ=255)                  aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATc                                       >  1:1380233/1‑69 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:252944/70‑1 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:3119064/70‑1 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:3209848/70‑1 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:4121897/70‑1 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:2201445/70‑1 (MQ=255)                        gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAtt                                <  1:795465/70‑1 (MQ=255)                                            aCCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAAt           >  1:294773/1‑71 (MQ=255)                                            aCCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAAt           >  1:2438637/1‑71 (MQ=255)                                            aCCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAAt           >  1:2356072/1‑71 (MQ=255)                                                     gagaCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAat  >  1:3397335/1‑71 (MQ=255)                                                        |                                                                    TGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAAT  >  W3110S.gb/1109909‑1110031

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