Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,827,934	G→A	A109V (GCC→GTC)	astB ←	succinylarginine dihydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,827,934	0	G	A	100.0%	89.6 / NA	27	A109V (GCC→GTC)	astB	succinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/11);  total (16/11)

AGCGATTCAGTGACGGGCGCTTCCA‑‑GCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGA  >  NC_000913/1827817‑1828059                                                                                                                        |                                                                                                                              aGCGATTCAGTGACGGGCGCTTCCA‑‑GCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCgg                                                                                                          >  1:1083935/1‑139 (MQ=255)   cGATTCAGTGACGGGCGCTTCCAGCGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCg                                                                                                                >  2:960574/1‑133 (MQ=255)   cGATTCAGTGACGGGCGCTTCCAGCGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCg                                                                                                                <  1:960574/133‑1 (MQ=255)                ggCGCTTCCA‑‑GCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCAGATTgg                                                                                                                  >  1:673264/1‑116 (MQ=255)                ggCGCTTCCA‑‑GCGAACGGCGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTgg                                                                                                                  <  2:673264/116‑1 (MQ=255)                ggCGCTTCCA‑‑GCGAACGGCGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGaaa                                                                                           >  1:19018/1‑139 (MQ=255)                                 cGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCt                                                                            >  2:253574/1‑139 (MQ=255)                                   gTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACTCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTgt                                                                          >  2:857080/1‑138 (MQ=255)                                   gTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGAc                                                                                                   <  1:308912/114‑1 (MQ=255)                                   gTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGAc                                                                                                   >  2:308912/1‑114 (MQ=255)                                         tttATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCa                                                                    <  2:319309/139‑1 (MQ=255)                                              tctcaccgTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACttt                                                               >  2:846397/8‑139 (MQ=255)                                                 tCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTc                                                            <  1:568110/139‑1 (MQ=255)                                                 tCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTc                                                            <  2:463030/139‑1 (MQ=255)                                                  cAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGACGTGCAACTTTTTcc                                                           >  1:368618/1‑139 (MQ=255)                                                          cAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTg                                                                      >  2:24189/1‑120 (MQ=255)                                                          cAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTg                                                                      <  1:24189/120‑1 (MQ=255)                                                                  agaTGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGc                                           <  1:253574/139‑1 (MQ=255)                                                                        aCTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCt                                                                            <  2:191855/100‑1 (MQ=255)                                                                        aCTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCt                                                                            >  1:191855/1‑100 (MQ=255)                                                                                 tCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCa                            >  2:354336/1‑139 (MQ=255)                                                                                  ccAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCaa                           >  1:841614/1‑139 (MQ=255)                                                                                      cGTATCGGCAGATGGCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGa                       >  2:367566/1‑139 (MQ=255)                                                                                                   ggCGCGATCGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGaa          >  2:404169/1‑139 (MQ=255)                                                                                                           cGTTGCCGCATTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGTAATAAACGGa  <  1:857080/139‑1 (MQ=255)                                                                                                                    aTTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACt                                                                 <  2:805893/67‑1 (MQ=255)                                                                                                                    aTTGACTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACt                                                                 >  1:805893/1‑67 (MQ=255)                                                                                                                        |                                                                                                                              AGCGATTCAGTGACGGGCGCTTCCA‑‑GCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGCCAGTGCGGTGCCTGGCGTGCAACTTTTTCCAGTACCTGCTCATCGCTACCACTGAATCCCAACTGACGCAGCACCGGAATAAACGGA  >  NC_000913/1827817‑1828059

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