BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,191	A→C	9.6%	R366R (AGA→CGA) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,191	0	A	C	9.6%	275.3 / 28.1	115	R366R (AGA→CGA) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (51/53); new base C (7/4); total (58/57)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 121 total aligned reads displayed.

CATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571121‑2571256

cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                   <  1:626791/71‑1 (MQ=255)

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

   cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                                 >  1:3994567/1‑70 (MQ=255)

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

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

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

      cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                               <  1:4080024/69‑1 (MQ=255)

      cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             <  1:3668258/71‑1 (MQ=255)

         ccAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                                                          <  1:80103/71‑1 (MQ=255)

            ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                   <  1:2728731/59‑1 (MQ=255)

             gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                                                      >  1:1095014/1‑71 (MQ=255)

             gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                                                      >  1:210843/1‑71 (MQ=255)

             gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      >  1:3945463/1‑71 (MQ=255)

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

                gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                               >  1:1617535/1‑59 (MQ=255)

                gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    <  1:952344/70‑1 (MQ=255)

                 cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                                  >  1:2814993/1‑71 (MQ=255)

                 cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                                  <  1:239614/71‑1 (MQ=255)

                  ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc                                                 <  1:3009714/71‑1 (MQ=255)

                  ggTTTACTCCGACCATTTCAGCCATGTTACGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:3312450/66‑1 (MQ=255)

                    tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                   <  1:1068004/51‑1 (MQ=255)

                    tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                   <  1:4142337/51‑1 (MQ=255)

                     ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             >  1:3380717/1‑56 (MQ=255)

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

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

                      tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                               >  1:3680293/1‑69 (MQ=255)

                      tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                             >  1:1855529/1‑71 (MQ=255)

                       aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                             <  1:253048/70‑1 (MQ=255)

                       aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                             <  1:248070/70‑1 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                  >  1:2717872/1‑45 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             >  1:2593160/1‑50 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                                       >  1:3403223/1‑56 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCCGCTAAGAATGACCAGCTGGATGCATTATCtt                                         >  1:2423530/1‑70 (MQ=255)

                               cATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                >  1:251483/1‑43 (MQ=255)

                               cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                                >  1:930626/1‑59 (MQ=255)

                                aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                                       <  1:1667984/67‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                  <  1:2863215/39‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:2049599/51‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                  <  1:679027/71‑1 (MQ=255)

                                   tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa                                 <  1:2342399/70‑1 (MQ=255)

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

                                     aGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                              <  1:2335195/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                                       <  1:457302/60‑1 (MQ=255)

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

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

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

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

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

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

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

                                        cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                         >  1:4082641/1‑57 (MQ=255)

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

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

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

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

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

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

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:295266/42‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:2223298/42‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:1544684/42‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                             <  1:2066690/67‑1 (MQ=255)

                                          tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                        <  1:2911709/56‑1 (MQ=255)

                                          tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                        <  1:2338841/56‑1 (MQ=255)

                                           gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                                        <  1:623439/39‑1 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                                       >  1:485237/1‑39 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:2596392/1‑42 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                    <  1:2604808/58‑1 (MQ=255)

                                            tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                                      >  1:2445816/3‑56 (MQ=255)

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

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

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                             >  1:3582595/1‑47 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            >  1:3999649/1‑64 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                             >  1:3111456/1‑63 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:594621/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:1448316/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:3462416/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:3936350/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:3833739/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      <  1:1567541/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg                       >  1:1463137/1‑69 (MQ=255)

                                              tCGGAGATGCGGAAAGTGCTGCTAAGAATGACCAg                                                         >  1:2288865/1‑35 (MQ=255)

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

                                                ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:3888278/1‑38 (MQ=255)

                                                  tGATGCGGAAAGTGCTGCTAAGAATGACCATCTGGATGCATTATCTTTATTTTTATCa                              <  1:2245708/58‑1 (MQ=255)

                                                   gATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagggg                >  1:1510288/1‑71 (MQ=255)

                                                     tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                              <  1:2250476/55‑1 (MQ=255)

                                                     tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                         <  1:409267/60‑1 (MQ=255)

                                                      gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATc                                           >  1:1592872/1‑41 (MQ=255)

                                                      gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                              >  1:3914078/1‑54 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                              <  1:3005181/35‑1 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt          >  1:3659699/1‑71 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt          >  1:295657/1‑71 (MQ=255)

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

                                                           aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            <  1:1692103/51‑1 (MQ=255)

                                                              tgctgcTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                               <  1:1206929/45‑1 (MQ=255)

                                                                ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa   >  1:1196068/1‑71 (MQ=255)

                                                                ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa   >  1:3857652/1‑71 (MQ=255)

                                                                 tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                           <  1:3705773/46‑1 (MQ=255)

                                                                    tAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                               <  1:785489/39‑1 (MQ=255)

                                                                     aaGAATGACCAGCT‑‑ATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAc  >  1:2990063/1‑65 (MQ=255)

CATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571121‑2571256

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