Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,785,519	G→A	E204K (GAA→AAA)	ECB_03523 →	TrbC‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,785,519	0	G	A	100.0%	96.5 / NA	29	E204K (GAA→AAA)	ECB_03523	TrbC‑like protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/17);  total (12/17)

CTGGCCCTGCGCTTCGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTGAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGCC  >  REL606/3785447‑3785584                                                                         |                                                                  cTGGCCCTGCGCTTCGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAAc                                                                >  1:430710/1‑76 (MQ=255)  tGGCCCTGCGCTTCGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACt                                                               <  1:606256/76‑1 (MQ=255)   ggCCCTGCGCTTCGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTg                                                              >  2:709688/1‑76 (MQ=255)               cGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGc                                                  <  1:795979/76‑1 (MQ=255)               cGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGc                                                  <  1:489167/76‑1 (MQ=255)                   cGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACa                                              >  2:245259/1‑76 (MQ=255)                     tGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACaaa                                            >  2:486113/1‑76 (MQ=255)                       aaGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATc                                          <  2:341602/76‑1 (MQ=255)                        aGATTATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATcc                                         >  1:525182/1‑76 (MQ=255)                         gatgatGTGGAGTTCATCAACTATATGACCGGCGGACTCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATccc                                        <  1:709688/76‑1 (MQ=255)                          atgatgTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCt                                       >  2:410697/1‑76 (MQ=255)                          atgatgTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCt                                       <  1:547596/76‑1 (MQ=255)                                 ggAGTTGATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGaa                                >  1:146394/1‑76 (MQ=255)                                  gAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAAt                               <  2:180746/76‑1 (MQ=255)                                    gTTCATCCACTATATGAACGGCGGACGATCCCGCAGTAAACCGCTGCATGCCGGCCACAATTCCCTGCCGGAATcc                             >  1:862020/1‑76 (MQ=255)                                       catcaACTATATGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAAc                          >  1:302692/1‑76 (MQ=255)                                             cTATATGACCAGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTc                    >  1:468737/1‑76 (MQ=255)                                              tataTGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCa                   <  2:930279/76‑1 (MQ=255)                                                 atGACCGGCGGACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATc                <  1:451151/76‑1 (MQ=255)                                                     cgggcgggcgCTCCCGCAGTAGACTGCTCGATGCCGGCGATAAATCCCTGCCGGATTCCAACTCCTTCAATCCCtt            <  1:110109/74‑1 (MQ=255)                                                      cggcggACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCttt           >  1:596002/1‑76 (MQ=255)                                                      cggcggACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCttt           <  2:304217/76‑1 (MQ=255)                                                       ggcggACGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTg          <  2:468737/76‑1 (MQ=255)                                                            aCGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTAc     <  1:287908/76‑1 (MQ=255)                                                             cGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACg    >  1:263538/1‑76 (MQ=255)                                                             cGCTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACg    <  1:614621/76‑1 (MQ=255)                                                               ctaccGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGcc  <  1:486113/73‑1 (MQ=255)                                                               cccccGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGcc  <  1:245259/74‑1 (MQ=255)                                                               cTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGcc  <  1:281163/76‑1 (MQ=255)                                                               cTCCCGCAGTAAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGcc  >  2:963520/1‑76 (MQ=255)                                                                         |                                                                  CTGGCCCTGCGCTTCGGTCGTGAAGATGATGTGGAGTTCATCAACTATATGACCGGCGGACGCTCCCGCAGTGAACTGCTCCATGCCGGCGACAAATCCCTGCCGGAATCCAACTCCTTCAATCCCTTTGTCTACGCC  >  REL606/3785447‑3785584

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