Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,489,824	T→A	A119A (GCA→GCT)	argD ←	N‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,489,824	0	T	A	100.0%	141.0 / NA	40	A119A (GCA→GCT)	argD	N‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (23/17);  total (23/17)

CGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGTTGCCT  >  NC_000913/3489746‑3489897                                                                               |                                                                          cgagcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGc                                                                         <  2:666415/81‑1 (MQ=255)  gagcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGc                                                                         >  2:666458/1‑80 (MQ=255)     cgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGt                                                                     <  2:666416/81‑1 (MQ=255)     cgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGt                                                                     >  1:666472/1‑81 (MQ=255)         cATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTa                                                                  >  2:666479/1‑80 (MQ=255)               aaGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCgg                                                           >  1:666475/1‑81 (MQ=255)               aaGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCgg                                                           >  2:666453/1‑81 (MQ=255)                aGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGt                                                          >  1:666466/1‑81 (MQ=255)                aGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGt                                                          <  1:666417/81‑1 (MQ=255)                  cGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGttt                                                        <  2:666418/81‑1 (MQ=255)                   gTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTc                                                       >  1:666467/1‑81 (MQ=255)                   gTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTc                                                       >  2:666494/1‑81 (MQ=255)                     tATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGt                                                     >  1:666460/1‑81 (MQ=255)                      aTGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGt                                                     <  2:666419/80‑1 (MQ=255)                       tGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGtt                                                    <  2:666420/80‑1 (MQ=255)                          aaGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCt                                                >  2:666457/1‑81 (MQ=255)                            ggCAATAATCTTTGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTc                                              <  2:666421/81‑1 (MQ=255)                                       tGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGTGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTTCCGGAGt                                   >  1:666482/1‑81 (MQ=255)                                         gTTTTGAACGGGCTATGACGCACACAGGCGTGATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTc                                 <  2:666422/81‑1 (MQ=255)                                         gTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGtt                                  <  1:666423/80‑1 (MQ=255)                                         gTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGtt                                  >  1:666480/1‑80 (MQ=255)                                            ttGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCAt                               <  1:666424/80‑1 (MQ=255)                                             tGAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATa                              >  2:666477/1‑80 (MQ=255)                                              gAACGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATaaa                            >  1:666487/1‑81 (MQ=255)                                                 cGGGCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAacac                         >  1:666471/1‑81 (MQ=255)                                                  gggCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAacac                         >  1:666463/1‑80 (MQ=255)                                                  gggCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAAcacc                        <  2:666425/81‑1 (MQ=255)                                                   ggCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAAcacca                       >  2:666492/1‑81 (MQ=255)                                                    gCTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAAcaccac                      <  1:666426/81‑1 (MQ=255)                                                     cTATGACGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCAcg                     >  1:666464/1‑81 (MQ=255)                                                          aCGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGtt                 <  2:666427/80‑1 (MQ=255)                                                          aCGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTc                <  1:666428/81‑1 (MQ=255)                                                           cGCACACAGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCg               >  2:666468/1‑81 (MQ=255)                                                                acaGGCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCaaa          >  2:666491/1‑81 (MQ=255)                                                                    gCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAATCACCACGCGTTCGGCAAACGtt      <  1:666429/81‑1 (MQ=255)                                                                    gCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGtt      >  1:666470/1‑81 (MQ=255)                                                                    gCGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGtt      <  1:666430/81‑1 (MQ=255)                                                                     cGTAATGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGTTg     >  1:666476/1‑81 (MQ=255)                                                                         aTGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGTTGCCt  <  2:666431/80‑1 (MQ=255)                                                                         aTGGCGAGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGTTGCCt  >  1:666489/1‑80 (MQ=255)                                                                               |                                                                          CGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTTCCGTGCCGGAGTTCATAAACACCACGCGTTCGGCAAACGTTGCCT  >  NC_000913/3489746‑3489897

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