Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	692,900	T→C	A144A (GCT→GCC)	hyaD →	protein involved in processing of HyaA and HyaB proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	692,900	0	T	C	100.0%	209.6 / NA	56	A144A (GCT→GCC)	hyaD	protein involved in processing of HyaA and HyaB proteins
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/28);  total (28/28)

CAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAT  >  minE/692835‑692960                                                                  |                                                             cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCgg                                                          >  1:3588778/1‑70 (MQ=255) cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCgg                                                          >  1:1007355/1‑70 (MQ=255) cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCgg                                                          >  1:2676690/1‑70 (MQ=255) cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCgg                                                          >  1:797755/1‑70 (MQ=255) cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCgg                                                          >  1:590687/1‑70 (MQ=255) cAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGCGAGCAACTGCCCGCCGCgg                                                          >  1:996247/1‑70 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAggcggc                                                >  1:436194/1‑70 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:4085756/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:318867/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2684216/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2654468/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2560902/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2381224/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2300468/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:1238391/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:1908114/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:701078/1‑71 (MQ=255)           tGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGgcggcg                                               >  1:2162526/1‑71 (MQ=255)                    acTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:1140867/71‑1 (MQ=255)                    acTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:2828636/71‑1 (MQ=255)                    acTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:2633202/71‑1 (MQ=255)                    acTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:2628071/71‑1 (MQ=255)                     cTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:917242/70‑1 (MQ=255)                     cTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAg                                      <  1:3629627/70‑1 (MQ=255)                        cggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCtt                                   <  1:3757519/70‑1 (MQ=255)                        cggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCtt                                   <  1:3511258/70‑1 (MQ=255)                        cggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCtt                                   <  1:3435563/70‑1 (MQ=255)                        cggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCtt                                   <  1:3006936/70‑1 (MQ=255)                        cggcggTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCtt                                   <  1:2397176/70‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:3736892/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:999283/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:426562/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:2925407/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:1976977/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:2178860/71‑1 (MQ=255)                                     aGCGAACTGGCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAAt                     <  1:2359207/71‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCTTGGGGAATTGTGCCGc             <  1:3450779/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTTCCGc             <  1:247554/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc             <  1:1257848/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc             <  1:82835/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc             <  1:2889631/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc             <  1:2256115/70‑1 (MQ=255)                                              gCACGGGAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGc             <  1:403073/70‑1 (MQ=255)                                                   ggAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCgg       >  1:3364071/1‑71 (MQ=255)                                                    gAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc      >  1:4254601/1‑71 (MQ=255)                                                    gAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc      >  1:2533672/1‑71 (MQ=255)                                                    gAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc      >  1:614883/1‑71 (MQ=255)                                                    gAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc      >  1:2796474/1‑71 (MQ=255)                                                    gAGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGc      >  1:2801499/1‑71 (MQ=255)                                                     aGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCt     <  1:2917319/71‑1 (MQ=255)                                                     aGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCt     <  1:1838352/71‑1 (MQ=255)                                                     aGCAACTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCt     <  1:1815915/71‑1 (MQ=255)                                                         aCTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAt  >  1:2433759/1‑70 (MQ=255)                                                         aCTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAt  >  1:4112316/1‑70 (MQ=255)                                                         aCTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAt  >  1:848108/1‑70 (MQ=255)                                                         aCTGCCCGCCGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAt  >  1:3110874/1‑70 (MQ=255)                                                                  |                                                             CAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAACCGGCTAAT  >  minE/692835‑692960

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