BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,187	G→A	24.8%	L364L (CTG→CTA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,187	0	G	A	24.8%	267.0 / 88.1	153	L364L (CTG→CTA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (50/65); new base A (20/18); total (70/83)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 163 total aligned reads displayed.

ATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571122‑2571256

aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                                                     >  1:177315/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

              ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTggg                                                    >  1:1178487/1‑70 (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:1867771/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   <  1:4417888/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                         tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAttaa                          <  1:5176130/70‑3 (MQ=255)

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

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

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

                                            ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   >  1:4929438/1‑58 (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:875924/1‑46 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                   aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            <  1:2620025/58‑1 (MQ=255)

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

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

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

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

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

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

                                                             tgctgcTAAGAATGACCAGCTGGATGCATTATCttt                                        <  1:5405463/36‑1 (MQ=255)

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

                                                               ctgctAAGAATGACCAGCTGGATGCATTATCTTTa                                       >  1:5485272/1‑35 (MQ=255)

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

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

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

ATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571122‑2571256

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