Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,562,085	C→T	A133T (GCA→ACA)	ddpA ←	putative D,D‑dipeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,562,085	0	C	T	100.0%	231.2 / NA	57	A133T (GCA→ACA)	ddpA	putative D,D‑dipeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/28);  total (29/28)

CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562038‑1562131                                                |                                               cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              >  1:1963933/1‑50 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:492391/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:876665/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:1622050/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:1440489/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             >  1:2875321/1‑50 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             >  1:1715156/1‑50 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             >  1:1153370/1‑50 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:4439583/50‑1 (MQ=255)    cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           <  1:4346690/50‑1 (MQ=255)       aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                        >  1:2398258/1‑50 (MQ=255)         gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                      >  1:723029/1‑50 (MQ=255)         gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                      >  1:1631469/1‑50 (MQ=255)         gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                      >  1:4300835/1‑50 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     >  1:2713983/1‑50 (MQ=255)            cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                   >  1:41815/1‑50 (MQ=255)            cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGccc                                    <  1:3840372/49‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:3336314/50‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:2634406/50‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:1719451/50‑1 (MQ=255)               cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                >  1:2257566/1‑50 (MQ=255)                gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                               <  1:2261900/50‑1 (MQ=255)                   aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                            >  1:1715258/1‑50 (MQ=255)                    gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                           >  1:3450651/1‑50 (MQ=255)                     cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                          >  1:2058473/1‑50 (MQ=255)                      aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         >  1:743926/1‑50 (MQ=255)                      aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         >  1:1145865/1‑50 (MQ=255)                      aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         <  1:1713698/50‑1 (MQ=255)                       tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:1852081/1‑50 (MQ=255)                       tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        <  1:3553139/50‑1 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttttt                       <  1:3351349/49‑1 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:2505058/1‑48 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:4049445/1‑48 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:3612689/1‑48 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                      <  1:2498922/50‑1 (MQ=255)                          aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                     <  1:2995876/50‑1 (MQ=255)                             tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  >  1:1772304/1‑50 (MQ=255)                              ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 >  1:1645544/1‑50 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  <  1:308546/48‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                   <  1:712580/47‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagt                    <  1:4067155/46‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 <  1:4068327/49‑1 (MQ=255)                               tAAATCTTTGGAAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  <  1:551951/48‑1 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                <  1:3483609/49‑1 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               >  1:880610/1‑50 (MQ=255)                                  aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt             >  1:2750333/1‑50 (MQ=255)                                    cTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc            >  1:1351327/1‑49 (MQ=255)                                      ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:253871/1‑49 (MQ=255)                                      ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:4431581/1‑49 (MQ=255)                                       tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaa         <  1:362566/49‑1 (MQ=255)                                       tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaa         <  1:30756/49‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa        <  1:1168118/49‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa        <  1:2378111/49‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:2307942/1‑47 (MQ=255)                                         ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      <  1:2529678/50‑1 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaa     <  1:853909/50‑1 (MQ=255)                                           aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa    <  1:1086495/50‑1 (MQ=255)                                            aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag   >  1:689433/1‑50 (MQ=255)                                            aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag   <  1:2430417/50‑1 (MQ=255)                                            aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag   <  1:371043/50‑1 (MQ=255)                                             aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  <  1:1178196/50‑1 (MQ=255)                                             aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  >  1:2695100/1‑50 (MQ=255)                                                |                                               CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562038‑1562131

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