Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,646,699	C→G	Q1415H (CAG→CAC)	yfhM ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,646,699	0	C	G	100.0%	198.0 / NA	57	Q1415H (CAG→CAC)	yfhM	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (21/36);  total (21/36)

TCGTACGGGCAGCCAGGAACAACGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGA  >  W3110S.gb/2646632‑2646766                                                                    |                                                                    tCGTACGGGCAGCCAGGAACAACGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCctggtg                                                                   >  1:4041403/1‑70 (MQ=255) tCGTACGGGCAGCCAGGAACAACGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCtggtgc                                                                  >  1:406403/1‑71 (MQ=255)              cAGGAACAACGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGtt                                                       >  1:1017868/1‑69 (MQ=255)                  aacaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:4173313/1‑70 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:1846910/1‑69 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:585517/1‑69 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:634518/1‑69 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:737827/1‑69 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:935880/1‑69 (MQ=255)                   acaacGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCa                                                  >  1:1006011/1‑69 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCTTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:4300880/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:797382/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:516458/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:460126/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:997773/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:966264/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:4048888/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:3860822/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:3801524/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:3586493/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:3229554/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:3212024/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:1492791/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:1715721/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:283688/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:2190806/71‑1 (MQ=255)                      acGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:2518781/71‑1 (MQ=255)                       cgcgTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtgt                                            <  1:2042211/71‑1 (MQ=255)                       cgcgTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:1199811/70‑1 (MQ=255)                           ttAGTTTGCTGGGTCGATAGCCAGCGTTCACGAAACGCCTGGTGCGAAAGTGTGTTCAGCATAgtgt                                            <  1:2396061/67‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGaa                                                               <  1:736371/48‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGaa                                                               <  1:1675269/48‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGaa                                                               <  1:938885/48‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGaa                                                               <  1:2663567/48‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGaa                                                               <  1:944327/48‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtgt                                            <  1:2502495/67‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:4222183/66‑1 (MQ=255)                           ttACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtg                                             <  1:317511/66‑1 (MQ=255)                            tACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAgtgt                                            <  1:346958/66‑1 (MQ=255)                            tACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTAc                                          >  1:978141/1‑68 (MQ=255)                             aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGc                                       >  1:1501098/1‑70 (MQ=255)                             aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCt                                      >  1:3722500/1‑71 (MQ=255)                             aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCt                                      >  1:2800025/1‑71 (MQ=255)                             aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCt                                      >  1:2837393/1‑71 (MQ=255)                             aCTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCt                                      >  1:2842774/1‑71 (MQ=255)                                             aGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGtt                       >  1:3916963/1‑70 (MQ=255)                                             aGCCAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGtt                       >  1:1823461/1‑70 (MQ=255)                                                cAGCGTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTAtt                   <  1:413253/71‑1 (MQ=255)                                                    gTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTttct               <  1:4544212/71‑1 (MQ=255)                                                    gTTCACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTttct               >  1:3274269/1‑71 (MQ=255)                                                       cACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTg                                        >  1:120724/1‑43 (MQ=255)                                                       cACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTAt                    >  1:1153252/1‑63 (MQ=255)                                                       cACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTttcttc             >  1:2637997/1‑70 (MQ=255)                                                       cACCAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTttcttc             >  1:3163153/1‑70 (MQ=255)                                                          cAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc         <  1:894986/71‑1 (MQ=255)                                                          cAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc         <  1:519572/71‑1 (MQ=255)                                                          cAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc         <  1:3003196/71‑1 (MQ=255)                                                          cAAACGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGc         <  1:28776/71‑1 (MQ=255)                                                           aaaCGCCTGGTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCa        <  1:1411636/71‑1 (MQ=255)                                                                 ctggtgCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGgaga  <  1:4610075/71‑1 (MQ=255)                                                                    |                                                                    TCGTACGGGCAGCCAGGAACAACGCGTTACTTTCCTGCGTCGATAGCCAGCGTTCACCAAACGCCTGCTGCGAAAGTGTGTTCAGCAAAGTGTACTGCTCATCGGGTAGCAGTTTATTTTCTTCCAGCAAGGAGA  >  W3110S.gb/2646632‑2646766

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