Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	REL606	1,930,575	C→A	100%	F137L (TTC→TTA)	yecE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,930,575	0	C	A	100.0%	93.1 / NA	29	F137L (TTC→TTA)	yecE	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/18);  total (11/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGT  >  REL606/1930428‑1930718                                                                                                                                                    |                                                                                                                                                gTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGc                                                                                                                                               <  2:49134/150‑1 (MQ=255)              ccGCGCATTGGGCAATACTGGCTGCAACTGCCTTCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                                                                                                                  <  1:554064/150‑1 (MQ=255)              ccGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                                                                                                                  <  2:353352/150‑1 (MQ=255)               cgcgCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAgg                                                                                                                                 <  2:242909/150‑1 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCTACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  1:365997/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  2:327632/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  1:643012/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  1:313817/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  2:130811/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  1:80861/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  1:64214/1‑150 (MQ=255)                       gggCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                                                                         >  2:327181/1‑150 (MQ=255)                                 ggCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcgcg                                                                                                               >  2:554082/1‑150 (MQ=255)                                          tgcctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc                                                                                                      <  2:64405/150‑1 (MQ=255)                                          tgcctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc                                                                                                      <  2:367058/150‑1 (MQ=255)                                          tgcctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc                                                                                                      <  1:377057/150‑1 (MQ=255)                                          tgcctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc                                                                                                      <  1:498757/150‑1 (MQ=255)                                                         gctgaCGGGCGCCGCCTGCGCTTTGGCATTTTCTCTATTCTCTTCCTGGCCAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATc                                                                                       <  1:362598/146‑1 (MQ=255)                                                                       cctgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGAc                                                                         <  1:326705/150‑1 (MQ=255)                                                                       cctgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGAc                                                                         <  2:314707/150‑1 (MQ=255)                                                                                      tttCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCAt                                                          >  1:670151/1‑150 (MQ=255)                                                                                                     ccTTGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCcacac                                           <  1:84764/150‑1 (MQ=255)                                                                                                     ccTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCcacac                                           <  1:130434/150‑1 (MQ=255)                                                                                                       tgtcGAATTTAATTATGGCGTGGAAGTCGGCCCTCCACAGTTTTTGGCCAACGGAGAAGAGGAACAAACGCTTGATCGCGGTTCACAGCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCGGGTTCATGCAGCACGTCCACACAg                                         <  1:287715/147‑1 (MQ=255)                                                                                                       tGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAg                                         <  2:585775/150‑1 (MQ=255)                                                                                                       tGGCGAATTTAATTATGGCGTGGAAGTCCGCCAGCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAACCGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAg                                         <  1:555385/150‑1 (MQ=255)                                                                                                             ctttAATTATGGCGGGGAAGTGCGCCATCCACAGTTTTTAGCCAAAGGAGAATAGGAACAAACGCTTAATCGCGGGTTACATCAGCGCGGCGTTAATCGGGGGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGc                                   <  1:383223/149‑1 (MQ=255)                                                                                                             aTTTGATTATGGCGCGCAAGTCCGCCATCCACAGTTTTTATCCAAAGGAGAAGAGGAAGAAACGCTTAGTCGAGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGc                                   <  1:386416/150‑1 (MQ=255)                                                                                                             aTTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGc                                   <  2:81095/150‑1 (MQ=255)                                                                                                               ttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCATGTCCACACAGTGAAGCTa                                 <  1:552483/150‑1 (MQ=255)                                                                                                                                              gTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGAAAGCCGCCCGGTTCATGCATCACGTCCACCCAGTGAAGCTATTCGCGACGCTCAACGAAAAACACCTTAAGt  >  1:665263/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                GTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGT  >  REL606/1930428‑1930718

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