Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,475,341	T→C	V126V (GTA→GTG)	tusD ←	sulfurtransferase complex subunit TusD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,475,341	0	T	C	100.0%	40.7 / NA	13	V126V (GTA→GTG)	tusD	sulfurtransferase complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/6);  total (7/6)

GACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCG  >  NC_000913/3475222‑3475469                                                                                                                        |                                                                                                                                 gaCAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTcac                                                                                                                     >  1:435652/1‑133 (MQ=255) gaCAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTcac                                                                                                                     >  2:287380/1‑133 (MQ=255) gaCAGCCAGATCGTCAGTGAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCATCTTATCAGAACAGCACGACCCTGTcac                                                                                                                     >  2:247311/1‑133 (MQ=255)      ccAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTc                                                                                                                >  2:4430/1‑133 (MQ=255)         gATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAg                                                                                                                   <  1:350946/127‑1 (MQ=255)         gATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAg                                                                                                                   >  2:350946/1‑127 (MQ=255)                                gCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACGGGTCAGCGAGGATTACGCCAGAGCCCCAAg                                                                                      >  1:246231/1‑133 (MQ=255)                                       aaGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCtt                                                                               <  1:4430/133‑1 (MQ=255)                                                       ccGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCctgct                                                               >  2:283100/1‑133 (MQ=255)                                                                 gCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGa                                                     <  1:283100/133‑1 (MQ=255)                                                                             ggTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGccc                                         <  2:498055/133‑1 (MQ=255)                                                                                 cggTAGAAAATACAAGCGCAATTCGTCTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGt                                     <  2:246231/131‑1 (MQ=255)                                                                                                acGCAATTCGTTTCATCAGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGtt                      >  2:366189/1‑133 (MQ=255)                                                                                                             catcaGAACTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGc         <  1:305560/133‑1 (MQ=255)                                                                                                                    aCTGCACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCg  <  1:287380/133‑1 (MQ=255)                                                                                                                        |                                                                                                                                 GACAGCCAGATCGTCAGTTAATGCGGAAGTTGCCAGTAAAGCATCTAAACCTTCCCGGCCTGCGGCTGTACCATGAGGTGCAGTAGAAAAAACAAACGCAATTCGTTTCATCAGAACTGTACCACCCTGTCACAGGTCAGCGAGGCTTCCGCCAGCGCCCCAAGTCCGCTTAAGGTAAATCCCTGCTGAAGGTTTGACGAAGCCAGCCCCAGTCTTCCGGCCTCCGTTTCATCAACAACGCCACGGCG  >  NC_000913/3475222‑3475469

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