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%	364.9 / NA	91	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 (43/48);  total (43/48)

CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562038‑1562131                                                |                                               cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              >  1:796775/1‑50 (MQ=255) cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              >  1:1335295/1‑50 (MQ=255)  aacGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             <  1:1081442/48‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:971173/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:1652332/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             <  1:3769036/50‑1 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:1106731/50‑1 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:2390649/50‑1 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            >  1:4096460/1‑50 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:3692456/50‑1 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            >  1:1847651/1‑50 (MQ=255)    cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           >  1:624049/1‑50 (MQ=255)    cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           <  1:372124/50‑1 (MQ=255)    cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           >  1:1464893/1‑50 (MQ=255)     gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                          <  1:4493094/50‑1 (MQ=255)       aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                        >  1:5981920/1‑50 (MQ=255)       aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                        >  1:638671/1‑50 (MQ=255)        tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                        <  1:5414686/49‑1 (MQ=255)        tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                        <  1:4688768/49‑1 (MQ=255)        tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                       <  1:1870186/50‑1 (MQ=255)        tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                       >  1:5767308/1‑50 (MQ=255)        tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                       <  1:3457773/50‑1 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     <  1:2690155/50‑1 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     >  1:6026472/1‑50 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     <  1:3427937/50‑1 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     <  1:1581030/50‑1 (MQ=255)            cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGccc                                    <  1:3528234/49‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                   <  1:1167818/49‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:3284351/50‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:3529624/50‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:2598975/50‑1 (MQ=255)             gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                  <  1:4763459/50‑1 (MQ=255)              tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                 <  1:6042329/50‑1 (MQ=255)              tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                 <  1:3104622/50‑1 (MQ=255)              tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                 >  1:5256007/1‑50 (MQ=255)              tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                 <  1:754421/50‑1 (MQ=255)              tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                 <  1:4338521/50‑1 (MQ=255)               cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                >  1:3068899/1‑50 (MQ=255)                gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                <  1:3724459/49‑1 (MQ=255)                gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                               >  1:2235828/1‑50 (MQ=255)                  gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                             <  1:3966029/50‑1 (MQ=255)                  gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                             <  1:4897154/50‑1 (MQ=255)                  gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                             <  1:5785096/50‑1 (MQ=255)                   aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                            >  1:3350989/1‑50 (MQ=255)                   aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                            >  1:2725897/1‑50 (MQ=255)                   aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                            >  1:591549/1‑50 (MQ=255)                    gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                           >  1:4719227/1‑50 (MQ=255)                    gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                           >  1:1222294/1‑50 (MQ=255)                     cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                          >  1:3402837/1‑50 (MQ=255)                      aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         >  1:2109850/1‑50 (MQ=255)                      aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         <  1:4366795/50‑1 (MQ=255)                       tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:3455935/1‑50 (MQ=255)                       tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        <  1:2766528/50‑1 (MQ=255)                       tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:1889364/1‑50 (MQ=255)                        cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:5518148/1‑49 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:2090473/1‑48 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:186953/1‑48 (MQ=255)                          aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                     >  1:3605253/1‑50 (MQ=255)                             tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  >  1:958734/1‑50 (MQ=255)                              ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 >  1:3028716/1‑50 (MQ=255)                              ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 >  1:4841357/1‑50 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  <  1:5334056/48‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  <  1:1453505/48‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  <  1:1296994/48‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                   <  1:1080143/47‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                   <  1:655766/47‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                   <  1:2583609/47‑1 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                >  1:5880868/1‑50 (MQ=255)                               tAAATCTTTGGGAAATGTTTCTGCTGGCCCATGGCCGATTTTtagta                   <  1:5207441/47‑1 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               >  1:5197288/1‑50 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               >  1:689408/1‑50 (MQ=255)                                 aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               <  1:1207602/49‑1 (MQ=255)                                 aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               <  1:1236312/49‑1 (MQ=255)                                 aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               <  1:3720359/49‑1 (MQ=255)                                  aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt             <  1:5840632/50‑1 (MQ=255)                                    cTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc            >  1:5482992/1‑49 (MQ=255)                                      ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:5606227/1‑49 (MQ=255)                                       tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa           <  1:2414803/47‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa        <  1:2425449/49‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa        <  1:599225/49‑1 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:2269890/1‑47 (MQ=255)                                         ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag       >  1:1700663/1‑49 (MQ=255)                                         ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag       <  1:2536736/49‑1 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaa     <  1:2024971/50‑1 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaa     <  1:2461116/50‑1 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:2730026/1‑49 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:4443507/1‑49 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:447541/1‑49 (MQ=255)                                           aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa    <  1:5267113/50‑1 (MQ=255)                                           aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:1023435/1‑48 (MQ=255)                                           aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:3549615/1‑48 (MQ=255)                                           aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:5247691/1‑48 (MQ=255)                                            aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag   >  1:492803/1‑50 (MQ=255)                                             aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  >  1:2294935/1‑50 (MQ=255)                                             aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  >  1:4912730/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: ‑