BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	494,006	T→A	D126E (GAT→GAA)	recR →	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	494,006	0	T	A	100.0%	200.8 / NA	70	D126E (GAT→GAA)	recR	gap repair protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (51/19); total (51/19)

TATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCAC  >  W3110S.gb/493941‑494068

tATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCt                                                           >  1:2316023/1‑71 (MQ=255)

 aTTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTg                                                          >  1:3079630/1‑71 (MQ=255)

 aTTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTg                                                          >  1:3794593/1‑71 (MQ=255)

    ttGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGa                                                        >  1:3157304/1‑70 (MQ=255)

                    tCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGt                                                             <  1:3238227/49‑1 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAAc                                                      >  1:3883294/1‑54 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGt                                                 >  1:1975707/1‑59 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGt                                                 >  1:1358020/1‑59 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAgg                                       >  1:3623290/1‑69 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:3887151/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:3824908/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:3723315/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:3589963/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:906276/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:1513462/1‑70 (MQ=255)

                      tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGa                                      >  1:1472013/1‑70 (MQ=255)

                       gTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaa                                     >  1:4477361/1‑70 (MQ=255)

                        tCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaa                                    >  1:3382630/1‑70 (MQ=255)

                          aCCGCTGGACTGCATCGGTCCGGATGATATCGGGCTTGAACg                                                              <  1:9135/42‑1 (MQ=38)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCag                                          >  1:1807659/1‑62 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                 >  1:207676/1‑71 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                 >  1:418326/1‑71 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                 >  1:2965206/1‑71 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaaa                                 >  1:2958624/1‑71 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                  <  1:1834939/70‑1 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                  >  1:3371697/1‑70 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                  >  1:3204424/1‑70 (MQ=255)

                          aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                  <  1:44305/70‑1 (MQ=255)

                            cGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGGCTGGAACAGCGTCTGGCAGAGGAAAAAATc                               <  1:3728666/71‑1 (MQ=255)

                             gCTGGACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCa                              >  1:3007952/1‑71 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                           >  1:422317/1‑70 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                           <  1:3751929/70‑1 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                           >  1:3598660/1‑70 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                           <  1:627767/70‑1 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg                           >  1:3593293/1‑70 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGa                          <  1:848976/71‑1 (MQ=255)

                                 gACGGCATCGGTCCGGATGATATCGGGCTTGAACGGCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGa                          <  1:3018627/71‑1 (MQ=255)

                                   cGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAg                        <  1:4399477/71‑1 (MQ=255)

                                   cGGCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAg                        >  1:23154/1‑71 (MQ=255)

                                    ggCATCGGTCCGGATGATATAGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGt                       >  1:2461820/1‑71 (MQ=255)

                                     gCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                      <  1:2682963/71‑1 (MQ=255)

                                     gCATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                      <  1:4265262/71‑1 (MQ=255)

                                      cATCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                     <  1:1905465/71‑1 (MQ=255)

                                       aTCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa                     >  1:3651849/1‑70 (MQ=255)

                                       aTCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGAt                    >  1:2316693/1‑71 (MQ=255)

                                        tCGGTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc                   <  1:4460168/71‑1 (MQ=255)

                                          ggTCCGTATGATATCGGGCTTGAACGTCTGGAACAGAGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:3451/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGTCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:1102389/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:552760/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:505765/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:2242730/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:1859657/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:2602863/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:292741/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:1655216/1‑71 (MQ=255)

                                          ggTCCGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                 >  1:2399667/1‑71 (MQ=255)

                                              cGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGaaaaa                                  >  1:2358725/1‑50 (MQ=255)

                                              cGGATGATATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGc              >  1:2519118/1‑70 (MQ=255)

                                                 atgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGccacc          <  1:3816205/71‑1 (MQ=255)

                                                 atgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGccac           <  1:2988000/70‑1 (MQ=255)

                                                  tgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCAc                             >  1:2886599/1‑51 (MQ=255)

                                                  tgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGccacc          >  1:2103282/1‑70 (MQ=255)

                                                  tgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGccacc          <  1:2378121/70‑1 (MQ=255)

                                                  tgatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGccacc          >  1:3120164/1‑70 (MQ=255)

                                                   gatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCa                              >  1:685038/1‑49 (MQ=255)

                                                   gatATCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGaa                         >  1:2998453/1‑54 (MQ=255)

                                                       tCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGa                          <  1:2362817/49‑1 (MQ=255)

                                                       tCGGGCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAccc     >  1:3257672/1‑70 (MQ=255)

                                                         gggCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCAc  >  1:32001/1‑71 (MQ=255)

                                                           gCTTGAACGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg                      <  1:3467150/49‑1 (MQ=255)

TATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCAC  >  W3110S.gb/493941‑494068

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