BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,193	A→C	25.2%	R366S (AGA→AGC) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,193	0	A	C	25.2%	202.9 / 109.2	151	R366S (AGA→AGC) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (52/61); new base C (20/18); total (72/79)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 161 total aligned reads displayed.

CTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTT  >  minE/2571124‑2571258

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

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

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

   cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGGTAAGaa                                                                  <  1:10340/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                   tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                                                       >  1:5479730/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)

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

                          atcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCACTACCttt                                          <  1:1088495/67‑1 (MQ=255)

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

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

                           tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTg                                                    <  1:1535765/57‑1 (MQ=255)

                           tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                                       <  1:3609129/70‑2 (MQ=255)

                           tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                                       <  1:5010097/70‑2 (MQ=255)

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

                            cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAcc                                                             >  1:1493427/1‑48 (MQ=255)

                            cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAcc                                                             >  1:4489153/1‑48 (MQ=255)

                            cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAcc                                                             >  1:3192962/1‑48 (MQ=255)

                            cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtct                                      >  1:2350337/1‑69 (MQ=255)

                             aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                        <  1:258812/52‑1 (MQ=255)

                             aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                     <  1:1718755/71‑1 (MQ=255)

                             aTTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTaac                                       >  1:5161912/1‑67 (MQ=255)

                             aTTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtctc                                     >  1:1381078/1‑68 (MQ=255)

                              tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                  >  1:3620312/1‑41 (MQ=255)

                              tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                      >  1:5674979/1‑53 (MQ=255)

                              tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                      >  1:2274517/1‑53 (MQ=255)

                              tttCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                                                            <  1:3494986/47‑1 (MQ=255)

                                tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                  >  1:472049/1‑71 (MQ=255)

                                 cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                  <  1:956006/70‑1 (MQ=255)

                                  aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                                <  1:1321618/71‑1 (MQ=255)

                                  aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                                <  1:2616552/71‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                       <  1:377226/62‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                       <  1:3158190/62‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              >  1:4362758/1‑71 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              <  1:1570369/71‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              <  1:1733253/71‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              <  1:5384325/71‑1 (MQ=255)

                                    ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                               >  1:3999890/1‑70 (MQ=255)

                                     cATGTTTCGGTGATTCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                             <  1:5186671/71‑1 (MQ=255)

                                     cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                             <  1:5332099/71‑1 (MQ=255)

                                     cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                             <  1:5103990/71‑2 (MQ=255)

                                      aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                        <  1:4807220/43‑1 (MQ=255)

                                      aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                        <  1:4100801/43‑1 (MQ=255)

                                      aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              >  1:603299/1‑69 (MQ=255)

                                      aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                              >  1:3875052/1‑69 (MQ=255)

                                      aTGTTTCGGTGATGCGGAAAGTCCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                     >  1:1980416/1‑62 (MQ=255)

                                       tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                    <  1:4604329/62‑1 (MQ=255)

                                       tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                    <  1:4432622/62‑1 (MQ=255)

                                       tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                           <  1:3632032/71‑1 (MQ=255)

                                        gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                             <  1:2647854/36‑1 (MQ=255)

                                        gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                            >  1:2894390/1‑69 (MQ=255)

                                        gTTTAGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                             <  1:2705178/36‑1 (MQ=255)

                                          ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                     >  1:924756/1‑58 (MQ=255)

                                          ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                            >  1:665949/1‑67 (MQ=255)

                                           tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                      >  1:2804004/1‑40 (MQ=255)

                                           tCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCAtt                                                >  1:875924/1‑46 (MQ=255)

                                            cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                       <  1:3658348/70‑1 (MQ=255)

                                            cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                        >  1:4349559/1‑69 (MQ=255)

                                             ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                        <  1:2341009/36‑1 (MQ=255)

                                             ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                        <  1:2767349/36‑1 (MQ=255)

                                             ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                      <  1:42438/70‑1 (MQ=255)

                                              gTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagg                    <  1:5671021/71‑1 (MQ=255)

                                                 ctGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                            <  1:5495212/43‑1 (MQ=255)

                                                 aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                  <  1:3269735/54‑1 (MQ=255)

                                                 aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                  <  1:1691149/54‑1 (MQ=255)

                                                   gcgcAAAGTGCTATTACGCATGACCAGCTGGGTGCAt                                                 >  1:3719673/1‑37 (MQ=255)

                                                   gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                      <  1:5496116/48‑1 (MQ=255)

                                                       aaaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt            >  1:616329/1‑70 (MQ=255)

                                                        aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa          >  1:249509/1‑71 (MQ=255)

                                                        aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa          >  1:2522013/1‑71 (MQ=255)

                                                            gctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                       >  1:4843058/1‑54 (MQ=255)

                                                             ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGa      >  1:4544095/1‑70 (MQ=255)

                                                              tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAc    >  1:4026501/1‑71 (MQ=255)

                                                               gctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                <  1:724465/42‑1 (MQ=255)

CTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTT  >  minE/2571124‑2571258

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