Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,802	A→G	L242P (CTG→CCG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,802	0	A	G	100.0%	166.7 / NA	46	L242P (CTG→CCG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/28);  total (18/28)

GTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTG  >  NC_000913/4558666‑4558936                                                                                                                                         |                                                                                                                                       gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGgcg                                                                                                                                     <  1:639561/140‑1 (MQ=255) gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGgcg                                                                                                                                     <  1:855976/140‑1 (MQ=255)         cGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGgcgc                                                                                                                                    >  2:885421/1‑133 (MQ=255)         cGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGgcgc                                                                                                                                    <  1:885421/133‑1 (MQ=255)          ggAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCaa                                                                                                                            >  1:796047/1‑140 (MQ=255)          ggAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCaa                                                                                                                            <  1:1264746/140‑1 (MQ=255)           gAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCaaa                                                                                                                           <  1:1334301/140‑1 (MQ=255)           gAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCaaa                                                                                                                           <  1:42648/140‑1 (MQ=255)                   ggTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTa                                                                                                                   <  1:716880/140‑1 (MQ=255)                      gACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGt                                                                                                                <  2:344905/140‑1 (MQ=255)                      gACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGt                                                                                                                <  2:925961/140‑1 (MQ=255)                                  cGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGt                                                                                                    <  2:796047/140‑1 (MQ=255)                                      aTACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGt                                                                                                <  2:1064467/140‑1 (MQ=255)                                      aTACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGt                                                                                                <  2:217084/140‑1 (MQ=255)                                          ccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGgc                                                                                            >  1:1293502/1‑140 (MQ=255)                                            gCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGgcag                                                                                          >  1:158329/1‑140 (MQ=255)                                               tGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGcagcag                                                                                       >  1:1183487/1‑140 (MQ=255)                                                   cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTg                                                                                   <  2:98526/140‑1 (MQ=255)                                                   cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTg                                                                                   >  1:121140/1‑140 (MQ=255)                                                   cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTg                                                                                   <  2:1183487/140‑1 (MQ=255)                                                   cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTg                                                                                   >  1:299472/1‑140 (MQ=255)                                                      cgcgGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTg                                                                                >  2:651564/1‑140 (MQ=255)                                                             aaTACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGca                                                                         <  2:299472/140‑1 (MQ=255)                                                                aCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACAt                                                                      <  2:121140/140‑1 (MQ=255)                                                                     cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAg                                                                 >  1:80918/1‑140 (MQ=255)                                                                     cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAg                                                                 >  1:1085936/1‑140 (MQ=255)                                                                     cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAg                                                                 >  1:1251266/1‑140 (MQ=255)                                                                     cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAg                                                                 >  2:419162/1‑140 (MQ=255)                                                                           gggCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTcc                                                           >  1:49084/1‑140 (MQ=255)                                                                              cGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCcagc                                                        <  2:759710/140‑1 (MQ=255)                                                                              cGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCcagc                                                        <  1:482415/140‑1 (MQ=255)                                                                              cGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCcagc                                                        <  2:530643/140‑1 (MQ=255)                                                                                        gTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATaaa                                              <  2:1293502/140‑1 (MQ=255)                                                                                                  ctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTaa                                    >  1:610314/1‑140 (MQ=255)                                                                                                    ggTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACg                                  <  2:287395/140‑1 (MQ=255)                                                                                                           tCGATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCttttt                           <  1:419162/140‑1 (MQ=255)                                                                                                             gATGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTAc                         <  2:158329/140‑1 (MQ=255)                                                                                                              aTGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACt                        <  2:552326/140‑1 (MQ=255)                                                                                                               tGGTGCCGCCTTTGCGCGCGAACTCCGGCGCCTACTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTg                       >  2:397281/1‑140 (MQ=255)                                                                                                                  tGCCGCCTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCg                    >  1:469360/1‑140 (MQ=255)                                                                                                                    ccgccTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGcc                  >  1:1232969/1‑140 (MQ=255)                                                                                                                    ccgccTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGcc                  >  1:566707/1‑140 (MQ=255)                                                                                                                        cTTTGCGCGCGAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCAtg              <  2:469360/140‑1 (MQ=255)                                                                                                                             cgcgcgAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGaa         <  2:794693/140‑1 (MQ=255)                                                                                                                             cgcgcgAACTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGaa         <  1:1290578/140‑1 (MQ=255)                                                                                                                                    aCTCCGGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCATGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTg  <  2:80918/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       GTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTG  >  NC_000913/4558666‑4558936

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