Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,392,351	G→T	A52S (GCA→TCA)	yfcN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,392,351	0	G	T	100.0%	112.0 / NA	34	A52S (GCA→TCA)	yfcN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/19);  total (15/19)

GGACACGATTGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGGCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGCCCGGTGAAA  >  REL606/2392278‑2392416                                                                          |                                                                  ggACACGATTGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCa                                                                 >  1:182906/1‑76 (MQ=255)        aTTGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgcca                                                          <  1:1097252/76‑1 (MQ=255)         ttGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgccag                                                         <  2:812725/76‑1 (MQ=255)          tGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgccagc                                                        >  1:1068845/1‑76 (MQ=255)          tGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgccagc                                                        >  1:407941/1‑76 (MQ=255)           gTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgccagcc                                                       >  2:1001636/1‑76 (MQ=255)           gTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATgccagcc                                                       >  1:1359244/1‑76 (MQ=255)              cACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCatt                                                    <  2:1396868/76‑1 (MQ=255)              cACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCatt                                                    >  2:1344217/1‑76 (MQ=255)                ccgaccgCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCattat                                                  <  1:1045316/76‑1 (MQ=255)                 cgaccgCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCattatt                                                 >  2:467971/1‑76 (MQ=255)                  gaccgCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTAttt                                                >  2:1530290/1‑76 (MQ=255)                  gaccgCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTAttt                                                >  1:1196894/1‑76 (MQ=255)                     cgCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTtctc                                             <  2:606268/76‑1 (MQ=255)                        aGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGa                                          >  2:1625201/1‑76 (MQ=255)                           gTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGa                                       <  2:1587246/76‑1 (MQ=255)                              aGAAAATCAGCGAAGTGCCGGTTAAGGGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGtt                                    <  2:1199289/76‑1 (MQ=255)                                aaaaTCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTc                                  <  2:733872/76‑1 (MQ=255)                                aaaaTCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTc                                  <  2:1357743/76‑1 (MQ=255)                                 aaaTCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCa                                 <  1:341051/76‑1 (MQ=255)                                        cGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTa                          >  1:737285/1‑76 (MQ=255)                                         gAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTAt                         >  1:380919/1‑76 (MQ=255)                                         gAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTAt                         <  1:952237/76‑1 (MQ=255)                                          aaGTGCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTAtt                        <  1:1305098/76‑1 (MQ=255)                                              gCCGGTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAAt                    <  1:625025/76‑1 (MQ=255)                                                  gTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCg                <  1:1202857/76‑1 (MQ=255)                                                  gTTAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCg                <  2:684752/76‑1 (MQ=255)                                                   ttAAGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGa               <  2:331701/76‑1 (MQ=255)                                                     aaGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAg             <  2:1485980/76‑1 (MQ=255)                                                     aaGCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAg             <  2:1159096/76‑1 (MQ=255)                                                       gCGGTTGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGc           >  1:1713881/1‑76 (MQ=255)                                                           ttGATTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGCCCgg       >  2:1077815/1‑76 (MQ=255)                                                              aTTCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGCCCGGTGa    <  2:796013/76‑1 (MQ=255)                                                                tCAGGAACAGTCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGCCCGGTGaaa  >  1:1018541/1‑76 (MQ=255)                                                                          |                                                                  GGACACGATTGTCCACCGACCGCAGCGTAAGAAAATCAGCGAAGTGCCGGTTAAGCGGTTGATTCAGGAACAGGCAGATGCCAGCCATTATTTCTCAGATGAGTTTCAGCCGTTATTGAATACCGAAGGCCCGGTGAAA  >  REL606/2392278‑2392416

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