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%	105.9 / NA	35	G77G (GGC→GGT)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

AGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGT  >  NC_000913/1559570‑1559814                                                                                                               |                                                                                                                                       aGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCa                                                                                                            >  2:1176454/1‑139 (MQ=255) aGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCa                                                                                                            >  2:1014056/1‑139 (MQ=255)    aGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGTTACCAAACCAGTGCGCCGCAGAAGGCGGCaaaa                                                                                                         >  1:351287/1‑139 (MQ=255)         aCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAggcgg                                                                                                    <  2:351287/139‑1 (MQ=255)                      gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                       <  1:1014056/139‑1 (MQ=255)                      gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                       <  1:1283226/139‑1 (MQ=255)                        aaTGGCAACCACGACTAATCCGGCGAGAGTTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCa                                                                                                                                     <  2:825859/91‑1 (MQ=255)                        aaTGGCAACCACGACTAATCCGGCGAGAGTTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCa                                                                                                                                     >  1:825859/1‑91 (MQ=255)                        aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                                                     >  1:765421/1‑139 (MQ=255)                        aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                                                     >  1:1177234/1‑139 (MQ=255)                            gCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                                 >  1:346776/1‑139 (MQ=255)                            gCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGAACAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                                 <  2:853603/139‑1 (MQ=255)                                          tCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCATTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTc                                                                   >  2:1221597/1‑139 (MQ=255)                                          tCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTg                                                                                                                             >  1:639825/1‑81 (MQ=255)                                          tCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTg                                                                                                                             <  2:639825/81‑1 (MQ=255)                                              gCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgcca                                                               >  1:1036193/1‑139 (MQ=255)                                                gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgccagc                                                             <  1:1176454/139‑1 (MQ=255)                                                  gaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGCCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgccagcca                                                           <  2:228981/139‑1 (MQ=255)                                                  gaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCATCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCggcagaca                                                           <  2:1177234/139‑1 (MQ=255)                                                    attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGt                                                                                               <  1:1216126/101‑1 (MQ=255)                                                    attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGt                                                                                               >  2:1216126/1‑101 (MQ=255)                                                      tgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGt                                                       <  1:796447/139‑1 (MQ=255)                                                        attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCg                                                                                  <  1:305959/110‑1 (MQ=255)                                                        attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCg                                                                                  >  2:305959/1‑110 (MQ=255)                                                        attgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGa                                                     >  2:942567/1‑139 (MQ=255)                                                             tGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGaaaaaa                                                <  1:942567/139‑1 (MQ=255)                                                                      aCCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                       <  1:876491/91‑1 (MQ=255)                                                                      aCCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAACGCGGGCGGTTAAATCaa                                                                                       >  2:876491/1‑91 (MQ=255)                                                                             cgcgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt                                <  2:346776/139‑1 (MQ=255)                                                                                  cTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt                                >  2:1004271/1‑134 (MQ=255)                                                                                  cTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt                                <  1:1004271/134‑1 (MQ=255)                                                                                           tCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAAtgatga                  <  2:765421/139‑1 (MQ=255)                                                                                                    aCTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAAt                            <  2:15004/120‑1 (MQ=255)                                                                                                    aCTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAAt                            >  1:15004/1‑120 (MQ=255)                                                                                                           cGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGt  <  2:1036193/139‑1 (MQ=255)                                                                                                               |                                                                                                                                       AGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGT  >  NC_000913/1559570‑1559814

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