Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,797	G→T	S172* (TCG→TAG)	argI ←	ornithine carbamoyltransferase ArgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,797	0	G	T	100.0%	91.4 / NA	27	S172* (TCG→TAG)	argI	ornithine carbamoyltransferase ArgI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (11/16);  total (11/16)

AGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGC  >  NC_000913/4477664‑4477929                                                                                                                                      |                                                                                                                                     agtcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGc                                                                                                                                <  1:60175/140‑1 (MQ=255)   tcagCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGccc                                                                                                                              <  2:80580/140‑1 (MQ=255)                      aTTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCAc                                                                                                           <  2:113705/140‑1 (MQ=255)                          tGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGc                                                                                                       >  2:217128/1‑140 (MQ=255)                          tGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTCGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGc                                                                                                       >  2:328499/1‑140 (MQ=255)                               tgCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCTGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                                                                  >  1:251859/1‑140 (MQ=255)                               tgCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                                                                  >  2:251559/1‑140 (MQ=255)                                 ccAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGAcc                                                                                                >  2:269251/1‑140 (MQ=255)                                                ccGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                                                                 <  2:64321/140‑1 (MQ=255)                                                          cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCGCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                                       <  1:169631/140‑1 (MQ=255)                                                          cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                                       >  1:104521/1‑140 (MQ=255)                                                          cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                                       >  1:237270/1‑140 (MQ=255)                                                          cgcAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                                       >  2:292161/1‑140 (MQ=255)                                                                ttCCGGCCAGCACGCTCGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                                 <  1:328499/140‑1 (MQ=255)                                                                         gCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCt                                                        >  1:75744/1‑140 (MQ=255)                                                                          cACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                                       >  1:34991/1‑140 (MQ=255)                                                                                         aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                     >  2:25368/1‑127 (MQ=255)                                                                                         aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                     <  1:25368/127‑1 (MQ=255)                                                                                                 cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagca                                <  2:242384/140‑1 (MQ=255)                                                                                                    aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                             <  2:37994/140‑1 (MQ=255)                                                                                                       cAGACCGGTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                          <  1:34410/140‑1 (MQ=255)                                                                                                              gTAAGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                   <  1:292161/140‑1 (MQ=255)                                                                                                                aaGCGCCGCAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                 <  1:217128/140‑1 (MQ=255)                                                                                                                     ccgcAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa            <  1:92179/140‑1 (MQ=255)                                                                                                                     ccgcAGCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa            <  1:184377/140‑1 (MQ=255)                                                                                                                          gCTTCGAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt       <  1:251559/140‑1 (MQ=255)                                                                                                                               gAGCATCTAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc  <  1:269251/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     AGTCAGCGTAATATTCCCACCATTTTGCTGTGCCAGGGCGCGGCATTCCGTAACCAGCGCAGCTTCCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGC  >  NC_000913/4477664‑4477929

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: ‑