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%	156.2 / NA	53	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (27/26); total (27/26)
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:350991/71‑1 (MQ=255)

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

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

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

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

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

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

          agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        <  1:3114842/71‑1 (MQ=255)

          agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        <  1:1846395/71‑1 (MQ=255)

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

              tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                 >  1:381449/1‑58 (MQ=255)

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

                cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   >  1:1593073/1‑70 (MQ=255)

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

                  gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:3792057/1‑56 (MQ=255)

                  gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:3516193/1‑56 (MQ=255)

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

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

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

                          caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                        <  1:931641/71‑1 (MQ=255)

                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    >  1:1329684/1‑71 (MQ=255)

                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                           >  1:2380212/1‑47 (MQ=255)

                               cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    >  1:2375479/1‑70 (MQ=255)

                                aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   <  1:1635391/54‑1 (MQ=255)

                                  ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                                 >  1:695097/1‑70 (MQ=255)

                                   tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                          >  1:3979028/1‑60 (MQ=255)

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

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                          >  1:1112306/1‑59 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                               <  1:490527/70‑1 (MQ=255)

                                       cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                           >  1:3484800/1‑71 (MQ=255)

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

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

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

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

                                           agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgca                       <  1:86172/71‑2 (MQ=255)

                                           agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgc                        >  1:114261/1‑70 (MQ=255)

                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:3329728/52‑1 (MQ=255)

                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                      <  1:1281930/71‑1 (MQ=255)

                                            ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                      <  1:2811536/71‑1 (MQ=255)

                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:3398383/1‑53 (MQ=255)

                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:948681/57‑1 (MQ=255)

                                             gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                   <  1:2874819/57‑1 (MQ=255)

                                               cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                             >  1:2975988/1‑45 (MQ=255)

                                                aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:3403876/40‑1 (MQ=255)

                                                 tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                           >  1:791268/1‑61 (MQ=255)

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

                                                    aTAGGCTGAAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                          >  1:2864818/1‑43 (MQ=38)

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

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

                                                        gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:662133/1‑66 (MQ=255)

                                                        gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:903761/1‑66 (MQ=255)

                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:2618166/1‑59 (MQ=255)

                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:3885103/1‑59 (MQ=255)

                                                              aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:4017671/1‑60 (MQ=255)

                                                                gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCggg   >  1:1811603/1‑70 (MQ=255)

                                                                gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  >  1:1545609/1‑71 (MQ=255)

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656

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