Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,490,419	A→T	D2E (GAT→GAA)	glpE ←	thiosulfate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,490,419	0	A	T	100.0%	94.0 / NA	25	D2E (GAT→GAA)	glpE	thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (15/10);  total (15/10)

GAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAA  >  REL606/3490346‑3490491                                                                          |                                                                         gAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTc                                                                        <  1:247081/76‑1 (MQ=255)  aaTATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTcc                                                                       >  2:118405/1‑76 (MQ=255)         aCCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGtc                                                                >  2:463010/1‑76 (MQ=255)               aCCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCtctt                                                          >  1:892210/1‑76 (MQ=255)                    ctctTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTctttcttt                                                     <  1:1032989/76‑1 (MQ=255)                     tctTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTa                                                    >  1:696984/1‑76 (MQ=255)                     tctTGTTCCGGCAACTTCTGGTAGGCGTCCGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCtctttcttca                                                    <  1:1002971/76‑3 (MQ=255)                       tttttCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTaca                                                  <  1:560115/76‑1 (MQ=255)                          ttCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTacaaac                                               >  2:92098/1‑76 (MQ=255)                            ccTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTacaaacaa                                             <  2:520049/76‑1 (MQ=255)                                cAACTTCTNGCGCGCGTCGNCAACGTTAATACATTCGAACTGNTCCATGCGTCTCTCTTTCTTTACAAACAAGTgg                                         <  2:970242/76‑1 (MQ=255)                                    ttCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCaa                                     >  2:656886/1‑76 (MQ=255)                                       tGGTGCGCGTCGNCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAAtt                                  >  1:16549/1‑76 (MQ=255)                                              cgTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGca                           >  2:265994/1‑76 (MQ=255)                                               gTCGGCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGcac                          >  1:643558/1‑76 (MQ=255)                                                  tgCAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGt                       >  2:706141/2‑76 (MQ=255)                                                    cAACGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGttt                     >  2:597735/1‑76 (MQ=255)                                                     aaCGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTa                    <  1:990746/76‑1 (MQ=255)                                                      aCGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTAc                   <  1:119445/76‑1 (MQ=255)                                                       cGTTAATACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACg                  >  2:701830/1‑76 (MQ=255)                                                           aaTACATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGa              <  2:643558/76‑1 (MQ=255)                                                              aCATTCGAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGc           >  1:922473/1‑76 (MQ=255)                                                                    gAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGa     >  1:706609/1‑76 (MQ=255)                                                                    gAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGa     >  2:325174/1‑76 (MQ=255)                                                                    gAACTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGa     <  1:183912/76‑1 (MQ=255)                                                                       cTGTTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATaa  <  1:463010/76‑1 (MQ=255)                                                                          |                                                                         GAATATCGACCAGCACCGCCTCTTTTTCCTGCAACTTCTGGTGCGCGTCGGCAACGTTAATACATTCGAACTGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAA  >  REL606/3490346‑3490491

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