BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,860	G→T	8.2%	E255D (GAG→GAT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,860	0	G	T	8.2%	389.9 / 11.6	121	E255D (GAG→GAT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (54/57); new base T (4/6); total (58/63)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 126 total aligned reads displayed.

AGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCG  >  Exported/2570790‑2570929

aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                                       >  1:2878984/1‑71 (MQ=255)

 ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                                       >  1:3890918/1‑70 (MQ=255)

  gTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCa                                                                     <  1:2235105/71‑1 (MQ=255)

    aCTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCag                                                                    >  1:5457882/1‑69 (MQ=255)

     cTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCaga                                                                   <  1:1211790/70‑3 (MQ=255)

     cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                                  <  1:3691402/71‑1 (MQ=255)

     cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:3172040/70‑1 (MQ=255)

     cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:2417283/70‑1 (MQ=255)

     cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:1549848/70‑1 (MQ=255)

     cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:1625402/70‑1 (MQ=255)

       gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                                     >  1:1108039/1‑66 (MQ=255)

        aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                                 <  1:5462102/69‑1 (MQ=255)

        aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                                                <  1:4349864/70‑1 (MQ=255)

         tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACgg                                                              >  1:4744885/1‑71 (MQ=255)

           cGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                                            <  1:2631997/71‑1 (MQ=255)

             cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                                            <  1:1483664/69‑1 (MQ=255)

                 ggTACATACGCAGCTGCCTTGCCGCTGCGTACGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                                                      >  1:1046851/1‑71 (MQ=255)

                  gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                                  <  1:2733484/58‑1 (MQ=255)

                  gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTc                                                     <  1:1125899/71‑1 (MQ=255)

                    aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                                   >  1:4283257/1‑71 (MQ=255)

                       tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGc                                                                      <  1:2466572/49‑1 (MQ=255)

                       tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:5627644/52‑1 (MQ=255)

                       tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:1697807/52‑1 (MQ=255)

                        aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                                               >  1:3306124/1‑71 (MQ=255)

                        aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                                               >  1:5118750/1‑71 (MQ=255)

                         cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                              >  1:2123252/1‑71 (MQ=255)

                         cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                              >  1:4691553/1‑71 (MQ=255)

                          gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                   <  1:4283010/49‑1 (MQ=255)

                           cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            >  1:1563921/1‑71 (MQ=255)

                           cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            <  1:2527354/71‑1 (MQ=255)

                           cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            <  1:2737867/71‑1 (MQ=255)

                            aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                           <  1:447287/71‑1 (MQ=255)

                            aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                           <  1:5590133/71‑1 (MQ=255)

                              cTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTctg                                         <  1:5453895/71‑3 (MQ=255)

                                 ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCgg                                         >  1:5018281/1‑68 (MQ=255)

                                 ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                      <  1:4446137/71‑1 (MQ=255)

                                  gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAAc                                                                >  1:3872822/3‑44 (MQ=12)

                                  gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                      <  1:729192/68‑2 (MQ=255)

                                  gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                      <  1:4087924/68‑2 (MQ=255)

                                    tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            >  1:3171260/1‑62 (MQ=255)

                                    tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            >  1:1634008/1‑62 (MQ=255)

                                    tGCCGCTGCGTCCGGTGATGGTCTTCAT‑CCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                            >  1:1560527/1‑61 (MQ=255)

                                     gCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                                                      >  1:2903263/1‑51 (MQ=255)

                                       cGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTcc                                <  1:1330329/71‑2 (MQ=255)

                                        gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                                        <  1:1895200/46‑1 (MQ=255)

                                        gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGc                               <  1:774436/71‑1 (MQ=255)

                                         cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                              >  1:3564061/1‑71 (MQ=255)

                                            cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                                <  1:4251836/66‑1 (MQ=255)

                                            cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                                <  1:3251621/66‑1 (MQ=255)

                                              tCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                         <  1:3943/71‑1 (MQ=255)

                                               ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                         <  1:5338924/70‑1 (MQ=255)

                                                cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                       >  1:406426/1‑71 (MQ=255)

                                                cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                       >  1:2216233/1‑71 (MQ=255)

                                                  gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                     <  1:2653371/71‑1 (MQ=255)

                                                  gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                     <  1:5055378/71‑1 (MQ=255)

                                                   tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGc                                        <  1:2582970/51‑1 (MQ=255)

                                                   tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                     >  1:3074010/1‑70 (MQ=255)

                                                   tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTg                    <  1:5162461/71‑1 (MQ=255)

                                                     aTGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                      >  1:155038/1‑51 (MQ=255)

                                                     aTGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGa                   <  1:2546995/70‑1 (MQ=255)

                                                      tGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGaaa                 >  1:3445468/1‑71 (MQ=255)

                                                       ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc                >  1:852979/1‑71 (MQ=255)

                                                       ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc                >  1:1499124/1‑71 (MQ=255)

                                                       ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc                >  1:964811/1‑71 (MQ=255)

                                                        gTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAcc               <  1:4823359/71‑1 (MQ=255)

                                                         tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt              >  1:5240482/1‑71 (MQ=255)

                                                          cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                    >  1:211175/1‑48 (MQ=255)

                                                          cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTg             >  1:5380065/1‑71 (MQ=255)

                                                          cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTg             >  1:1113757/1‑71 (MQ=255)

                                                           ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCa                                                <  1:660700/35‑1 (MQ=255)

                                                           ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCa                                                <  1:4878030/35‑1 (MQ=255)

                                                           ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTCCGGCaa                                      <  1:733199/45‑1 (MQ=255)

                                                            tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                      >  1:3492808/1‑60 (MQ=255)

                                                            tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                      >  1:4280645/1‑60 (MQ=255)

                                                             cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                    >  1:260624/1‑45 (MQ=255)

                                                             cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                         <  1:3973375/56‑1 (MQ=255)

                                                             cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                         <  1:359928/56‑1 (MQ=255)

                                                             cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGACTGAAACCTgg            >  1:2235055/1‑69 (MQ=255)

                                                              aTCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCg          <  1:2946237/70‑1 (MQ=255)

                                                              aTCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCg          <  1:2788367/70‑1 (MQ=255)

                                                               tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                       >  1:3639720/1‑56 (MQ=255)

                                                               tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        <  1:3762874/71‑1 (MQ=255)

                                                               tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        <  1:1789805/71‑1 (MQ=255)

                                                               tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        <  1:2224427/71‑1 (MQ=255)

                                                                ccccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        >  1:5520658/1‑70 (MQ=255)

                                                                 cccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgatga      >  1:4876718/1‑71 (MQ=255)

                                                                 cccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgatg       <  1:3959626/70‑1 (MQ=255)

                                                                 cccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        >  1:1885707/1‑69 (MQ=255)

                                                                  ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                    >  1:5564383/1‑40 (MQ=255)

                                                                  ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                    >  1:3451386/1‑40 (MQ=255)

                                                                  ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                    >  1:662394/1‑40 (MQ=255)

                                                                  ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat        <  1:2358549/68‑1 (MQ=255)

                                                                  ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCCCTGTTACGTCTGAAACCTGGCgat        <  1:5389464/68‑1 (MQ=255)

                                                                   cGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                      <  1:4974498/37‑1 (MQ=255)

                                                                   cGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGc                               <  1:4690872/44‑1 (MQ=255)

                                                                    gAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGc   <  1:4336821/71‑1 (MQ=255)

                                                                     aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg  >  1:5579984/1‑71 (MQ=255)

                                                                     aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg  >  1:4182288/1‑71 (MQ=255)

                                                                     aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg  >  1:3725461/1‑71 (MQ=255)

                                                                      gCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg  >  1:401519/1‑70 (MQ=255)

AGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCG  >  Exported/2570790‑2570929

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