BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,589	T→A	100%	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,589	0	T	A	100.0%	136.3 / NA	48	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (26/22); total (26/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563523‑4563654

gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:3747503/71‑1 (MQ=255)

 tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:4068751/71‑1 (MQ=255)

 tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:4341473/71‑1 (MQ=255)

 tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:4346712/71‑1 (MQ=255)

   ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            <  1:1626527/71‑1 (MQ=255)

    gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                           <  1:4028250/71‑1 (MQ=255)

        cagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                       >  1:3409840/1‑71 (MQ=255)

           cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:118323/71‑1 (MQ=255)

            ttGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   >  1:4059117/1‑71 (MQ=255)

                tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:3350658/1‑66 (MQ=255)

                  aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:439960/69‑1 (MQ=255)

                  aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:2340433/1‑71 (MQ=255)

                          aTCTCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              <  1:595755/46‑1 (MQ=39)

                            cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   >  1:4357963/1‑71 (MQ=255)

                            cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                   >  1:3617940/1‑71 (MQ=255)

                             gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            >  1:1883247/1‑45 (MQ=255)

                             gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                            >  1:3719553/1‑45 (MQ=255)

                               aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                <  1:1155371/71‑1 (MQ=255)

                                gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:1601509/1‑71 (MQ=255)

                                 ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                              >  1:491162/1‑71 (MQ=255)

                                 ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                              >  1:3811078/1‑71 (MQ=255)

                                 ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                              >  1:357655/1‑71 (MQ=255)

                                  tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                             >  1:1684814/1‑71 (MQ=255)

                                    tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:849996/70‑1 (MQ=255)

                                      cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:1719248/49‑1 (MQ=255)

                                        gcagGTCATAAATAGGCTGTAACGCCATTTTACTTTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                       <  1:3422740/71‑1 (MQ=255)

                                        gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:1543136/1‑49 (MQ=255)

                                          agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:193739/71‑1 (MQ=255)

                                             tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:2599658/1‑37 (MQ=255)

                                              cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:469688/1‑70 (MQ=255)

                                               aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                <  1:2613439/71‑1 (MQ=255)

                                                tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACTCCAGGTTTACCGCCGAGCa                <  1:2756297/70‑1 (MQ=255)

                                                tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                           <  1:575061/59‑1 (MQ=255)

                                                tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgc                      <  1:3117606/64‑1 (MQ=255)

                                                 aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                       >  1:320666/1‑62 (MQ=255)

                                                 aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                       >  1:4147594/1‑62 (MQ=255)

                                                 aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa              >  1:836102/1‑71 (MQ=255)

                                                  aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:1381034/1‑53 (MQ=255)

                                                  aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             >  1:2336379/1‑71 (MQ=255)

                                                   aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:2961066/71‑1 (MQ=255)

                                                     aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGat                                             >  1:290769/1‑35 (MQ=38)

                                                      ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:3194708/71‑1 (MQ=255)

                                                         tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCTCCGAGCAAACCGCCAAc      >  1:4072401/1‑71 (MQ=255)

                                                           tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    <  1:204641/71‑1 (MQ=255)

                                                           tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:1350958/1‑71 (MQ=255)

                                                           tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:2794903/1‑71 (MQ=255)

                                                           tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcg     >  1:4273623/1‑70 (MQ=255)

                                                             aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:1840704/71‑1 (MQ=255)

                                                             aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:1327518/71‑1 (MQ=255)

GTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563523‑4563654

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