BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,316	T→C	7.7%	V646V (GTA→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,316	0	T	C	7.7%	179.1 / 5.4	65	V646V (GTA→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (27/33); new base C (2/3); total (29/36)
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.12e-01

AGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGAC  >  minE/608249‑608386

agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACt                                                                     >  1:2131970/1‑71 (MQ=255)

  aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                                    >  1:1123783/1‑70 (MQ=255)

       gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                                              >  1:953217/1‑71 (MQ=255)

        tacgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCttt                                                             <  1:1753908/68‑1 (MQ=255)

          ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                                            <  1:2401660/70‑1 (MQ=255)

              aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                                       <  1:1654394/71‑1 (MQ=255)

                aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                                       >  1:1331851/1‑69 (MQ=255)

                  cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGa                                                   >  1:720036/1‑71 (MQ=255)

                   aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACATTTCTGGTCACGAc                                                  >  1:3260619/1‑71 (MQ=255)

                    tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAc                                                  >  1:1051499/1‑70 (MQ=255)

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

                        tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGTACCTTTCTGGTCACGACCGTg                                              >  1:87529/1‑70 (MQ=255)

                        tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                                 <  1:2037412/67‑1 (MQ=255)

                           cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                                          <  1:3250035/71‑1 (MQ=255)

                           cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTAACGACCGTGc                                             <  1:913851/68‑1 (MQ=255)

                           cGTAAGAAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                                          <  1:3164314/71‑1 (MQ=255)

                             tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                                            <  1:3559758/51‑1 (MQ=255)

                             tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                        >  1:1019022/1‑71 (MQ=255)

                             tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAATg                                                                    <  1:355408/43‑1 (MQ=255)

                               aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                                   >  1:2251034/1‑42 (MQ=255)

                               aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCATCTTTCTGGTCACGACCGTGCATCGGGt                                      <  1:3696793/71‑1 (MQ=255)

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

                               aGCAAACGGCAGTTTAGCAACGGAAGGCAGCGCGGCTACTGCACCTTTCTGGTCACGACCGTGCATCggg                                       <  1:3344432/70‑1 (MQ=255)

                                gtgaaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGtt                                     <  1:2517230/68‑1 (MQ=255)

                                gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                        <  1:1104843/52‑1 (MQ=255)

                                gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                        >  1:1086037/1‑68 (MQ=255)

                                gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                     <  1:3646045/71‑1 (MQ=255)

                                   aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAc                                                                      >  1:1135848/1‑35 (MQ=255)

                                   aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                                  >  1:1794823/1‑39 (MQ=255)

                                   aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                                          <  1:1799456/63‑1 (MQ=255)

                                    aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCa                                 <  1:552843/71‑1 (MQ=255)

                                     cGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcgc                                >  1:2072388/1‑69 (MQ=255)

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

                                         aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                            <  1:539325/71‑1 (MQ=255)

                                         aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                            <  1:1867519/71‑1 (MQ=255)

                                          gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCCTGCATCGGGTTAGCGCCCggc                           <  1:3263540/71‑2 (MQ=255)

                                           tttAGCATCGGAAGTCTGAGAGGCTACTGCACCTTTc                                                            <  1:2397951/37‑1 (MQ=255)

                                           tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                                 <  1:2470539/48‑1 (MQ=255)

                                           tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc                          >  1:3031956/1‑71 (MQ=255)

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

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

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

                                             tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc                          >  1:2297856/1‑69 (MQ=255)

                                              aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                               >  1:323299/1‑63 (MQ=255)

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

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

                                                  aCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACgg                   >  1:3186160/1‑71 (MQ=255)

                                                  aCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg                    <  1:3366397/70‑1 (MQ=255)

                                                   cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                  >  1:1922/1‑71 (MQ=255)

                                                    ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcg                 >  1:3165444/1‑71 (MQ=255)

                                                      aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                                  <  1:751043/52‑1 (MQ=255)

                                                      aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                                  <  1:1505466/52‑1 (MQ=255)

                                                       aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc              >  1:779293/1‑71 (MQ=255)

                                                        gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCg                             <  1:2161746/55‑1 (MQ=255)

                                                           agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                     >  1:546803/1‑44 (MQ=255)

                                                           agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCa            <  1:1961746/69‑1 (MQ=255)

                                                           agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc          <  1:2220336/71‑1 (MQ=255)

                                                             agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac        <  1:2923959/71‑1 (MQ=255)

                                                              gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                     <  1:506275/41‑1 (MQ=255)

                                                              gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                      <  1:3905117/56‑1 (MQ=255)

                                                              gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac        >  1:696155/1‑70 (MQ=255)

                                                              gagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac        >  1:1172352/1‑70 (MQ=255)

                                                               agGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                      <  1:2689728/39‑1 (MQ=255)

                                                                ggCTACTGCACCTTTCTGGACACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacgac     <  1:3200847/71‑1 (MQ=255)

                                                                 gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac        <  1:836964/67‑1 (MQ=255)

                                                                   tACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac  >  1:3728350/1‑71 (MQ=255)

AGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGAC  >  minE/608249‑608386

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