Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,456,196	G→A	A385V (GCG→GTG)	yfaQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,456,196	0	G	A	100.0%	192.5 / NA	62	A385V (GCG→GTG)	yfaQ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/36);  total (26/36)

ATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGG  >  minE/1456135‑1456258                                                              |                                                               aTGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTg                                                        <  1:4249142/70‑1 (MQ=255)   gATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTc                                                     >  1:352622/1‑71 (MQ=255)   gATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTc                                                     >  1:2854438/1‑71 (MQ=255)    aTAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCg                                                    >  1:1399211/1‑71 (MQ=255)    aTAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCg                                                    >  1:3818155/1‑71 (MQ=255)    aTAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCg                                                    >  1:7583/1‑71 (MQ=255)     tAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGc                                                   <  1:2941973/71‑1 (MQ=255)              ggATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGc                                          <  1:3135314/71‑1 (MQ=255)              ggATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGc                                          <  1:4597983/71‑1 (MQ=255)              ggATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGc                                          <  1:2906082/71‑1 (MQ=255)               gATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGc                                          <  1:2578007/70‑1 (MQ=255)                 tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:2050712/71‑1 (MQ=255)                 tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:3941089/71‑1 (MQ=255)                 tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:4174513/71‑1 (MQ=255)                 tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:1974335/71‑1 (MQ=255)                 tCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:1731974/71‑1 (MQ=255)                  cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       >  1:1290199/1‑70 (MQ=255)                  cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       <  1:3505267/70‑1 (MQ=255)                  cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       >  1:4462294/1‑70 (MQ=255)                  cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       >  1:2088006/1‑70 (MQ=255)                  cGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGa                                       >  1:1365375/1‑70 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAacac                                  >  1:4061893/1‑71 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAacac                                  >  1:3810752/1‑71 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAacac                                  >  1:1728785/1‑71 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAacac                                  >  1:3215506/1‑71 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAaca                                   >  1:1976188/1‑70 (MQ=255)                      tGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAaca                                   >  1:3612641/1‑70 (MQ=255)                       gCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACg                                 >  1:825065/1‑71 (MQ=255)                       gCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACg                                 >  1:3146552/1‑71 (MQ=255)                                      tGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGaa                                      <  1:4623098/51‑1 (MQ=255)                                      tGCGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGaa                                      <  1:2597873/51‑1 (MQ=255)                                        cGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGctgctg                 <  1:3108112/70‑1 (MQ=255)                                        cGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGctgctg                 <  1:2331102/70‑1 (MQ=255)                                        cGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGctgctg                 <  1:3949534/70‑1 (MQ=255)                                        cGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGctgctg                 <  1:3873129/70‑1 (MQ=255)                                        cGTCCAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGctgctg                 <  1:3786585/70‑1 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:344821/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:1954239/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:4765926/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            <  1:2282877/71‑1 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            <  1:2617066/71‑1 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            <  1:2697422/71‑1 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:3662322/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            <  1:3633368/71‑1 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:3500012/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            >  1:3578608/1‑71 (MQ=255)                                            cAGGCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCggg            <  1:3404460/71‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:4369442/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:589336/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:687045/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:4879381/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:3502310/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:4599595/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:1175869/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:4232626/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:2991762/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:3649998/62‑1 (MQ=255)                                               gCGGTCATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:2705308/62‑1 (MQ=255)                                                    cATCGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCgctgct                  <  1:3695276/57‑1 (MQ=255)                                                       cGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAgg  >  1:586066/1‑70 (MQ=255)                                                       cGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAgg  >  1:2084750/1‑70 (MQ=255)                                                       cGTTCGCACCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAgg  >  1:2095608/1‑70 (MQ=255)                                                              |                                                               ATGATAATGAACCGGATCGCCTGCATAGATGAGATCCTGCGTCCAGGCGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGG  >  minE/1456135‑1456258

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