Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,890,781	G→A	A373V (GCC→GTC)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,890,781	0	G	A	100.0%	22.4 / NA	10	A373V (GCC→GTC)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (3/7);  total (3/7)

GGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCG  >  NC_000913/2890593‑2891005                                                                                                                                                                                             |                                                                                                                                                                                                                                 ggACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCGCGTCTTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGATGCGATTGAATACATACGCTGCGTCAGCGGGCGCAGCAGATTAATTATCGCTTCGATATCAATCTGTTCCGGGGAGAAACTCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATga                                                                                                                                                                    <  1:159079/251‑2 (MQ=255)                                         ggAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGGAGAAACGCACCATGGCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAgtg                                                                                                                               <  1:160822/248‑1 (MQ=255)                                                                          gtcgtAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcgtc                                                                                                                                                                                   <  1:117041/163‑1 (MQ=255)                                                                                              cGGTGACGTGTACGTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGGAGAAACTCACTATGTCAACCATCGGCTTCGGCGCGGCGTAATTGTGTAACTTCGCGTTATCGTCCACCAGCGGCAGCAGTGTTTCTGGGCGGGTAAGTTTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTc                                                                       <  1:163145/251‑1 (MQ=255)                                                                                                    cGTGTACTGCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGACGGCAGCGGGCGCAGCAGATGAATTAGCGCTTCGACATCAAGCTGGGCCCGAGAGAAACGCACCATGCCAACAATCGGCGTCGGCGCGGCGTAATGCTGGACCTTCGCTTTAGCGCCCTCCAGCGGCAGCACTGTGTAACTGCGTGTTAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCATTTc                                                                  <  1:185700/250‑1 (MQ=255)                                                                                                            tCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGTAGAAACGCACCATGTCAACAAGCGGCGGCGTCGCGGCGTAATGCTGTATCTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGc                                                           <  1:69070/249‑1 (MQ=255)                                                                                                                 ctcGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATAAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGGGTTGACGGTCAGTTCGAAGTGCCACTGt                                                     >  1:109350/1‑250 (MQ=255)                                                                                                                                cGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCagag                                     <  1:64534/251‑1 (MQ=255)                                                                                                                                                    gACGCGGGGTCAGCGGGCGAAGCAGATTAATTAGAGCTTCGACATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACg                 >  1:173235/1‑251 (MQ=255)                                                                                                                                                                   ngcgcAGCAGATTAATTAGCGCTTCGACATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACACTCGGCGTCGACGCGGCGTACTGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCg  >  1:26997/2‑251 (MQ=255)                                                                                                                                                                                             |                                                                                                                                                                                                                                 GGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCG  >  NC_000913/2890593‑2891005

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