Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,051	A→G	T162A (ACC→GCC)	iaaA →	isoaspartyl dipeptidase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,051	0	A	G	100.0%	121.2 / NA	33	T162A (ACC→GCC)	iaaA	isoaspartyl dipeptidase proenzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

AAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAC  >  NC_000913/866915‑867186                                                                                                                                         |                                                                                                                                        aaGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCg                                                                                                                                      <  2:390163/140‑1 (MQ=255) aaGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCg                                                                                                                                      >  1:792752/1‑140 (MQ=255)      cGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTc                                                                                                                                 >  1:695127/1‑140 (MQ=255)       gCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCg                                                                                                                                >  1:642108/1‑140 (MQ=255)        cATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGa                                                                                                                               >  1:784327/1‑140 (MQ=255)           gtgATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCa                                                                                                                            <  1:164193/140‑1 (MQ=255)              atgatgatTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGTAGCCGTCCTCGACCATAg                                                                                                                         >  2:594362/1‑140 (MQ=255)                  tgatTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGt                                                                                                                     >  1:259277/1‑140 (MQ=255)                   gatTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTg                                                                                                                    <  2:185837/140‑1 (MQ=255)                         cGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCAc                                                                                                              >  1:378852/1‑140 (MQ=255)                               ggCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATg                                                                                                        <  2:632458/140‑1 (MQ=255)                                  aGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGaaa                                                                                                     <  2:403953/140‑1 (MQ=255)                                     aaaTTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTACTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAAc                                                                                                  <  2:256036/140‑1 (MQ=255)                                               tttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATggg                                                                                        <  2:469229/140‑1 (MQ=255)                                               tttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATggg                                                                                        <  2:259277/140‑1 (MQ=255)                                                ttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGc                                                                                       <  1:779142/140‑1 (MQ=255)                                                ttGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGc                                                                                       <  1:249516/140‑1 (MQ=255)                                                                   cgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGTCGTGGCGt                                                                    <  2:642108/140‑1 (MQ=255)                                                                   cgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGt                                                                    <  2:105600/140‑1 (MQ=255)                                                                                                tGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGcc                                       >  1:227331/1‑140 (MQ=255)                                                                                                         aacaacTACTGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCcaca                              >  2:614653/1‑140 (MQ=255)                                                                                                                 ctGGCAGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAAt                      >  1:397645/1‑140 (MQ=255)                                                                                                                      aGCGCGCAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCa                 >  1:781860/1‑140 (MQ=255)                                                                                                                        cgcgcAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAAt               >  1:144648/1‑140 (MQ=255)                                                                                                                        cgcgcAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAAt               >  2:529156/1‑140 (MQ=255)                                                                                                                        cgcgcAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAAt               >  2:727146/1‑140 (MQ=255)                                                                                                                         gcgcAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa              >  2:414027/1‑140 (MQ=255)                                                                                                                         gcgcAAGGAAGGGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa              >  1:607082/1‑140 (MQ=255)                                                                                                                              aggaaggGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAAt                      <  2:563025/127‑1 (MQ=255)                                                                                                                              aggaaggGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAAt                      >  1:563025/1‑127 (MQ=255)                                                                                                                                 aaggGGCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTAccc      <  2:784327/140‑1 (MQ=255)                                                                                                                                     ggCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGGCGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAc  >  1:814416/1‑140 (MQ=255)                                                                                                                                     ggCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAc  <  1:810530/140‑1 (MQ=255)                                                                                                                                     ggCAGCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAc  >  2:810530/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        AAGCCCGCATGTGATGATGATTGGCGAAGGGGCAGAAAATTTTGCGTTTGCTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAC  >  NC_000913/866915‑867186

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