BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,121	C→T	7.7%	G342G (GGC→GGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,121	0	C	T	7.7%	356.7 / 13.2	117	G342G (GGC→GGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (66/42); new base T (6/3); total (72/45)
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.42e-01

Only 100 of 121 total aligned reads displayed.

GGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAA  >  minE/2571053‑2571191

ggCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCAt                                                                      >  1:3283887/1‑71 (MQ=255)

 gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:920613/1‑71 (MQ=255)

 gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:778516/1‑71 (MQ=255)

 gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     <  1:484090/71‑1 (MQ=255)

 gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                     >  1:1084205/1‑71 (MQ=255)

      gacggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGc                                                                >  1:1474017/1‑71 (MQ=255)

       acggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGcc                                                               >  1:2809320/1‑71 (MQ=255)

        cggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGccc                                                              <  1:2890056/71‑1 (MQ=255)

             tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTg                                                         >  1:1394229/1‑71 (MQ=255)

             ggCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTg                                                         <  1:971848/70‑1 (MQ=255)

               cTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTgg                                                        <  1:2814524/70‑1 (MQ=255)

                tGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGcc                                                               >  1:413449/1‑62 (MQ=255)

                tGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCg                                                      >  1:1924455/1‑71 (MQ=255)

                 gggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     >  1:1207196/1‑71 (MQ=255)

                  ggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:3121236/1‑57 (MQ=255)

                   gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   <  1:3886982/71‑1 (MQ=255)

                    aaCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:1040565/1‑70 (MQ=255)

                     aCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTgg                                                        <  1:800904/64‑1 (MQ=255)

                      ccATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                                                        >  1:2862358/1‑47 (MQ=255)

                      ccATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    <  1:215910/51‑1 (MQ=255)

                       cATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                               <  1:1072581/71‑1 (MQ=255)

                        aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:1949305/1‑51 (MQ=255)

                        aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:831212/1‑51 (MQ=255)

                        aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              <  1:1564331/71‑1 (MQ=255)

                        aTTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              <  1:2012768/71‑1 (MQ=255)

                         ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  <  1:2137484/50‑1 (MQ=255)

                         ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:2924826/1‑65 (MQ=255)

                         ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                             >  1:2782638/1‑71 (MQ=255)

                         ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              <  1:2857808/70‑1 (MQ=255)

                         ttGAACTGGACCATAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                   >  1:2257528/1‑65 (MQ=255)

                          tGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGc                                                       >  1:657469/1‑60 (MQ=255)

                           gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTAc                                                                  >  1:2853522/1‑48 (MQ=255)

                            aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCa                                                             <  1:3041690/52‑1 (MQ=255)

                            aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                     >  1:3641952/1‑60 (MQ=255)

                             aCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAcc                                         <  1:1789316/71‑1 (MQ=255)

                                 gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:1893495/1‑71 (MQ=255)

                                 gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAt                                       <  1:181142/69‑1 (MQ=255)

                                 gACCCTAACGATGATATCGATGCCTTAACCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGaa                                          >  1:1334145/1‑65 (MQ=255)

                                   cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:2810103/1‑69 (MQ=255)

                                   cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                     >  1:2571257/1‑69 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    >  1:3737068/1‑37 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    >  1:1573982/1‑37 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGcc                                                               >  1:1864418/1‑42 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGttt                                                  <  1:2134404/55‑1 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                  >  1:3908733/1‑71 (MQ=255)

                                    ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                  >  1:1721624/1‑71 (MQ=255)

                                      tAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                    >  1:2301272/1‑35 (MQ=255)

                                      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:1274267/1‑57 (MQ=255)

                                      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCg                                            <  1:593424/59‑1 (MQ=255)

                                      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                 >  1:2467528/1‑70 (MQ=255)

                                      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                 >  1:3405967/1‑70 (MQ=255)

                                       aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               >  1:2523851/1‑71 (MQ=255)

                                       aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               >  1:3696805/1‑71 (MQ=255)

                                        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                           >  1:3162418/1‑42 (MQ=255)

                                        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                           >  1:1069144/1‑42 (MQ=255)

                                        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                              >  1:2296564/1‑71 (MQ=255)

                                         cGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCa                                   >  1:828070/1‑65 (MQ=255)

                                           atgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                              >  1:613427/1‑68 (MQ=255)

                                            tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                          <  1:1883885/71‑1 (MQ=255)

                                             gatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                          >  1:1149055/1‑70 (MQ=255)

                                              atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                        >  1:1771985/1‑71 (MQ=255)

                                                atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:3848879/62‑1 (MQ=255)

                                                atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:2926085/62‑1 (MQ=255)

                                                atCGATGCCTTACACGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      >  1:2142560/1‑71 (MQ=255)

                                                 tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:504025/61‑1 (MQ=255)

                                                 tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      <  1:1753831/70‑1 (MQ=255)

                                                  cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:3213189/1‑48 (MQ=255)

                                                  cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATg                             >  1:3206653/1‑62 (MQ=255)

                                                  cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    >  1:3917161/1‑71 (MQ=255)

                                                   gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    <  1:1608881/70‑1 (MQ=255)

                                                   gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                    <  1:2318931/70‑1 (MQ=255)

                                                   gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                   <  1:1329041/71‑1 (MQ=255)

                                                    aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGcaa                               <  1:3331794/58‑3 (MQ=255)

                                                    aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                               <  1:3485451/58‑1 (MQ=255)

                                                     tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:1718118/1‑42 (MQ=255)

                                                     tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:353665/1‑45 (MQ=255)

                                                     tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGa                                           >  1:194779/1‑45 (MQ=255)

                                                        tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:1826468/2‑36 (MQ=25)

                                                        tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:1066394/2‑36 (MQ=25)

                                                        tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTa                                                 >  1:122854/2‑36 (MQ=25)

                                                        tttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgc               >  1:2862143/2‑69 (MQ=255)

                                                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTc                                              >  1:3900741/1‑39 (MQ=255)

                                                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                      <  1:2898723/63‑1 (MQ=255)

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

                                                         ttACCCGACGGCATCTACGCCCACGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa              >  1:3533788/1‑70 (MQ=255)

                                                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa              >  1:3348979/1‑69 (MQ=255)

                                                          tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca              <  1:3022781/69‑3 (MQ=255)

                                                           aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa             <  1:2904217/69‑1 (MQ=255)

                                                           aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg           >  1:21810/1‑71 (MQ=255)

                                                            cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGTTGAAAg           <  1:508836/70‑1 (MQ=255)

                                                            cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg           <  1:3333110/70‑1 (MQ=255)

                                                               gACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc        >  1:1792418/1‑70 (MQ=255)

                                                               gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct       >  1:1197432/1‑71 (MQ=255)

                                                                aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                              >  1:2463630/1‑47 (MQ=255)

                                                                 cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat     <  1:2438973/71‑3 (MQ=255)

                                                                 cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctgc     >  1:652245/1‑71 (MQ=255)

                                                                   gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct    <  1:203597/70‑1 (MQ=255)

                                                                   gCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGCtctgct    <  1:30981/70‑1 (MQ=255)

                                                                    cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct    >  1:2786236/1‑69 (MQ=255)

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

GGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAA  >  minE/2571053‑2571191

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