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%	97.7 / NA	30	A69A (GCC→GCT)	tsuA	thiosulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/13);  total (17/13)

AAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTA  >  NC_000913/2085081‑2085531                                                                                                                                                                                                                                                |                                                                                                                                                                                                                     aaGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACAATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGAAACGCCCCAGCTTCGTAAGTc                                                                                                                                                                                                          >  1:542543/1‑251 (MQ=255)                  ttATGTTAAGTACTGTAGTGCGGCAAGGTTTGATTTACACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCAGATCAATACCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGACAAATTCCCAGCCCGACGATATACCAACCTATAACAGTACCAAGCCACGGGAAAGCCCCAGCTTCGTCAGTCAGTAGGCCCGACTGAAt                                                                                                                                                                                         >  2:359206/1‑251 (MQ=255)                             aCTGTAGTGCTGCAAGGTTGGAGTTAAACCACTGGAATGGGGAGAACCCACCACCGCACCCATCACCATAAAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATCCCAGGTCCCGTTGGCACATCC‑CCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCaaaaa                                                                                                                                                                              <  1:293538/251‑1 (MQ=255)                                  aGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGCTACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCCCGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGCATTAACGCAAAAAccccc                                                                                                                                                                         >  1:444754/1‑251 (MQ=255)                                     gctgcAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGCTACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCCGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACg                                                                                                                                                                      >  1:432080/1‑251 (MQ=255)                                         cAAGGTTTGATTTAAGCCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAACGCGATCCAACTGCCGATCAAACACTCACCCGCGCGATACCAGGGCTCGGCGGCACATCA‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATCCCCACCTATAACAGTAACAACCCCCGGAAACGCACCAGCTTCGTAAGCCTGTAGGCCAGCCTGCATCAACGCAAAAAACCCCaagattt                                                                                                                                                                  >  1:588276/1‑245 (MQ=255)                                                         ccACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAAt                                                                                                                                                  >  2:314768/1‑251 (MQ=255)                                                                   gtgGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGc                                                                                                                                         <  1:314768/250‑1 (MQ=255)                                                                           cGCACCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACGGCCGATCAAGCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACAt                                                                                                                                <  1:196392/251‑1 (MQ=255)                                                                                         catccccaaataaggaaaaaGCGCCATCCAACTGCCGATCAATCCCTCGCCCGCGCGATACCAGGTCCCGGTGGCCCATCC‑ACCCGCCATAACAATTCCCAGCCCGAATATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCATCAAGGCGTAAAACATGCGATTAttttg                                                                                                                    <  1:489521/233‑2 (MQ=255)                                                                                             aCCATATAAGTGTAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGAGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGc                                                                                                                                                                                                   >  2:302976/1‑166 (MQ=255)                                                                                             aCCATATAAGTGAAAAGCGCGACCCAACTGCCGATCAATCCCTCACCCTCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCGCCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGc                                                                                                                                                                                                   <  1:302976/166‑1 (MQ=255)                                                                                                     agaaaaacgccctatacAACGGCCGACCAACCCCTCACCCGCGCGATACCAGGTCCCGGTGGCAAAGCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCTCGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATGTTTCACGATATACatat                                                                                                      <  2:302698/235‑1 (MQ=255)                                                                                                       tGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcg                                                                                                    >  1:484635/1‑251 (MQ=255)                                                                                                       gaaaacgcgcattccacctcccGATCAACCCCCCCCCCGCGCGATACCAGGTCCCGGGGGCACACCC‑CCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcg                                                                                                    <  2:420337/231‑1 (MQ=255)                                                                                                                      aaccgccgcccaaccccTCACCCCCGCGATACCAGGTCCCG‑TTGCACCTCCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCGTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                     <  2:444754/237‑1 (MQ=255)                                                                                                                                            cgcgcgATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCAGCAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGACTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAACAGCCCAAGCTGc                                                                >  2:448677/1‑250 (MQ=255)                                                                                                                                                cgATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGAc                                                           >  1:34872/1‑251 (MQ=255)                                                                                                                                                      cAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATcc                                                     >  1:99888/1‑251 (MQ=255)                                                                                                                                                                           aCCCGCCAGAACAATGCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCTCGCTTTGTACCGAAATTGCAATCAGCAAGGCTTAAAACATGCGATTATTTTTCACGATATACATATCTCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTa                                 <  1:579536/251‑1 (MQ=255)                                                                                                                                                                               gCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGccc                             <  2:484635/251‑1 (MQ=255)                                                                                                                                                                                 cAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAAAGCACCACAAATTAGCCCGc                           >  2:167656/1‑251 (MQ=255)                                                                                                                                                                                       aaTTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAat                      >  1:251844/1‑250 (MQ=255)                                                                                                                                                                                        aTTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCGTAatat                    <  1:471539/251‑1 (MQ=255)                                                                                                                                                                                          tCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATc                   >  2:130742/1‑250 (MQ=255)                                                                                                                                                                                          tCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATc                   <  2:99888/250‑1 (MQ=255)                                                                                                                                                                                             cAGCCCGAAAATATCCCCTCCTATAACAGTCCCAAGCCCCGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAAGCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTg               <  2:467951/251‑1 (MQ=255)                                                                                                                                                                                                 ccGAAGATATACACACCTATAACAGTACCAAGCCACGGGAACGCCCCAGATTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGATTTGTACCGCAATTGCAATAAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATAGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGACTGACAAATCCAAGCAGAGCAACACAAATTAGCCCGCTTAATATCATTGaaa            >  1:466899/1‑250 (MQ=255)                                                                                                                                                                                                      gATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATa       >  1:531375/1‑250 (MQ=255)                                                                                                                                                                                                           aCCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGACAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaatta  >  2:182375/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                     AAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCC‑ACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTA  >  NC_000913/2085081‑2085531

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