BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,178	G→C	25.7%	R361R (CGG→CGC)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,178	0	G	C	25.7%	309.0 / 115.8	179	R361R (CGG→CGC)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (51/82); new base C (27/19); total (78/101)
Fisher's exact test for biased strand distribution p-value = 2.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 183 total aligned reads displayed.

CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGG  >  minE/2571109‑2571240

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

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

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

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

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

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

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

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

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

      agACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                                         <  1:1571938/70‑1 (MQ=255)

       gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCAGAAAGtgc                                                         <  1:1849558/70‑1 (MQ=255)

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

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

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

               cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                                             >  1:3105923/1‑58 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                           ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                                     >  1:4524765/1‑70 (MQ=255)

                           ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                    >  1:3191083/1‑71 (MQ=255)

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

                              ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGc                                 <  1:740455/71‑1 (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:5544141/1‑70 (MQ=255)

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

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

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

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

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

                                     tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                      <  1:696623/59‑1 (MQ=255)

                                      ccGACCATTTCAGCCATTTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                         >  1:1245131/1‑71 (MQ=255)

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

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

                                       cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                         <  1:1151396/70‑1 (MQ=255)

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

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

                                          tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                     <  1:3609129/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)

                                           cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                                          >  1:887306/1‑49 (MQ=255)

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

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

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

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

                                               tCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGAt                >  1:5009621/1‑71 (MQ=255)

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

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

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

                                                   ccATGTTTCGGTGATGCGGAAATTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt            <  1:4346763/71‑1 (MQ=255)

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

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

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

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

                                                   ccATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtt            <  1:2614570/71‑1 (MQ=255)

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

                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                      <  1:3886935/44‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                        tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                      <  1:1476266/40‑1 (MQ=255)

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

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

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

                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga          >  1:4006483/1‑66 (MQ=255)

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

                                                           cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                          <  1:5112401/49‑1 (MQ=255)

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

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

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

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

                                                               gATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg       >  1:1555389/1‑64 (MQ=255)

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

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

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

                                                                 tGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta           <  1:4467704/58‑2 (MQ=255)

                                                                    ggAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga          <  1:4302414/56‑1 (MQ=255)

CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGG  >  minE/2571109‑2571240

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