Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,887,602	C→T	R237R (AGG→AGA)	cysH ←	phosphoadenosine phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,887,602	0	C	T	100.0%	45.2 / NA	17	R237R (AGG→AGA)	cysH	phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/10);  total (7/10)

GCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAT  >  NC_000913/2887367‑2887816                                                                                                                                                                                                                                            |                                                                                                                                                                                                                       gCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGATTCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAg                                                                                                                                                                                                         <  1:134423/251‑1 (MQ=255) gCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGGTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTGTAATTTGGCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAg                                                                                                                                                                                                         <  1:292218/251‑1 (MQ=255)    ggCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAGCTGGAACCTCGCGAATGCTTAACGCCAGACATCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAAGGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTGTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGc                                                                                                                                                                                                              <  1:157749/243‑1 (MQ=255)                        atCGCCTGGGGCTTTTGTCCATCTGGAACATCGCGAATCCTTAACGCCAGACAGCATGAAGAACCCGACGCCATTCTATACCAGATAATAATTCCAGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAACTAAAAACGCCACCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAAAGCGTTTCTTCTTCTGCCAt                                                                                                                                                                                  >  1:293986/1‑251 (MQ=255)                                           cATCTGGAACCGCGCGAATGCTTAACGCCAGACAGCCTCAAGCCCCCGACGCCATTCTAGACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAAGCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCAGCATTTCTGACAGAGGCGTTTAATGTGTCCGGCAATATTTAGCCGTCGTGGAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGcc                                                                                                                                                                                    <  1:164359/230‑1 (MQ=255)                                            aTCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATGTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGaaan                                                                                                                                                                                                     <  1:199496/212‑2 (MQ=255)                                                    ncTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGACATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAACAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACggggaatcgg                                                                                                                                                       >  1:66412/2‑243 (MQ=255)                                                                               cTCAAGCACCCGACGCCATTCTATCCCTGATAATTCTTCACACGGCGTTCACGTTTCAAATAAAGAATGCTACTCGCAAAAAAAACACGCCCCAATCATTTCTGACAGAGGCGATTAATTTGTCAGGCAATCATTACCCTTCGCGAAACACACATTCTCTTTTTAAGCCAAAGAAACGCGTTTATTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGAATCGCACACCGATAAATATACTTc                                                                                                                            >  1:218644/1‑250 (MQ=255)                                                                                            cGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATgggc                                                                                                              <  1:269791/251‑2 (MQ=255)                                                                                               naTTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCCAAATAAAAACGCCACCATCATTTATGACAGAGGCGTTTCATTTGTACGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGat                                                                                                           >  1:37754/2‑251 (MQ=255)                                                                                                      taCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAAGAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAAGTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTCTCTGCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCGTCATCCCATAATgg                                                                                                                <  1:115913/239‑1 (MQ=255)                                                                                                                                                         cGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGc                                                                                                                                                                                                              <  1:145921/94‑1 (MQ=255)                                                                                                                                                             aaaTAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                             <  1:26423/251‑1 (MQ=255)                                                                                                                                                                 aaaaaCGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagca                                         <  1:151848/251‑1 (MQ=255)                                                                                                                                                                     ccGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACCTTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATc                                                             >  1:87430/2‑227 (MQ=255)                                                                                                                                                                                                        ntCCGGCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTATTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAt  >  1:61008/2‑251 (MQ=255)                                                                                                                                                                                                            ggCAATATTTACCCTTCGTGTAACCCACATTCTCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATANTggnnn                                                                                                             >  2:205455/1‑137 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                       GCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAT  >  NC_000913/2887367‑2887816

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