BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,212,852	C→G	5.5%	I4M (ATC→ATG)	yobH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,212,852	0	C	G	5.5%	293.7 / 4.8	90	I4M (ATC→ATG)	yobH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (34/51); new base G (3/2); total (37/53)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 102 total aligned reads displayed.

TTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGTTTT  >  Exported/1212783‑1212917

ttgcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATtcatca                                                                  >  1:2718569/1‑71 (MQ=255)

ttgcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATtcatca                                                                  >  1:1682448/1‑71 (MQ=255)

 tgcttCAGCGCATATTGACGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                 >  1:3156554/1‑71 (MQ=255)

  gcttCAGCGCATATTGGCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                 <  1:3277791/70‑1 (MQ=255)

  gcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATtcatc                                                                   <  1:882500/68‑1 (MQ=255)

  gcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                 <  1:2525969/70‑1 (MQ=255)

  gcttCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                 <  1:2389046/70‑1 (MQ=255)

      cAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                            >  1:523720/1‑71 (MQ=255)

      cAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                            >  1:1701447/1‑71 (MQ=255)

       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGcac                                                           <  1:2620396/71‑1 (MQ=255)

       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                            >  1:3447645/1‑70 (MQ=255)

       aGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGca                                                            >  1:1169397/1‑70 (MQ=255)

       aGCGCATATTGCCGCTAAGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGAATGATGAATCGcat                                                           <  1:1251357/71‑2 (MQ=255)

        gcgcATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGAAATGCGATTCATCATTCgaa                                                            >  1:678781/1‑67 (MQ=255)

          gcATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATtcatca                                                                  <  1:3315856/61‑1 (MQ=255)

             tatTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                     <  1:1491297/71‑1 (MQ=255)

              atTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                     <  1:3568419/70‑1 (MQ=255)

                tGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                  >  1:3044916/1‑71 (MQ=255)

                 gCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                  <  1:429803/70‑1 (MQ=255)

                 gCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                  <  1:478818/70‑1 (MQ=255)

                   cGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                               <  1:896662/71‑1 (MQ=255)

                   cGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                               <  1:443524/71‑1 (MQ=255)

                    gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                  <  1:1584533/67‑1 (MQ=255)

                    gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGAt                                               <  1:1345729/70‑1 (MQ=255)

                        cGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCAc                                          <  1:713195/71‑1 (MQ=255)

                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:3144332/1‑70 (MQ=255)

                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:2579033/1‑70 (MQ=255)

                            tAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:2463037/1‑70 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTcatcat                                                                 <  1:18404/43‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:1521084/70‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:126628/70‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:2026097/70‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:3463506/70‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:3490180/70‑1 (MQ=255)

                             aaGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:171884/70‑1 (MQ=255)

                             aaGCGAACATAAAAACAGAAGAGGTTGTAATGCGATGCATCATTCGCACAGTTATGCTGATAGCACTGGt                                      <  1:2440365/70‑1 (MQ=255)

                              aGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTa                                     <  1:681767/70‑1 (MQ=255)

                                  aaCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGAt                                <  1:375387/71‑1 (MQ=255)

                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                     >  1:1311326/1‑49 (MQ=255)

                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTAt                                                     >  1:934189/1‑49 (MQ=255)

                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                 >  1:2195428/1‑69 (MQ=255)

                                   aCATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                 >  1:888625/1‑69 (MQ=255)

                                    cATAAAAAGAGAAGATGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAg                              <  1:2618588/71‑1 (MQ=255)

                                    cATAAAAAGAGAAGAGGTTGTAATGCGATtcatca                                                                  >  1:1632999/1‑35 (MQ=255)

                                    cATAAAAAGAGAAGAGGTTGTAATGCGATtcatca                                                                  >  1:1832148/1‑35 (MQ=255)

                                    cATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAg                              >  1:705242/1‑71 (MQ=255)

                                      tAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGa                            <  1:571527/71‑1 (MQ=255)

                                       aaaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                           >  1:2184047/1‑71 (MQ=255)

                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                           >  1:528763/1‑70 (MQ=255)

                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                           >  1:40880/1‑70 (MQ=255)

                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                           >  1:111163/1‑70 (MQ=255)

                                        aaaaGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGat                           >  1:1843003/1‑70 (MQ=255)

                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                        >  1:2751613/1‑71 (MQ=255)

                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                        <  1:1842050/71‑1 (MQ=255)

                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattat                        <  1:2496561/71‑1 (MQ=255)

                                          aagagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTAATAGCACTGGTATGGATAGtattat                        <  1:52290/71‑1 (MQ=255)

                                          aagagaagagGTTGTAATGCGAATGATGAATCGcat                                                           >  1:2090826/1‑35 (MQ=9)

                                          aagagaagagGTTGTAATGCGAATGATGAATCGcat                                                           >  1:559331/1‑35 (MQ=9)

                                          aagagaagagGTTGTAATGCGAATGATGAATCGca                                                            >  1:2551730/1‑35 (MQ=12)

                                           agagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGattatt                       <  1:659146/71‑1 (MQ=255)

                                            gagaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGAttatta                      <  1:2558826/71‑1 (MQ=255)

                                             agaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTAc                     >  1:2550229/1‑71 (MQ=255)

                                              gaagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCt                                                  <  1:978167/41‑1 (MQ=255)

                                               aagagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACt                    <  1:3134944/70‑1 (MQ=255)

                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                        >  1:791758/1‑65 (MQ=255)

                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                        >  1:1774031/1‑65 (MQ=255)

                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTTGTATGGATAGGattat                        >  1:1383582/1‑65 (MQ=255)

                                                agagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACt                    >  1:3594056/1‑69 (MQ=255)

                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                          <  1:106768/62‑1 (MQ=255)

                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                          <  1:2101532/62‑1 (MQ=255)

                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGatt                          <  1:2105977/62‑1 (MQ=255)

                                                 gagGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAg                 <  1:728624/71‑1 (MQ=255)

                                                   ggTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCg               <  1:2855628/71‑1 (MQ=255)

                                                       gTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGata           >  1:907658/1‑71 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATTGATAGGATTATTACTCAGCGGatat          <  1:2131645/71‑1 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                 <  1:1879091/48‑1 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGa                                 <  1:3235325/48‑1 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat          <  1:3319188/71‑1 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat          <  1:1886657/71‑1 (MQ=255)

                                                        tAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat          <  1:1933167/71‑1 (MQ=255)

                                                         aaTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGatat          <  1:87869/70‑1 (MQ=255)

                                                          gtGCGAATGATGAATCGCACAGTTATGCTGATAGCACTGGTATg                                   <  1:3482005/43‑1 (MQ=18)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:618864/1‑40 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:2079543/1‑40 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:2024988/1‑40 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:809168/1‑40 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       >  1:1356725/1‑40 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTg                                        >  1:2711318/1‑39 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg         <  1:2201186/70‑1 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg         <  1:936376/70‑1 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg         <  1:418561/70‑1 (MQ=255)

                                                          aTGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATg         <  1:1840887/70‑1 (MQ=255)

                                                           tGCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       <  1:935766/39‑1 (MQ=255)

                                                            gCGATTCATCATTCGCACAGTTATGCTGATAGCACTgg                                       <  1:3285753/38‑1 (MQ=255)

                                                            gCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCg      >  1:526194/1‑71 (MQ=255)

                                                             cGATTGATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGt     <  1:1054421/71‑1 (MQ=255)

                                                               aTTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGttt   <  1:692002/71‑1 (MQ=255)

                                                               aTTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGttt   <  1:2265151/71‑1 (MQ=255)

                                                                ttCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGtttt  >  1:2748453/1‑71 (MQ=255)

                                                                    tcatTCGCACAGTTATGCTGATAGCACTGTTATGGATAGGATTATTACTCAGCGGat            <  1:3572123/57‑1 (MQ=255)

TTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTGTAATGCGATTCATCATTCGCACAGTTATGCTGATAGCACTGGTATGGATAGGATTATTACTCAGCGGATATGGCGTTTT  >  Exported/1212783‑1212917

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