Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,891,771	A→G	V43A (GTA→GCA)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,891,771	0	A	G	100.0%	81.9 / NA	24	V43A (GTA→GCA)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/12);  total (12/12)

ATCACGTAATGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATGCG  >  NC_000913/2891647‑2891901                                                                                                                             |                                                                                                                                   atcaCGTAATGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGc                                                                                                                     >  2:186513/1‑140 (MQ=255)   caCGTAATGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCa                                                                                                                   >  1:173748/1‑140 (MQ=255)    aCGTAATGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAg                                                                                                                  <  2:281655/140‑1 (MQ=255)          tGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAAccc                                                                                                            >  1:517456/1‑140 (MQ=255)            cTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAg                                                                                                          >  1:534644/1‑140 (MQ=255)                       cGGCGCGCTTGGCCGGTTTGGGGGGGGGAGAAAATAGTTATCCCCGGCTTTTCTGGGGCTGGCGCGGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCCCCCCAAAAATAGCCAAAAACCCAGGCAAGCggggt                                                                                               <  2:517456/140‑1 (MQ=255)                       cGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCtgggt                                                                                               <  2:373693/140‑1 (MQ=255)                       cGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCtgggt                                                                                               <  1:462692/140‑1 (MQ=255)                       cGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCtgggt                                                                                               <  1:539312/140‑1 (MQ=255)                                           cGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGGGGGGGTTAAATCGGTCGTgg                                                                           >  2:252445/1‑140 (MQ=255)                                             aGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGgcc                                                                         >  1:361633/1‑140 (MQ=255)                                                          aGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTg                                                            <  2:361633/140‑1 (MQ=255)                                                             tataCCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCa                                                         >  1:440590/1‑140 (MQ=255)                                                               taCCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGt                                                       >  1:410908/1‑140 (MQ=255)                                                                   cGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAg                                                        <  1:91907/135‑1 (MQ=255)                                                                   cGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAg                                                        >  2:91907/1‑135 (MQ=255)                                                                   cGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCt                                                   <  2:415851/140‑1 (MQ=255)                                                                     gCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCGcc                                                 >  2:333272/1‑140 (MQ=255)                                                                                           cgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCgg                           <  2:410908/140‑1 (MQ=255)                                                                                                                 tgctgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCAt     <  1:160971/140‑1 (MQ=255)                                                                                                                 tgctgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCAt     >  2:32781/1‑140 (MQ=255)                                                                                                                  gctgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATg    <  2:534644/140‑1 (MQ=255)                                                                                                                    tgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGCAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATGcg  <  2:478195/140‑1 (MQ=255)                                                                                                                    tgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATGcg  >  2:174408/1‑140 (MQ=255)                                                                                                                    tgATTGAGTGCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATGcg  <  1:333272/140‑1 (MQ=255)                                                                                                                             |                                                                                                                                   ATCACGTAATGCTTCAGCAACCCGGCGCGCATTGCCGGTTTGCGAGGCGGAGATAATAGTTATACCCGGCATTTCTGCGGCTGGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAAGGTGGGACCTGTGTCGTCATGCG  >  NC_000913/2891647‑2891901

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