BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,206	A→G	18.7%	M371V (ATG→GTG)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,206	0	A	G	18.7%	254.2 / 50.5	123	M371V (ATG→GTG)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (50/50); new base G (11/12); total (61/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 128 total aligned reads displayed.

GGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAC  >  minE/2571136‑2571274

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

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

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

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

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

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

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

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

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

       tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTa                                                               >  1:1495673/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                      >  1:5674979/1‑53 (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)

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

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

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

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

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

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

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

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

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

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

                         cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                                             <  1:5103990/71‑2 (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:4432622/62‑1 (MQ=255)

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

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

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

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

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

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

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

                              ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                            >  1:665949/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)

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

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

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

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

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

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

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

                                     aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                  <  1:1691149/54‑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:1695461/71‑1 (MQ=255)

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

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

                                            aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa                          >  1:249509/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)

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

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

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

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

                                                       aGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                  >  1:5504686/1‑36 (MQ=255)

                                                          aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc            >  1:5250554/1‑71 (MQ=255)

                                                          aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc            >  1:332961/1‑71 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                              <  1:3748805/36‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                              <  1:294048/36‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                              <  1:3713627/36‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                             <  1:5005324/37‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                             <  1:690733/37‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                             <  1:838098/37‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                             <  1:2055284/37‑1 (MQ=255)

                                                           tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                             <  1:4901180/37‑1 (MQ=255)

                                                            gACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                             >  1:3714040/1‑52 (MQ=255)

                                                            gACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                             >  1:2485115/1‑52 (MQ=255)

                                                             aCCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                      >  1:3707547/1‑42 (MQ=255)

                                                             aCCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                      >  1:3385383/1‑42 (MQ=255)

                                                               cAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                      >  1:1400839/1‑40 (MQ=255)

                                                               cAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                      >  1:4212884/1‑40 (MQ=255)

                                                               cAGCTGGATGCATTATCTTTATTTGGATCATTgag                                           >  1:5615586/1‑35 (MQ=255)

                                                                 gCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGc               >  1:615632/1‑61 (MQ=255)

                                                                     gATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAc  >  1:5389856/1‑70 (MQ=255)

GGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAC  >  minE/2571136‑2571274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`