Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,191,861	A→G	R42R (CGT→CGC)	purB ←	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,191,861	0	A	G	100.0%	63.3 / NA	19	R42R (CGT→CGC)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/10);  total (9/10)

TTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGATAATTCCATA  >  NC_000913/1191737‑1191987                                                                                                                             |                                                                                                                               ttttGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTc                                                                                                                        <  2:515855/133‑1 (MQ=255) ttttGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTc                                                                                                                        <  2:31806/133‑1 (MQ=255)      aTGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACtt                                                                                                                   >  1:537506/1‑133 (MQ=255)           cgccgcATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTaca                                                                                                              <  2:537506/133‑1 (MQ=255)                cATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGaa                                                                                                         >  1:131237/1‑133 (MQ=255)                        tcGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCa                                                                                                 >  2:309168/1‑133 (MQ=255)                                                             ttGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGAc                                                            >  1:448358/1‑133 (MQ=255)                                                                                     cTTCCTTGATCACTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGAc                                    <  2:7283/133‑1 (MQ=255)                                                                                        ccTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAgg                                 >  1:329270/1‑133 (MQ=255)                                                                                               cGCTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTATATTGCCGTAGCGTCCGTCGAGAGGGGAAAcg                          >  2:165195/1‑133 (MQ=255)                                                                                                 cTGCGTGCGCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAAcggc                        >  2:436542/1‑133 (MQ=255)                                                                                                  tgcgtgcgTGGCCAGATTATGCCGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAAcggcg                       <  1:165195/133‑1 (MQ=255)                                                                                                  tgcgtgcgCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAAcggcg                       >  2:55105/1‑133 (MQ=255)                                                                                                  tgcgtgcgCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGATGACTTTAACGCCGTACCGTCCATCTACAGAGGAAAcggcg                       >  2:199107/1‑133 (MQ=255)                                                                                                    cgtgcgCGGCCAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGt                     <  2:106331/133‑1 (MQ=255)                                                                                                             ccAGTTTTTGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGa            <  2:329270/133‑1 (MQ=255)                                                                                                             ccAGTTTTCGCAGCCGGCTTACCTCAACTTCTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGa            <  1:199107/133‑1 (MQ=255)                                                                                                                     tGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGATAATTCCa    <  2:102405/133‑1 (MQ=255)                                                                                                                     tGCAGCCAGCGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGATAATTCCa    <  2:448358/133‑1 (MQ=255)                                                                                                                       cagccagcGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAGACGGCGGTCAGTGAGGATAATTCCATa  <  2:71189/133‑1 (MQ=255)                                                                                                                             |                                                                                                                               TTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTAACCGATTGCGTCGGCAGCAAAAGCAGGAACTTCCTTGATCGCTGCGTGCGCGGCCAGTTTTTGCAGCCAACGTACTTCAACTTGTACACGGAATTTCAGCAAACCATATTCGCTGAAAATCCCGCGCAGCGCGCTGACTTTATCGCCGTAGCGTCCATCGACAGGGGAAACGGCGGTCAGTGAGGATAATTCCATA  >  NC_000913/1191737‑1191987

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