Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,270,070	A→G	K219K (AAA→AAG)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,070	0	A	G	100.0%	133.2 / NA	37	K219K (AAA→AAG)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/21);  total (16/21)

CGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAC  >  NC_000913/4269946‑4270206                                                                                                                             |                                                                                                                                         cGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCg                                                                                                                           >  1:40160/1‑140 (MQ=255)  gCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCgg                                                                                                                          <  2:21390/140‑1 (MQ=255)    aaTGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGt                                                                                                                        <  2:710423/140‑1 (MQ=255)    aaTGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGt                                                                                                                        <  1:853265/140‑1 (MQ=255)    aaTGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCATGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGt                                                                                                                        <  2:165683/140‑1 (MQ=255)         cTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGtt                                                                                                                   <  2:368573/140‑1 (MQ=255)         cTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGtt                                                                                                                   <  1:542375/140‑1 (MQ=255)                      atatCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTg                                                                                                      <  2:168017/140‑1 (MQ=255)                         tCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATc                                                                                                   >  1:1018686/1‑140 (MQ=255)                           cGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGt                                                                                                 <  2:40160/140‑1 (MQ=255)                           cGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGt                                                                                                 <  1:488134/140‑1 (MQ=255)                                   ttATGGCGATCCTGATAATGATATTTCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAg                                                                                         <  2:527063/140‑1 (MQ=255)                                     aTGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                                                                                       <  2:835097/140‑1 (MQ=255)                                           aTTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGa                                                                                 <  1:571029/140‑1 (MQ=255)                                           aTTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGa                                                                                 <  1:258982/140‑1 (MQ=255)                                             tCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACa                                                                               <  2:1018686/140‑1 (MQ=255)                                                           tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCgg                                                                 >  1:166285/1‑139 (MQ=255)                                                            aCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTg                                                                >  1:28220/1‑140 (MQ=255)                                                               gccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                             <  1:251255/140‑1 (MQ=255)                                                                ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                            >  2:864331/1‑140 (MQ=255)                                                                ccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                            >  2:204959/1‑140 (MQ=255)                                                                 cgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                           >  1:214266/1‑140 (MQ=255)                                                                 cgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCACCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                           >  2:923710/1‑140 (MQ=255)                                                                    aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCtctc                                                        >  1:441036/1‑140 (MQ=255)                                                                                          aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                  <  1:923710/140‑1 (MQ=255)                                                                                          aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                  <  2:874809/140‑1 (MQ=255)                                                                                          aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                  <  1:906629/140‑1 (MQ=255)                                                                                          aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                  <  1:204959/140‑1 (MQ=255)                                                                                                 ttCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                           >  1:256876/1‑140 (MQ=255)                                                                                                    tGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                        >  1:601829/1‑140 (MQ=255)                                                                                                      cgcgcCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                      >  1:507031/1‑140 (MQ=255)                                                                                                           cTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGc                 >  2:625816/1‑140 (MQ=255)                                                                                                               aaCTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc             <  1:864331/140‑1 (MQ=255)                                                                                                                     acctacAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg       >  2:29775/1‑140 (MQ=255)                                                                                                                     acctacAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg       >  2:996713/1‑140 (MQ=255)                                                                                                                         acAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa   <  2:601829/140‑1 (MQ=255)                                                                                                                          cAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc  >  1:928049/1‑140 (MQ=255)                                                                                                                             |                                                                                                                                         CGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAC  >  NC_000913/4269946‑4270206

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: ‑