Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	408,431	A→G	E85E (GAA→GAG)	ppnP →	nucleoside phosphorylase PpnP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	408,431	0	A	G	100.0%	86.7 / NA	25	E85E (GAA→GAG)	ppnP	nucleoside phosphorylase PpnP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/14);  total (11/14)

GCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGAT  >  NC_000913/408294‑408560                                                                                                                                          |                                                                                                                                  gCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGcc                                                                                                                                 <  1:219281/140‑1 (MQ=255) gCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGcc                                                                                                                                 >  2:219281/1‑140 (MQ=255) gCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGcc                                                                                                                                 >  2:130562/1‑140 (MQ=255)  cTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGccc                                                                                                                                >  1:103122/1‑140 (MQ=255)  cTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGccc                                                                                                                                >  2:8540/1‑140 (MQ=255)   tGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCa                                                                                                                               <  2:103122/140‑1 (MQ=255)           aaGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTAt                                                                                                                       <  1:233556/140‑1 (MQ=255)           aaGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTAt                                                                                                                       <  1:78937/140‑1 (MQ=255)              agaTGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCtg                                                                                                                    <  1:91763/140‑1 (MQ=255)                   aCGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCg                                                                                                               <  1:130562/140‑1 (MQ=255)                   aCGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCg                                                                                                               <  2:261841/140‑1 (MQ=255)                                gcgcgcTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTAt                                                                                                           <  2:200637/131‑1 (MQ=255)                                gcgcgcTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTAt                                                                                                           >  1:200637/1‑131 (MQ=255)                                                                       aGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAat                                                           <  1:8540/140‑1 (MQ=255)                                                                                   ccGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTa                                               >  1:91219/1‑140 (MQ=255)                                                                                    cGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAg                                              >  1:372051/1‑140 (MQ=255)                                                                                    cGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTCTTCTGAAATATCCTTTCTTTac                                              >  1:253145/1‑139 (MQ=255)                                                                                    cGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTCTTCTGAAATATCCTTTCTTTAg                                              <  2:253145/140‑1 (MQ=255)                                                                                     ggTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGc                                             >  1:237623/1‑140 (MQ=255)                                                                                      gTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGcc                                            <  2:395423/140‑1 (MQ=255)                                                                                          ggTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCata                                        <  1:312159/140‑1 (MQ=255)                                                                                                tAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCataataata                                  >  2:276964/1‑140 (MQ=255)                                                                                                                         cATCTGCAAGTTGCCGAGCCCACCTCATATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAAt         <  1:276964/140‑1 (MQ=255)                                                                                                                           tCTGCAAGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAAttt       >  1:285192/1‑140 (MQ=255)                                                                                                                                aaGTTGCCGAGCCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGat  <  2:285192/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                  GCTGAGCCGGAAGAGATGACGGTAATCAGTGGCGCGCTGAATGTGTTACTGCCTGACGCGACCGACTGGCAGGTGTATGAAGCCGGTTCGGTGTTTAATGTTCCCGGTCACAGTGAGTTTCATCTGCAAGTTGCCGAACCCACCTCTTATCTGTGCCGCTATCTGTAATTCCTCGCCTTCCCCTTGAACGGGAGGGCATTTTTCTGAAATATCCTTTCTTTAGCCCATAATAATATTTCCTTTGCTGCGATTTTTTCAATTTCCGAT  >  NC_000913/408294‑408560

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