Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,559,164	A→G	R121R (CGT→CGC)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,164	0	A	G	100.0%	148.8 / NA	43	R121R (CGT→CGC)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/19);  total (24/19)

ATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTC  >  NC_000913/4559033‑4559298                                                                                                                                    |                                                                                                                                       aTGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCAGGTCAGCACCCAAGCGCTGATGCCTTCTTCGTTGAGCGCGCGGGTCtt                                                                                                                                >  2:342030/1‑140 (MQ=255)   ggCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTgg                                                                                                                              >  1:646767/1‑140 (MQ=255)                    tCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCggg                                                                                                             <  1:112295/140‑1 (MQ=255)                    tCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCggg                                                                                                             >  2:768879/1‑140 (MQ=255)                    tCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCggg                                                                                                             <  1:93366/140‑1 (MQ=255)                         cgatcgatAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCg                                                                                                        <  2:483870/140‑1 (MQ=255)                               atAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGTGGGTCTTTGCGAGCAGGGATTACGGGTGGCGAGAGAt                                                                                                  >  1:403770/1‑140 (MQ=255)                                    cGCCACGTTTTTTTCCTCGGAACCCGTAATGTTGCTGGAAGGGACATGGTAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGt                                                                                             <  2:403770/140‑1 (MQ=255)                                    cGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAAGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGt                                                                                             >  2:595206/1‑140 (MQ=255)                                          gTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGc                                                                                                         <  1:99377/122‑1 (MQ=255)                                          gTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGc                                                                                                         >  2:99377/1‑122 (MQ=255)                                          gTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGc                                                                                       >  2:654094/1‑140 (MQ=255)                                               tttCCACGGAACCCGTGGTGGTGCGGGAAGGGAGACGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagc                                                                                  <  2:393045/140‑1 (MQ=255)                                                ttCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagca                                                                                 <  2:646767/140‑1 (MQ=255)                                                  ccACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAACATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGa                                                                               <  1:206858/140‑1 (MQ=255)                                                  ccACGGAACCCGTAATGGTGCGGGAAAGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGa                                                                               <  1:595206/140‑1 (MQ=255)                                                     cGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGCGAAAGAGTCGGTGCCCAGCAGacca                                                                            >  1:393045/1‑140 (MQ=255)                                                      ggAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGC‑‑GAGATAGAGTCGGTGCCCAGCAGaccaacc                                                                         >  1:196373/1‑140 (MQ=255)                                                                 tGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgtg                                                                <  1:768879/140‑1 (MQ=255)                                                                  ggTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgtga                                                               >  1:767093/1‑140 (MQ=255)                                                                    tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgtgacg                                                             <  1:654094/140‑1 (MQ=255)                                                                     gCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGc                                                            >  2:853412/1‑140 (MQ=255)                                                                      cGGGAAGGGACATGATAAGCCCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCTGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGcc                                                           >  1:471813/1‑140 (MQ=255)                                                                        ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCg                                                         <  1:176013/140‑1 (MQ=255)                                                                             ggACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCAGTGACCAGCAGACCAATCACTGACGTGACGCCCGCTTcc                                                    >  2:300722/1‑140 (MQ=255)                                                                              gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCg                                                   >  1:713184/1‑140 (MQ=255)                                                                                       aGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGAc                                          >  1:607720/1‑140 (MQ=255)                                                                                           ccGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAg                                      >  1:732698/1‑140 (MQ=255)                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCACCCACTGACGGGACGCCCGCTTCCCGCAGGCAACTTAgcgc                                   >  1:866057/1‑140 (MQ=255)                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgcgc                                   >  2:180182/1‑140 (MQ=255)                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgcgc                                   <  1:770954/140‑1 (MQ=255)                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAgcgc                                   <  1:180182/140‑1 (MQ=255)                                                                                                            gcTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgc                     <  2:713184/140‑1 (MQ=255)                                                                                                             cTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcg                    <  2:607720/140‑1 (MQ=255)                                                                                                             cTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcg                    <  1:733977/140‑1 (MQ=255)                                                                                                             cTGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCATgcg                    >  1:121682/1‑140 (MQ=255)                                                                                                              tGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                   >  1:569811/1‑140 (MQ=255)                                                                                                              tGATGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                   <  2:767093/140‑1 (MQ=255)                                                                                                               gATGCCTTCTTCATTGAGCGCGCGGGTCTGGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGcgcg                  <  1:300722/140‑1 (MQ=255)                                                                                                                aTGCCTTCTTCATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGcgt                 >  1:738434/1‑140 (MQ=255)                                                                                                                        ttcATTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAc         >  1:91706/1‑140 (MQ=255)                                                                                                                           aTTGAGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTg      >  1:11614/1‑140 (MQ=255)                                                                                                                               aGCGCGCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTc  >  2:73798/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       ATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTC  >  NC_000913/4559033‑4559298

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