Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,559,680	G→A	G77G (GGC→GGT)	ddpC ←	D,D‑dipeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,559,680	0	G	A	100.0%	88.0 / NA	28	G77G (GGC→GGT)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/16);  total (12/16)

ACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAA  >  NC_000913/1559545‑1559792                                                                                                                                        |                                                                                                                 accaAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCa                                                                                                               <  1:764448/139‑1 (MQ=255)         cacCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTg                                                                                                       <  2:174198/139‑1 (MQ=255)                tAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCa                                                                                                <  1:300196/139‑1 (MQ=255)                                                 aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                               >  1:359747/1‑139 (MQ=255)                                                  aTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                             >  2:753605/1‑124 (MQ=255)                                                  aTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                             <  1:753605/124‑1 (MQ=255)                                                            cGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                 >  2:499212/1‑126 (MQ=255)                                                            cGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                 <  1:499212/126‑1 (MQ=255)                                                                       gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgcca                                         >  2:928492/1‑139 (MQ=255)                                                                       gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgcca                                         >  1:891216/1‑139 (MQ=255)                                                                         gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgccagc                                       <  1:928492/139‑1 (MQ=255)                                                                           gaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgccagcca                                     <  2:359747/139‑1 (MQ=255)                                                                                 attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                             <  1:864405/93‑1 (MQ=255)                                                                                 attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                             >  2:864405/1‑93 (MQ=255)                                                                                           gCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAatcatc                     >  2:443265/1‑139 (MQ=255)                                                                                                   gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg             <  1:897855/139‑1 (MQ=255)                                                                                                   gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg             <  1:947981/139‑1 (MQ=255)                                                                                                      cgcgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt          <  2:891216/139‑1 (MQ=255)                                                                                                      cgcgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt          <  1:443265/139‑1 (MQ=255)                                                                                                             aaaCAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt               <  2:34364/127‑1 (MQ=255)                                                                                                             aaaCAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt               >  1:34364/1‑127 (MQ=255)                                                                                                              aaCAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCa           <  1:891835/130‑1 (MQ=255)                                                                                                              aaCAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACaa  >  1:295573/1‑139 (MQ=255)                                                                                                              aaCAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCa           >  2:891835/1‑130 (MQ=255)                                                                                                                     cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcg                                                                               <  2:586149/55‑1 (MQ=255)                                                                                                                     cgagTCCCACTACATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcg                                                                               >  1:586149/4‑55 (MQ=255)                                                                                                                        gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt               >  2:502454/1‑116 (MQ=255)                                                                                                                        gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt               <  1:502454/116‑1 (MQ=255)                                                                                                                                        |                                                                                                                 ACCAAGCACACCGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAA  >  NC_000913/1559545‑1559792

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