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%	162.6 / NA	56	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (33/23); total (33/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCC  >  W3110S.gb/4563521‑4563647

gggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatt                                                          >  1:1819786/1‑71 (MQ=255)

 ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                         <  1:994955/71‑1 (MQ=255)

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

        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                    >  1:4125166/1‑69 (MQ=255)

        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                  >  1:4119743/1‑71 (MQ=255)

        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                  <  1:1178466/71‑1 (MQ=255)

           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                               >  1:1601795/1‑71 (MQ=255)

           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                               >  1:3907179/1‑71 (MQ=255)

               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                             >  1:3735909/1‑69 (MQ=255)

                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                           <  1:599664/70‑1 (MQ=255)

                 cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATATGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                         <  1:510425/71‑1 (MQ=255)

                  tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                         >  1:4186158/1‑70 (MQ=255)

                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                            <  1:3031086/65‑1 (MQ=255)

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

                     tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                      >  1:3033468/1‑70 (MQ=255)

                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTTTCGCCCATGTGGAaca                                    >  1:3360180/1‑71 (MQ=255)

                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                   >  1:2876093/1‑71 (MQ=255)

                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                           <  1:4084290/62‑1 (MQ=255)

                        gCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                           <  1:3910644/62‑1 (MQ=255)

                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                <  1:1107532/71‑1 (MQ=255)

                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                          >  1:3002881/1‑60 (MQ=255)

                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACTCCatttt                                                        <  1:4560246/45‑1 (MQ=38)

                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                              <  1:4157413/71‑1 (MQ=255)

                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                               >  1:4236481/1‑51 (MQ=255)

                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                         >  1:3734977/1‑56 (MQ=255)

                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                         >  1:968032/1‑56 (MQ=255)

                                 aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                            >  1:2479332/1‑68 (MQ=255)

                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                        <  1:4599629/71‑1 (MQ=255)

                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                              >  1:653460/1‑48 (MQ=255)

                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                     >  1:109517/1‑39 (MQ=255)

                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                           >  1:414837/1‑49 (MQ=255)

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

                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                    <  1:4201964/70‑1 (MQ=255)

                                       ccAGCAGGGCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                    <  1:4519672/70‑1 (MQ=255)

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

                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                              >  1:3422259/1‑57 (MQ=255)

                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                >  1:2001/1‑71 (MQ=255)

                                           cagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                           >  1:1298365/1‑43 (MQ=255)

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

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

                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                      >  1:886745/1‑41 (MQ=255)

                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa        <  1:2706345/71‑1 (MQ=255)

                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa         <  1:2639762/70‑1 (MQ=255)

                                                    aaTAGGCTGTAACGCCATTTT‑CTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc     >  1:715312/1‑71 (MQ=255)

                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc     >  1:260311/1‑71 (MQ=255)

                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc     >  1:3229904/1‑71 (MQ=255)

                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                       <  1:1319402/52‑1 (MQ=255)

                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                  >  1:1340234/1‑39 (MQ=255)

                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCTGGTTTACCGCCGAGCAAACCGc   >  1:1186023/1‑70 (MQ=255)

                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACcccc  <  1:1967184/71‑1 (MQ=255)

                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc   >  1:3580203/1‑70 (MQ=255)

                                                         gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAc                             >  1:3133233/1‑43 (MQ=255)

                                                            gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc     <  1:2067982/64‑1 (MQ=255)

                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc     <  1:4012595/63‑1 (MQ=255)

                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                              >  1:2144819/1‑36 (MQ=255)

                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGAc                             <  1:1971306/37‑1 (MQ=255)

                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc   >  1:744457/1‑63 (MQ=255)

                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc   >  1:749759/1‑63 (MQ=255)

                                                                 gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                           >  1:554994/1‑37 (MQ=255)

                                                                 gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc  <  1:2514686/62‑1 (MQ=255)

GGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCC  >  W3110S.gb/4563521‑4563647

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