Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,156,962	T→C	V38A (GTG→GCG)	argH →	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,156,962	0	T	C	100.0%	62.6 / NA	19	V38A (GTG→GCG)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/7);  total (12/7)

AACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGA  >  NC_000913/4156841‑4157088                                                                                                                          |                                                                                                                               aaCAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCa                                                                                                                     <  1:473650/133‑1 (MQ=255) aaCAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCa                                                                                                                     <  2:486987/133‑1 (MQ=255)                cTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACgg                                                                                                      <  1:410092/133‑1 (MQ=255)                  ttGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTa                                                                                                    >  1:325386/1‑133 (MQ=255)                  ttGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTa                                                                                                    >  2:421545/1‑133 (MQ=255)                          ggcgTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCAGGTCACGGTCGGCGTGtt                                                                                            >  1:240558/1‑133 (MQ=255)                                             agATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGgcg                                                                         >  2:274618/1‑133 (MQ=255)                                                  aaCGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACt                                                                    >  1:85441/1‑133 (MQ=255)                                                              aaTTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCt                                                        <  2:445347/133‑1 (MQ=255)                                                                 tCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGaa                                                     <  2:93702/133‑1 (MQ=255)                                                                           aCTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTg                                           >  1:84842/1‑133 (MQ=255)                                                                                    tGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGtt                                  <  1:274618/133‑1 (MQ=255)                                                                                    tGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGtt                                  <  2:85441/133‑1 (MQ=255)                                                                                         aCCGTCTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCACCTGGAAGACGCGATGATCGTGTTGCGGGAAGATGTTCgcgc                             >  2:221158/1‑133 (MQ=255)                                                                                              cTGGCGGAGCAGGATATTGTTGGCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGc                        >  1:259275/1‑133 (MQ=255)                                                                                                       cAGGATATTGTTGGCTCTACGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACaaa               >  2:265480/1‑133 (MQ=255)                                                                                                                    gCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGa  >  2:410702/1‑133 (MQ=255)                                                                                                                    gCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGa  >  2:118385/1‑133 (MQ=255)                                                                                                                    gCTCTGCGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGa  >  1:394485/1‑133 (MQ=255)                                                                                                                          |                                                                                                                               AACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGA  >  NC_000913/4156841‑4157088

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