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

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656

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

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

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

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

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

  tCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                 >  1:4051525/1‑70 (MQ=255)

   cGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:61483/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                              >  1:4197943/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                              >  1:3888442/1‑71 (MQ=255)

      cagcagCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                            <  1:2234325/71‑1 (MQ=255)

           gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                       <  1:3625587/71‑1 (MQ=255)

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

              tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    >  1:1270009/1‑71 (MQ=255)

              tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                    >  1:3433425/1‑71 (MQ=255)

                cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                  <  1:2948004/71‑1 (MQ=255)

                cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   <  1:1290453/70‑1 (MQ=255)

                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                >  1:1189954/1‑54 (MQ=255)

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

                   aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               <  1:3463713/71‑1 (MQ=255)

                    tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               <  1:1241585/70‑1 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 >  1:4597133/1‑67 (MQ=255)

                      ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                            <  1:2379461/71‑1 (MQ=255)

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

                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:1795471/62‑1 (MQ=255)

                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:1313442/62‑1 (MQ=255)

                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:259352/70‑1 (MQ=255)

                           aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                 <  1:3272020/45‑1 (MQ=255)

                            tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:4746437/60‑1 (MQ=255)

                            tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:1737791/60‑1 (MQ=255)

                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:1199219/1‑45 (MQ=255)

                             cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:2265230/1‑45 (MQ=255)

                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:3475771/1‑44 (MQ=255)

                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:2513271/1‑44 (MQ=255)

                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:3642375/71‑1 (MQ=255)

                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:1907100/71‑1 (MQ=255)

                                 gTTTTCCAGCAGGTCATAAATAGGCTTTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 <  1:3384335/71‑1 (MQ=255)

                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                              <  1:2944127/41‑1 (MQ=255)

                                   tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                >  1:3108133/1‑38 (MQ=255)

                                   tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                >  1:570849/1‑38 (MQ=255)

                                   tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                            >  1:519679/1‑42 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:1180902/1‑62 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                              >  1:1089045/1‑71 (MQ=255)

                                      ccAGCAGGTCATAAATAGGCTATAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                             <  1:1168041/70‑1 (MQ=255)

                                          cagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgc                        >  1:1106802/1‑71 (MQ=255)

                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:1518751/43‑1 (MQ=255)

                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                     <  1:3244573/71‑1 (MQ=255)

                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:3374455/39‑1 (MQ=255)

                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                  >  1:826673/1‑70 (MQ=255)

                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                  >  1:1470428/1‑70 (MQ=255)

                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:634157/38‑1 (MQ=255)

                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACg                                           <  1:3269820/44‑1 (MQ=255)

                                                  aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:1117764/1‑48 (MQ=255)

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

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             >  1:32354/1‑41 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                              >  1:3776369/1‑56 (MQ=255)

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

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               <  1:3580552/71‑1 (MQ=255)

                                                     tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:3488935/43‑1 (MQ=255)

                                                     tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc              >  1:1449232/1‑70 (MQ=255)

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

                                                          tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                                  >  1:4798103/1‑45 (MQ=255)

                                                          tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc              <  1:1221569/65‑1 (MQ=255)

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

                                                            tAACGCAATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc      >  1:618383/1‑71 (MQ=255)

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

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

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

                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg     <  1:216916/70‑1 (MQ=255)

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

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

                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    >  1:1129086/1‑71 (MQ=255)

                                                                gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  <  1:3010436/71‑1 (MQ=255)

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656

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