BRESEQ :: Evidence

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%	178.6 / NA	46	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 (24/22); total (24/22)

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTT  >  NC_000913/1562039‑1562133

aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                <  1:3159591/49‑1 (MQ=255)

aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                <  1:227084/49‑1 (MQ=255)

aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                <  1:176334/49‑1 (MQ=255)

aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                               >  1:998790/1‑50 (MQ=255)

aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                               <  1:3583243/50‑1 (MQ=255)

  cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                             <  1:682690/50‑1 (MQ=255)

  cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                             <  1:349380/50‑1 (MQ=255)

  cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                             <  1:973208/50‑1 (MQ=255)

    tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                           <  1:1781202/50‑1 (MQ=255)

       gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                        >  1:952609/1‑50 (MQ=255)

        ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                       <  1:200914/50‑1 (MQ=255)

        ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                       <  1:3413557/50‑1 (MQ=255)

        ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                        <  1:3606046/49‑1 (MQ=255)

          cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                     >  1:2521812/1‑50 (MQ=255)

          cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGccc                                      <  1:2351825/49‑1 (MQ=255)

           gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                     <  1:1887695/49‑1 (MQ=255)

           gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                    <  1:3349912/50‑1 (MQ=255)

             cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                  >  1:676030/1‑50 (MQ=255)

             cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                  >  1:825892/1‑50 (MQ=255)

              gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                  <  1:2633357/49‑1 (MQ=255)

              gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                                 >  1:2112691/1‑50 (MQ=255)

               ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                                >  1:232222/1‑50 (MQ=255)

               ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                                <  1:2415112/50‑1 (MQ=255)

                   cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                            >  1:2771618/1‑50 (MQ=255)

                     tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  1:3314168/1‑50 (MQ=255)

                      cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  1:2674840/1‑49 (MQ=255)

                      cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  1:948913/1‑49 (MQ=255)

                      cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  1:2301312/1‑49 (MQ=255)

                       aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttttt                         <  1:354918/49‑1 (MQ=255)

                       aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                        <  1:3143744/50‑1 (MQ=255)

                        aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTGTtag                       <  1:2155425/50‑1 (MQ=255)

                         tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                       >  1:1085405/1‑49 (MQ=255)

                           tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    >  1:3539546/1‑50 (MQ=255)

                           tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    >  1:1766646/1‑50 (MQ=255)

                             tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                     <  1:3052946/47‑1 (MQ=255)

                             tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                     <  1:2362280/47‑1 (MQ=255)

                              aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg                 >  1:1910244/1‑50 (MQ=255)

                                aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               <  1:1068121/50‑1 (MQ=255)

                                   tttGTGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             >  1:823883/1‑49 (MQ=255)

                                   tttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             >  1:367419/1‑49 (MQ=255)

                                   tttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             >  1:1659444/1‑49 (MQ=255)

                                   tttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             >  1:2643303/1‑49 (MQ=255)

                                    ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa            >  1:3921914/1‑49 (MQ=255)

                                     tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa            <  1:3906183/48‑1 (MQ=255)

                                     tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             <  1:3553190/47‑1 (MQ=255)

                                         aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa      <  1:2707930/50‑1 (MQ=255)

                                         aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga        >  1:1274119/1‑48 (MQ=255)

                                           aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt    >  1:1989263/1‑50 (MQ=255)

                                           aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt    >  1:1565823/1‑50 (MQ=255)

                                              cTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGttt  <  1:3789436/49‑1 (MQ=255)

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTT  >  NC_000913/1562039‑1562133

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