BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,151	C→T	27.4%	D352D (GAC→GAT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,151	0	C	T	27.4%	238.3 / 100.9	153	D352D (GAC→GAT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (49/62); new base T (30/12); total (79/74)
Fisher's exact test for biased strand distribution p-value = 3.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 157 total aligned reads displayed.

AACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATT  >  minE/2571081‑2571221

aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                                                        <  1:2608988/71‑1 (MQ=255)

     gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                                                   <  1:887959/71‑1 (MQ=255)

        ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                                                        <  1:3077075/63‑1 (MQ=255)

            aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                             <  1:62860/70‑1 (MQ=255)

            aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                            >  1:1023888/1‑71 (MQ=255)

            aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                            >  1:4710826/1‑71 (MQ=255)

            aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                            >  1:18021/1‑71 (MQ=255)

            aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                            >  1:1613540/1‑71 (MQ=255)

                tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                        >  1:4241635/1‑71 (MQ=255)

                  atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                             >  1:200360/1‑64 (MQ=255)

                   tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                     <  1:3286039/71‑1 (MQ=255)

                      cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                  >  1:222200/1‑71 (MQ=255)

                       gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                                 <  1:2769342/71‑1 (MQ=255)

                       gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                                 >  1:4753995/1‑71 (MQ=255)

                       gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                                 >  1:3172167/1‑71 (MQ=255)

                        aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATg                                                >  1:5314680/1‑71 (MQ=255)

                         tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                        <  1:1158953/62‑1 (MQ=255)

                         tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                               <  1:2141511/71‑1 (MQ=255)

                           ctttACCAGACGGTATTTACGACCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                                              >  1:4595646/3‑70 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:2265140/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:3007160/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:4756545/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:109937/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:2489517/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:1555635/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:2004031/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:3948811/71‑1 (MQ=255)

                            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:4500377/71‑1 (MQ=255)

                             ttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            >  1:61714/1‑68 (MQ=255)

                              tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            >  1:2442065/1‑67 (MQ=255)

                              tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            >  1:2958551/1‑67 (MQ=255)

                              tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            <  1:4370586/69‑3 (MQ=255)

                               aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg                                         >  1:4422962/1‑71 (MQ=255)

                               aCCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCgg                                                     >  1:354773/1‑59 (MQ=255)

                                ccagACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                                                 >  1:3602603/4‑62 (MQ=255)

                                ccagACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                                                 >  1:2823070/4‑62 (MQ=255)

                                ccagACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                                                 >  1:3328451/4‑62 (MQ=255)

                                 ccGACGGCATCTACGCCCAGGTGTCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:3640913/58‑1 (MQ=255)

                                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:3453689/58‑1 (MQ=255)

                                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:101163/58‑1 (MQ=255)

                                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:4489800/58‑1 (MQ=255)

                                  cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                      >  1:3225277/1‑71 (MQ=255)

                                   gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:1387534/56‑1 (MQ=255)

                                     cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                   >  1:1634113/1‑69 (MQ=255)

                                     cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                   >  1:5093293/1‑69 (MQ=255)

                                     cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                            <  1:1933922/62‑1 (MQ=255)

                                         aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                  >  1:1785017/1‑68 (MQ=255)

                                         aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                               <  1:247959/71‑1 (MQ=255)

                                          tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                               >  1:3844818/1‑70 (MQ=255)

                                           cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                             >  1:3046096/1‑71 (MQ=255)

                                           cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                             >  1:5194102/1‑71 (MQ=255)

                                            tACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                   >  1:2071736/1‑62 (MQ=255)

                                            tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGTGGAAAGTGCTGCTAAGAAt                            >  1:5332413/1‑71 (MQ=255)

                                            tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                            >  1:1936429/1‑71 (MQ=255)

                                            tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                            >  1:1080029/1‑71 (MQ=255)

                                             aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                              >  1:5349075/1‑68 (MQ=255)

                                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                             <  1:4346162/68‑1 (MQ=255)

                                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                             <  1:4466867/68‑1 (MQ=255)

                                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                              <  1:3103857/67‑1 (MQ=255)

                                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                          >  1:3375289/1‑71 (MQ=255)

                                                cccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                        <  1:5510290/71‑1 (MQ=255)

                                                   aGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGt                                                    <  1:3003591/40‑1 (MQ=37)

                                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGa                                                  >  1:261139/1‑41 (MQ=37)

                                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                    <  1:2763793/71‑1 (MQ=255)

                                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                    >  1:622532/1‑71 (MQ=255)

                                                    ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                             <  1:3166395/46‑1 (MQ=255)

                                                    ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGt                                        >  1:3931399/1‑51 (MQ=255)

                                                     gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                   <  1:3993103/71‑1 (MQ=255)

                                                     gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                   <  1:4966286/71‑1 (MQ=255)

                                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:3559435/38‑1 (MQ=255)

                                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:2550300/38‑1 (MQ=255)

                                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                    <  1:2095646/38‑1 (MQ=255)

                                                       ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                          >  1:837368/1‑46 (MQ=255)

                                                       ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                          >  1:5106136/1‑46 (MQ=255)

                                                       ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTggg                 >  1:1178487/1‑70 (MQ=255)

                                                         cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                            <  1:5425119/58‑1 (MQ=255)

                                                          ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGc              <  1:740455/71‑1 (MQ=255)

                                                          ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGc              <  1:1867771/71‑1 (MQ=255)

                                                           gTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                            <  1:3804708/40‑3 (MQ=255)

                                                           gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                           <  1:4706375/41‑1 (MQ=255)

                                                            tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt            >  1:1102715/1‑71 (MQ=255)

                                                            tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt            >  1:5134461/1‑71 (MQ=255)

                                                            tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt            <  1:1072992/71‑1 (MQ=255)

                                                             ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                    >  1:4916584/1‑62 (MQ=255)

                                                             ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt            >  1:5544141/1‑70 (MQ=255)

                                                              tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                  >  1:4999038/1‑63 (MQ=255)

                                                              tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCAtt           <  1:2961591/70‑1 (MQ=255)

                                                              tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTa          >  1:1495673/1‑71 (MQ=255)

                                                              tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt           <  1:5006292/70‑1 (MQ=255)

                                                              tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt           <  1:3043775/70‑1 (MQ=255)

                                                               aCTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTac         >  1:2565042/1‑70 (MQ=255)

                                                               aCTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTac         >  1:478737/1‑70 (MQ=255)

                                                               aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt         >  1:682438/1‑71 (MQ=255)

                                                                ctctGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                      >  1:2373363/1‑41 (MQ=255)

                                                                  ccGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt       >  1:5092416/1‑70 (MQ=255)

                                                                   cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                     >  1:3770922/1‑55 (MQ=255)

                                                                    gACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa    <  1:3099256/71‑1 (MQ=255)

                                                                     aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                          >  1:589001/1‑48 (MQ=255)

                                                                      ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt  >  1:4091093/1‑71 (MQ=255)

                                                                      ccATTTCAGCCATGTTTCGGGGATGCGGAAAGTGCTGCTAAGAATGACCAGc                     <  1:2560631/52‑1 (MQ=255)

AACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATT  >  minE/2571081‑2571221

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