Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	866,902	T→C	L112P (CTG→CCG)	iaaA →	isoaspartyl dipeptidase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	866,902	0	T	C	100.0%	89.2 / NA	26	L112P (CTG→CCG)	iaaA	isoaspartyl dipeptidase proenzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/11);  total (15/11)

TAACGCCGGAATTGGCGCTGTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCTGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCG  >  NC_000913/866765‑867037                                                                                                                                          |                                                                                                                                        tAACGCCGGAATTGGCGCTGTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCccggccgg                                                                                                                                       >  1:429308/1‑140 (MQ=255)    cGCCGGAATTGGCGCTGTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGa                                                                                                                                    <  1:86457/140‑1 (MQ=255)         gAATTGGCGCTGTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAg                                                                                                                               >  1:545365/1‑140 (MQ=255)                    gTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACCCCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCa                                                                                                                    <  1:476596/140‑1 (MQ=255)                     tCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCAt                                                                                                                   <  2:545365/140‑1 (MQ=255)                         tACGCGTGATGAAACCCATGAATTGGACGGCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGtga                                                                                                               <  1:228545/140‑1 (MQ=255)                               tgatgaAACCCATGAGCTGGACCCCCGTGCGATGGATCGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGtgatgatga                                                                                                         <  1:325809/140‑1 (MQ=255)                                       cccATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGaa                                                                                                 <  2:464443/140‑1 (MQ=255)                                           tGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCTCAGGTGATGCTGATTGGCGAAgggg                                                                                             >  1:280621/1‑140 (MQ=255)                                                  gACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGaaaa                                                                                      >  1:445556/1‑140 (MQ=255)                                                  gACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGaaaa                                                                                      >  2:37769/1‑140 (MQ=255)                                                            tgATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCCTGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGttt                                                                            >  1:557100/1‑140 (MQ=255)                                                            tgATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGttt                                                                            >  1:656819/1‑140 (MQ=255)                                                                    ggTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTgg                                                                    >  2:597531/1‑140 (MQ=255)                                                                          acCCTGAAAGCAAGTGCCGTGGCGGGCGTTAGTCATCTGCGTAACCCGGTTCTTGCCGCCCGGGCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCACAAAATTTTGCGTTTGCTCGTGGCATGGa                                                              <  2:655528/140‑1 (MQ=255)                                                                           cccTGAAAGCCGGTGTGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGTCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAg                                                             >  2:626138/1‑140 (MQ=255)                                                                           cccTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAg                                                             >  2:621101/1‑140 (MQ=255)                                                                                      ggtgcggtgGCGGGCATTAGTCGTCAGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGcc                                                  <  1:310547/140‑1 (MQ=255)                                                                                       gtgcggtgGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCg                                                 <  2:494697/140‑1 (MQ=255)                                                                                       gtgcggtgGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCg                                                 <  2:557100/140‑1 (MQ=255)                                                                                             tgGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGAtt                                           <  1:37769/140‑1 (MQ=255)                                                                                                  ggcgTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTtctc                                      >  1:52259/1‑140 (MQ=255)                                                                                                             tCTGCGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGc                           >  2:401614/1‑140 (MQ=255)                                                                                                                 cGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTa                       >  1:451800/1‑140 (MQ=255)                                                                                                                 cGTAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTa                       >  2:560225/1‑140 (MQ=255)                                                                                                                   tAATCCGGTTCTTGCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATg                     <  1:83210/140‑1 (MQ=255)                                                                                                                           ttcttcCCGCCCGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAacta             >  2:40907/1‑140 (MQ=255)                                                                                                                                      cGGCCGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAgcgcg  >  2:69291/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        TAACGCCGGAATTGGCGCTGTCTTTACGCGTGATGAAACCCATGAACTGGACGCCTGTGTGATGGATGGTAACACCCTGAAAGCCGGTGCGGTGGCGGGCGTTAGTCATCTGCGTAATCCGGTTCTTGCCGCCCGGCTGGTGATGGAGCAAAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCG  >  NC_000913/866765‑867037

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