BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,571,148	C→T	27.7%	S351S (TCC→TCT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,571,148	0	C	T	27.7%	227.2 / 98.4	148	S351S (TCC→TCT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (44/63); new base T (29/12); total (73/75)
Fisher's exact test for biased strand distribution p-value = 1.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 155 total aligned reads displayed.

TTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  Exported/2571078‑2571217

ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                                                       >  1:2639266/1‑71 (MQ=255)

  gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCg                                                                      <  1:4911662/70‑1 (MQ=255)

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

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

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

        gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                                                <  1:3070300/70‑1 (MQ=255)

               aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTc                                                              <  1:53430/65‑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:1975749/1‑71 (MQ=255)

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

                      tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                  <  1:5151667/70‑1 (MQ=255)

                        tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                <  1:3017175/70‑1 (MQ=255)

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

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

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

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

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

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

                            tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                           >  1:4482994/1‑71 (MQ=255)

                              ctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                                          >  1:2597122/3‑70 (MQ=255)

                              ctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                                          >  1:2556208/3‑70 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                               cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:1480555/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:2708292/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:3328451/4‑62 (MQ=255)

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

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

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

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

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

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

                                      gACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                  >  1:2934453/1‑70 (MQ=255)

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

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

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

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

                                             tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                             <  1:2807924/52‑1 (MQ=255)

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

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

                                               tACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCa                                                         >  1:1875563/1‑38 (MQ=38)

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

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

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

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

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

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

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

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

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

                                                  gCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                     <  1:4215956/71‑1 (MQ=255)

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

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

                                                       ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                <  1:2196252/71‑1 (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:5104158/38‑1 (MQ=255)

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

                                                        gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg               <  1:252901/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

TTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  Exported/2571078‑2571217

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