BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,278	C→A	6.3%	E326* (GAA→TAA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,278	0	C	A	6.3%	305.0 / 4.1	94	E326* (GAA→TAA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (55/33); new base A (4/2); total (59/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 101 total aligned reads displayed.

GGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGC  >  minE/609210‑609348

ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                                      <  1:1293650/71‑1 (MQ=255)

  tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                                    <  1:3889893/71‑1 (MQ=255)

  tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                                    <  1:2079290/71‑1 (MQ=255)

      aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                                                >  1:1698968/1‑71 (MQ=255)

      aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                                                                 >  1:624818/1‑70 (MQ=255)

       ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                               >  1:4124624/1‑71 (MQ=255)

       ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                               >  1:572776/1‑71 (MQ=255)

         gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                             <  1:850734/71‑1 (MQ=255)

         gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                             <  1:1691527/71‑1 (MQ=255)

          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                             >  1:4205076/1‑70 (MQ=255)

             aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                                          >  1:4514107/1‑70 (MQ=255)

                 ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                                                      >  1:3484677/1‑70 (MQ=255)

                   aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                   >  1:4629968/1‑71 (MQ=255)

                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                  <  1:1712763/71‑1 (MQ=255)

                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                  <  1:3820930/71‑1 (MQ=255)

                      ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                                 <  1:3361070/70‑1 (MQ=255)

                      ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                                 <  1:821195/70‑1 (MQ=255)

                        accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                             >  1:378396/1‑56 (MQ=255)

                         ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                               >  1:5616054/1‑53 (MQ=255)

                          cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                             <  1:1657633/70‑1 (MQ=255)

                           gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                                                     >  1:503729/1‑45 (MQ=255)

                            aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                                        <  1:5490046/41‑1 (MQ=255)

                            aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                          >  1:4952182/1‑71 (MQ=255)

                             tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                               <  1:4114301/49‑1 (MQ=255)

                             tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                            <  1:2988351/68‑1 (MQ=255)

                              aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                              >  1:1460356/1‑65 (MQ=255)

                              aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                         >  1:5282517/1‑70 (MQ=255)

                              aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                          >  1:4756260/1‑69 (MQ=255)

                              aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                        >  1:399609/1‑71 (MQ=255)

                               gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                        <  1:4702818/70‑1 (MQ=255)

                               acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                                        >  1:5649674/3‑70 (MQ=255)

                                   cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg                                    >  1:4007674/1‑70 (MQ=255)

                                        gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                               >  1:1803539/1‑70 (MQ=255)

                                         cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                             <  1:5186223/71‑1 (MQ=255)

                                         cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                             <  1:3545860/71‑1 (MQ=255)

                                          ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                            >  1:3512965/1‑71 (MQ=255)

                                           cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                           <  1:2266684/71‑1 (MQ=255)

                                           cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                           <  1:3130415/71‑1 (MQ=255)

                                           cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                             <  1:3613448/53‑1 (MQ=255)

                                           cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                             <  1:42153/53‑1 (MQ=255)

                                           cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                             <  1:4897229/53‑1 (MQ=255)

                                            aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                           >  1:3455217/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:4891136/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:3794498/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:123839/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:3589419/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:2193704/1‑70 (MQ=255)

                                            aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                           >  1:4897799/1‑70 (MQ=255)

                                             gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                         >  1:3407185/1‑71 (MQ=255)

                                             gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                                            >  1:190192/1‑52 (MQ=255)

                                             gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGc                                                     <  1:2008814/43‑1 (MQ=255)

                                             gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                           >  1:5165173/1‑53 (MQ=255)

                                             gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                           >  1:2798119/1‑53 (MQ=255)

                                             gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                           >  1:4145257/1‑53 (MQ=255)

                                              aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                           >  1:5447180/1‑36 (MQ=255)

                                              aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                           >  1:440374/1‑36 (MQ=255)

                                              aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                        <  1:3607582/71‑1 (MQ=255)

                                               ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                         >  1:1024587/3‑69 (MQ=255)

                                               tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                                                    >  1:352507/1‑42 (MQ=255)

                                               tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                            >  1:1816350/1‑50 (MQ=255)

                                               tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                      >  1:3601350/1‑56 (MQ=255)

                                               tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                      >  1:23699/1‑56 (MQ=255)

                                               tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                       <  1:3888221/71‑1 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                         >  1:5177806/1‑52 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                         >  1:5637492/1‑52 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                         >  1:3362506/1‑52 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                         >  1:2447193/1‑52 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                              >  1:722803/1‑63 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                            >  1:5626028/1‑65 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                       <  1:5527717/70‑1 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                       <  1:2468755/70‑1 (MQ=255)

                                                cGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCg                                             >  1:5453611/1‑48 (MQ=255)

                                                  ggTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                        <  1:2814762/51‑1 (MQ=255)

                                                     cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                  >  1:4023666/1‑70 (MQ=255)

                                                      gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                             >  1:3567643/1‑42 (MQ=255)

                                                      gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                        >  1:4516847/1‑63 (MQ=255)

                                                      gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                >  1:3271539/1‑71 (MQ=255)

                                                      gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                >  1:5303577/1‑71 (MQ=255)

                                                        cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa              <  1:5031927/71‑1 (MQ=255)

                                                        cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                                <  1:228596/53‑1 (MQ=255)

                                                        cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa              <  1:278814/71‑1 (MQ=255)

                                                         agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                     >  1:1972603/1‑63 (MQ=255)

                                                          gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                      <  1:4348556/45‑1 (MQ=255)

                                                          gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat             >  1:2931471/1‑69 (MQ=255)

                                                            gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa          <  1:1885619/71‑1 (MQ=255)

                                                              aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt        >  1:1719359/1‑71 (MQ=255)

                                                              aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt        >  1:4089161/1‑71 (MQ=255)

                                                               cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt        >  1:2999037/1‑70 (MQ=255)

                                                                  ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                              <  1:2220146/45‑1 (MQ=255)

                                                                  ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                              <  1:4923374/45‑1 (MQ=255)

                                                                  ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                              <  1:64744/45‑1 (MQ=255)

                                                                  ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                              <  1:4537443/45‑1 (MQ=255)

                                                                   tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                 >  1:2543439/1‑41 (MQ=255)

                                                                   tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                 >  1:4645613/1‑41 (MQ=255)

                                                                    cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                            >  1:5633186/1‑45 (MQ=255)

                                                                    cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                            >  1:2829690/1‑45 (MQ=255)

                                                                    cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac  >  1:3684893/1‑69 (MQ=255)

                                                                    cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc  >  1:1981534/1‑71 (MQ=255)

                                                                    cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc  >  1:3693472/1‑71 (MQ=255)

                                                                    cATCGTATTCCGGAGTAAGCAGGAAGCGAACCATACGCAGTTTCa                            >  1:5174622/1‑45 (MQ=255)

GGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGC  >  minE/609210‑609348

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