BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,304	A→C	7.8%	T650T (ACT→ACG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,304	0	A	C	7.8%	158.2 / 9.5	66	T650T (ACT→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (31/29); new base C (2/3); total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCG  >  minE/608238‑608373

cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                                    <  1:778552/70‑1 (MQ=255)

 ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                  <  1:3146303/71‑1 (MQ=255)

 ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                  <  1:787453/71‑1 (MQ=255)

 ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                  <  1:1559083/71‑1 (MQ=255)

 ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                   <  1:1496844/70‑1 (MQ=255)

  gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                 <  1:2443807/71‑1 (MQ=255)

  gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                 <  1:3485917/71‑1 (MQ=255)

    aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                                 >  1:2310104/1‑67 (MQ=255)

    aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                               <  1:3595411/71‑1 (MQ=255)

    aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                               <  1:1848774/71‑1 (MQ=255)

     cGATAGAGAAGGTGTAGGAGATACCCTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                              >  1:3732737/1‑71 (MQ=255)

     cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                              >  1:2751276/1‑71 (MQ=255)

       aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                 >  1:3789780/1‑66 (MQ=255)

       aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                            >  1:1160041/1‑71 (MQ=255)

            gaAGGAGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                          >  1:116959/1‑68 (MQ=255)

             aaGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGc                                                      <  1:2397246/71‑1 (MQ=255)

              aGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                     >  1:66239/1‑71 (MQ=255)

                   tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                                >  1:1331584/1‑71 (MQ=255)

                        gaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                                    <  1:3062461/46‑1 (MQ=255)

                          tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                 <  1:3026541/47‑1 (MQ=255)

                          tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                           <  1:433874/69‑1 (MQ=255)

                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                   <  1:1565750/44‑1 (MQ=255)

                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                   <  1:2312224/44‑1 (MQ=255)

                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                        >  1:1997276/1‑71 (MQ=255)

                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                        >  1:3711575/1‑71 (MQ=255)

                           aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                        >  1:2276328/1‑71 (MQ=255)

                              aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                                                <  1:3950012/44‑1 (MQ=255)

                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCt                                                           <  1:612356/48‑1 (MQ=255)

                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                              <  1:57827/61‑1 (MQ=255)

                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                              <  1:931244/61‑1 (MQ=255)

                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                            <  1:2276715/63‑1 (MQ=255)

                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                    >  1:2824860/1‑71 (MQ=255)

                                 tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                  >  1:86378/1‑71 (MQ=255)

                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                 >  1:2871084/1‑70 (MQ=255)

                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                 >  1:2274216/1‑70 (MQ=255)

                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                 >  1:162429/1‑70 (MQ=255)

                                   tAGCGTAAGCAAACCGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                 >  1:3170732/1‑70 (MQ=255)

                                     gCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                       <  1:2989302/46‑1 (MQ=255)

                                       gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                            <  1:3083121/71‑1 (MQ=255)

                                       gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                            >  1:3206440/1‑71 (MQ=255)

                                        tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                            <  1:1532937/54‑1 (MQ=255)

                                          aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                         >  1:1485182/1‑71 (MQ=255)

                                          aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                         >  1:1218867/1‑71 (MQ=255)

                                          aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                         >  1:3070467/1‑71 (MQ=255)

                                              aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                        <  1:51290/68‑1 (MQ=255)

                                               aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                          >  1:1588358/1‑49 (MQ=255)

                                                cGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAc                   >  1:2351881/1‑71 (MQ=255)

                                                 ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                  <  1:1757453/71‑1 (MQ=255)

                                                   cAGTTTCGCCACCGAGGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCg                >  1:244307/1‑71 (MQ=255)

                                                     gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCggc              >  1:442481/1‑70 (MQ=255)

                                                      tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                      <  1:3418642/62‑1 (MQ=255)

                                                      tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                      <  1:3423367/62‑1 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                           >  1:2472862/1‑55 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc            >  1:1308062/1‑70 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc            >  1:816383/1‑70 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc            >  1:2686875/1‑70 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc            >  1:249141/1‑70 (MQ=255)

                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc            >  1:2777947/1‑70 (MQ=255)

                                                         aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                        <  1:3357585/57‑1 (MQ=255)

                                                          gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa          <  1:2898274/70‑1 (MQ=255)

                                                           cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCGCGACCGTGCATCgg                           <  1:2989599/52‑1 (MQ=255)

                                                            aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                           >  1:1076369/1‑51 (MQ=255)

                                                            aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg       <  1:1243177/71‑1 (MQ=255)

                                                            aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg       >  1:3648859/1‑71 (MQ=255)

                                                              cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACTAcc                                    >  1:3876436/1‑40 (MQ=255)

                                                                gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc     <  1:18356/69‑1 (MQ=255)

                                                                 aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg  >  1:3561611/1‑71 (MQ=255)

CGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCG  >  minE/608238‑608373

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