Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,085,319	G→A	A69A (GCC→GCT)	tsuA ←	thiosulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,085,319	0	G	A	100.0%	31.1 / NA	14	A69A (GCC→GCT)	tsuA	thiosulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/6);  total (8/6)

AAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGA  >  NC_000913/2085079‑2085557                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               aaaaGTCTCAGCAATAGAGTTATGTTCAGTACTGGAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGCGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGACCGAAGATATACCCACCTATAACCGTACCAAGCCACGGGAACGCCCCAGATTCGTa                                                                                                                                                                                                                                        >  1:376979/1‑249 (MQ=255)    aGTCTCAGCATAAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATCCCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCTCCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCa                                                                                                                                                                                                                                   <  1:333170/251‑1 (MQ=255)         aaGCAATCGCGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCACTCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAgg                                                                                                                                                                                                                              >  2:164127/2‑251 (MQ=255)           gCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACGGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGcc                                                                                                                                                                                                                            <  1:264082/251‑1 (MQ=255)                      atggtcaggagtGTAGGGCTGCAAGGTATGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAAATGCCGATCAATCCTTCCCCCGCGCGATAACAGGTACCGGTAGCACATCCACCCGCCAGAACAATTCCAAGCCCGAAGATATACCCACATATAATAGAACCAAGCCACGGGAACACCCCACCTTCGTAAGTCAGGCGGCCCGCCttacta                                                                                                                                                                                                                  >  2:226001/12‑245 (MQ=255)                        tttCAGTACTGTAATGCTGCACGGTTTGATTTAAACCACTGGCATGTGGAGAAAGCATCACAGCACTCATCACACTACCAGTGAAAAGCGCGATCCAACTGACGATCAATCCCTCACCCGCACGATACCAGGTCCCGGAGGCACATCAACCCGCCAGAACACTTCCCAGCCCGAAGATATACCAACCTATAACAGTACAAAGCCACGGAAACGCCCCAACTTAGTAAGTCTGTAGCCCCGCCTGAATTAAc                                                                                                                                                                                                               >  2:64164/2‑251 (MQ=255)                                    ggggctgcAAGGTTTGATTTAAACCGCGGGCATGTGGAGAACGCATCTCCGCACTCATCTCCAAAGAAGTGAAAAGCGCGATCCAACTGCCGATAAATCCCTCACCCGCGCTATACCAGGTCCGGGTGGCACATCCACCCGCCAGAACAATGCCCAGCCCGAAGATATACCCACCTATAACAGGACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTTCGTAGGCCCGCCTGAATTAACGCAAAAAccccc                                                                                                                                                                                                   <  1:206990/248‑1 (MQ=255)                                                   gATTTAAACCACTGGCATGTGGAGAACGCATCACCGCCCTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGAGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAAAGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAAGTAACGCAAAAACCCCCACGATTTGTAcagc                                                                                                                                                                                     >  2:110908/1‑247 (MQ=255)                                                                                                         tGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGCTACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcg                                                                                                                              >  2:116567/1‑251 (MQ=255)                                                                                                                        caCTGACGATCAATCCAGCACCCGCGCGAGACCAGGTCCCGGTGGAACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAAAAGTACCAAGCCACGGGCACGACCCAGCTTCGTAAGACAGTAAGCCCGCCAGAATTACCGCAACAACCCCCACGCTTTGTACCGAGATTGCAATCAGTAAGGCGTAAAACATGAGATTAGTTTTCACGAGATACATATCACGAAAACAACCAGTc                                                                                                               >  2:225760/2‑251 (MQ=255)                                                                                                                                                 gcgATCCCAGGTCCCGGTGGCACACCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGa                                                                                      <  1:277804/251‑1 (MQ=255)                                                                                                                                                                            acccgcatgAACAATTCCCAGCCCGATGATACCCCCACCTATACCCGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTa                                                           <  2:141780/243‑1 (MQ=255)                                                                                                                                                                                                                                    ccGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCAGTGAAAACATAAATGAATTAACCAGATGAACGTTAaggagg   <  2:79029/250‑1 (MQ=255)                                                                                                                                                                                                                                      gggAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAACACCCCCACGCTTTGTAACGCAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGa  >  1:367698/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               AAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGA  >  NC_000913/2085079‑2085557

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