Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	416,124	A→G	D207G (GAT→GGT)	ECB_00376 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	416,124	0	A	G	100.0%	93.0 / NA	24	D207G (GAT→GGT)	ECB_00376	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/10);  total (14/10)

CGTTACCCGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGATAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCATTCGTCGCCC  >  REL606/416052‑416196                                                                         |                                                                         cGTTACCCGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAg                                                                       >  2:647320/1‑76 (MQ=255) cGTTACCCGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAg                                                                       >  2:126691/1‑76 (MQ=255)       ccGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCtt                                                                 <  2:578710/76‑1 (MQ=255)        cGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTc                                                                <  2:380055/76‑1 (MQ=255)         gTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCg                                                               >  1:645414/1‑76 (MQ=255)          tGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGc                                                              >  1:614476/1‑76 (MQ=255)          tGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGc                                                              >  1:614487/1‑76 (MQ=255)             ttACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGccac                                                           >  2:173726/1‑76 (MQ=255)              tACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGccacc                                                          >  2:235521/1‑76 (MQ=255)                  gTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGccaccacca                                                      <  1:467972/76‑1 (MQ=255)                     aCCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATa                                                   <  1:74843/76‑1 (MQ=255)                      cccGTTCTGGGATGCGATATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATaa                                                  <  1:84428/76‑1 (MQ=255)                        cGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAAt                                                >  2:525106/1‑76 (MQ=255)                                   gCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTg                                     >  2:550067/1‑76 (MQ=255)                                      agaTAATTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGcc                                  >  2:88876/1‑76 (MQ=255)                                          aaTTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATc                              >  1:446882/1‑76 (MQ=255)                                          aaTTTTTTAGATAAAGGTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATc                              <  1:151652/76‑1 (MQ=255)                                                         ggTTTTTCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCa               <  2:650067/76‑1 (MQ=255)                                                             tttCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCAt           >  2:546072/1‑76 (MQ=255)                                                              ttCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCAtt          <  1:499833/76‑1 (MQ=255)                                                              ttCCTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCAtt          <  2:555951/76‑1 (MQ=255)                                                                 cTGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCATtcgt       >  1:203915/1‑76 (MQ=255)                                                                  tGCATGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCATtcgtc      >  2:206200/1‑76 (MQ=255)                                                                      tGGGTAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCATTCGTCGccc  <  1:593299/76‑1 (MQ=255)                                                                         |                                                                         CGTTACCCGTGATTACGGTTACCCGTTCTGGGATGCGAGATAATTTTTTAGATAAAGGTTTTTCCTGCATGGATAGCCTCTTCGCCACCACCAATAAATTAGGCGAAGTGGCCAATCTTCGCGTGACATTCAGCATTCGTCGCCC  >  REL606/416052‑416196

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