Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	REL606	226,159	C→A	100%	P163T (CCT→ACT)	gmhB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	226,159	0	C	A	100.0%	72.2 / NA	22	P163T (CCT→ACT)	gmhB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/3);  total (19/3)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTCCTGAG  >  REL606/226039‑226301                                                                                                                         |                                                                                                                                               tattTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAGTGGgtg                                                                                                                   >  1:666548/1‑150 (MQ=255)      gCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTaaa                                                                                                              <  2:668516/150‑1 (MQ=255)       gatatTGATATGGCAGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGAGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAAt                                                                                                             <  1:22137/149‑1 (MQ=255)          atTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGc                                                                                                          >  1:562430/1‑150 (MQ=255)           ttGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGcc                                                                                                         >  2:56209/1‑150 (MQ=255)                         tcttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCg                                                                                           >  1:227248/1‑150 (MQ=255)                          cttATATGGTGGGCGATATATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCTGATTGGGTGTTAAATAGCCTGACAGACCTGCCGc                                                                                          >  1:361351/1‑150 (MQ=255)                          cttATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGc                                                                                          >  2:426620/1‑150 (MQ=255)                          cttATATGGCGGGCGATAAATTAGAAGATATGAAGGCAGCGGCTGCGGCGAACGTTGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTAAGATCGAAGCAGAAAACGGAGAGGATTGGGCGTTAACTAGCCTGTCAAACCTGCAGc                                                                                          >  1:571435/1‑150 (MQ=255)                              tatGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGc                                                                                      <  1:56043/150‑1 (MQ=255)                                                    gaTATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCg                                                                >  1:280267/1‑150 (MQ=255)                                                    gaTATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAATGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCg                                                                >  1:108512/1‑150 (MQ=255)                                                       atGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCa                                                             >  2:528577/1‑150 (MQ=255)                                                                  ggCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTaaa                                                  >  2:335394/1‑150 (MQ=255)                                                                                   gAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAAt                                 >  2:193924/1‑150 (MQ=255)                                                                                                        cGGGTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTg            >  1:377591/1‑150 (MQ=255)                                                                                                          ggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTc          >  2:129121/1‑150 (MQ=255)                                                                                                          ggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTc          >  2:202954/1‑150 (MQ=255)                                                                                                          ggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTc          >  2:428339/1‑150 (MQ=255)                                                                                                          ggTAAACCTATTACGACTGAAGCAGAAAACGCGGCGGAGTGGGTGTTAAATAGCCTGGCAGATCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATAATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCAACTAGTc          >  1:662775/1‑150 (MQ=255)                                                                                                              aaCCTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGGCCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTcc      <  1:505243/150‑1 (MQ=255)                                                                                                                ccTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTCCTg    >  2:313423/1‑150 (MQ=255)                                                                                                                ccTATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTCCTg    >  1:277395/1‑150 (MQ=255)                                                                                                                  tATTACGACTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTCCTgag  >  2:617432/1‑150 (MQ=255)                                                                                                                         |                                                                                                                                               TATTTGCATATTGATATGGCCGCTTCTTATATGGTGGGCGATAAATTAGAAGATATGCAGGCAGCGGCTGCGGCGAACGTGGGAACAAAAGTGCTGGTGCGTACGGGTAAACCTATTACGCCTGAAGCAGAAAACGCGGCGGATTGGGTGTTAAATAGCCTGGCAGACCTGCCGCAAGCGATAAAAAAGCAGCAAAAACCGGCATAATGATTAAAAGATGAGCGATTGAAATAAAAATGCATTTTTCCACTTGTCTTCCTGAG  >  REL606/226039‑226301

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