Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,940,785	T→C	K473K (AAA→AAG)	bcsG ←	predicted inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,940,785	0	T	C	100.0%	113.4 / NA	34	K473K (AAA→AAG)	bcsG	predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/20);  total (14/20)

GCCAGGAAGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAATTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATACCTGCATTC  >  W3110S.gb/3940717‑3940852                                                                     |                                                                    gCCAGGAAGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAActt                                                                   <  1:451172/71‑1 (MQ=255) gCCAGGAAGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAActt                                                                   <  1:2531307/71‑1 (MQ=255)        aGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCAcccc                                                             <  1:3774016/70‑1 (MQ=255)        aGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCAcccc                                                             >  1:1464934/1‑70 (MQ=255)            gctCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACgg                                                        >  1:4352019/1‑71 (MQ=255)             ctCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACggg                                                       <  1:1127066/71‑1 (MQ=255)             ctCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACggg                                                       <  1:4341154/71‑1 (MQ=255)             ctCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACggg                                                       <  1:2403501/71‑1 (MQ=255)                       cGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGt                                              <  1:2031149/70‑1 (MQ=255)                          tGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:4085866/71‑1 (MQ=255)                          tGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:3448162/71‑1 (MQ=255)                          tGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:1758168/71‑1 (MQ=255)                          tGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:1730739/71‑1 (MQ=255)                                  ggTGCCCCCTGATGCGGGGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:1812384/63‑1 (MQ=255)                                   gTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCt                                  <  1:3390352/70‑1 (MQ=255)                                    tGCCCCCTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAg                                >  1:766369/1‑71 (MQ=255)                                        cccTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGat                            <  1:364884/71‑1 (MQ=255)                                          cTGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACggg                                                       <  1:2174514/42‑1 (MQ=255)                                           tGATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGatat                          >  1:1026399/1‑70 (MQ=255)                                            gATGCGGTGCCTTCATGCCGAAGAACTTCACCCCAACAGGGACGTCGGTGATAGACGGGCTAGGGATATCa                        >  1:730757/1‑71 (MQ=255)                                                     ccTTCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGc                 >  1:3088565/1‑69 (MQ=255)                                                       ttCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGt                                              <  1:4613343/38‑1 (MQ=255)                                                       ttCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATa                                          <  1:2329289/42‑1 (MQ=255)                                                       ttCATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGat                            <  1:86500/56‑1 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGa             <  1:2839390/69‑1 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGAt            <  1:992932/70‑1 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGAt            <  1:1745245/70‑1 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:282798/1‑71 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:255296/1‑71 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:1772029/1‑71 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:1415376/1‑71 (MQ=255)                                                         cATGCCGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:872677/1‑71 (MQ=255)                                                         cATGCCGAAGAACTTCACCCAAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATa           >  1:3160733/1‑71 (MQ=255)                                                             ccGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTa                    >  1:1110011/1‑58 (MQ=255)                                                             ccGAAGAACTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATACCTg       <  1:4280198/71‑1 (MQ=255)                                                                  gaaCTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATACCTGCATTc  >  1:2759828/1‑71 (MQ=255)                                                                     |                                                                    GCCAGGAAGCTGCTCGGTTGTTCGATGACAATCGGTGCCCCCTGATGCGGTGCCTTCATGCCGAAGAATTTCACCCCAACGGGGACGTCGGTGATAGACGGGCTAGGGATATCACGTAGGCCAGATACCTGCATTC  >  W3110S.gb/3940717‑3940852

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