BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,946	A→G	16.3%	R436R (CGT→CGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,946	0	A	G	16.3%	235.6 / 38.0	98	R436R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (44/38); new base G (5/11); total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 103 total aligned reads displayed.

AACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAT  >  minE/608877‑609016

aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                       <  1:1016534/71‑1 (MQ=255)

aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                       <  1:3040962/71‑1 (MQ=255)

aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                       <  1:2910320/71‑1 (MQ=255)

aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                       <  1:2681248/71‑1 (MQ=255)

 aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                      >  1:1329070/1‑71 (MQ=255)

  ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                                                     >  1:2043102/1‑71 (MQ=255)

    tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                                    >  1:24167/1‑70 (MQ=255)

      cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACgcgc                                                                   <  1:1895898/69‑1 (MQ=255)

          ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                               >  1:3654931/1‑69 (MQ=255)

          ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                               >  1:2381497/1‑69 (MQ=255)

          ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                               >  1:1494808/1‑68 (MQ=255)

            cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                            <  1:102482/70‑1 (MQ=255)

            cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                            <  1:27440/70‑1 (MQ=255)

            cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                                              <  1:1916967/68‑3 (MQ=255)

            cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                                              <  1:3865116/68‑3 (MQ=255)

            cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                               <  1:3618961/67‑2 (MQ=255)

            cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                            <  1:1401060/70‑1 (MQ=255)

            cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                            <  1:1703505/70‑1 (MQ=255)

             aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                          >  1:595649/1‑71 (MQ=255)

             aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                                           <  1:2728383/70‑1 (MQ=255)

             aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                                           <  1:3160660/70‑1 (MQ=255)

              gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                         <  1:2532648/71‑1 (MQ=255)

               tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt                                                        >  1:2008479/1‑71 (MQ=255)

                ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                       >  1:3491701/1‑71 (MQ=255)

                ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                       >  1:3852655/1‑71 (MQ=255)

                 tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                      >  1:1668574/1‑55 (MQ=255)

                 tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                       >  1:3339379/1‑70 (MQ=255)

                  ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                            >  1:229808/1‑64 (MQ=255)

                   tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                     <  1:2971132/70‑1 (MQ=255)

                    cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGATCTGCAt                                                    <  1:1270/70‑1 (MQ=255)

                     gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                    <  1:52892/69‑1 (MQ=255)

                      tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                                                 <  1:2163221/71‑1 (MQ=255)

                        aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                               >  1:1856739/1‑71 (MQ=255)

                         aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                              <  1:3815239/71‑1 (MQ=255)

                          ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                             <  1:2023686/69‑1 (MQ=255)

                          ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                             <  1:725150/69‑1 (MQ=255)

                           gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                        >  1:1648793/1‑43 (MQ=255)

                           gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                        >  1:2535825/1‑43 (MQ=255)

                           gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                                            >  1:2708355/1‑71 (MQ=255)

                            ccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                                                             >  1:2585727/1‑53 (MQ=255)

                              gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                                         >  1:128065/1‑71 (MQ=255)

                                cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                        >  1:3075564/1‑70 (MQ=255)

                                cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                       <  1:677044/71‑1 (MQ=255)

                                   tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                        >  1:2935016/1‑35 (MQ=255)

                                   tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                     >  1:657558/1‑70 (MQ=255)

                                    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                       >  1:2819634/1‑35 (MQ=255)

                                    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                   >  1:587614/1‑71 (MQ=255)

                                    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                   >  1:1852293/1‑71 (MQ=255)

                                    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                    <  1:4005748/70‑1 (MQ=255)

                                    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca                                   >  1:685199/1‑71 (MQ=255)

                                      ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                 <  1:1823578/71‑1 (MQ=255)

                                      ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                 <  1:910646/71‑1 (MQ=255)

                                        gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                                          >  1:3238564/1‑44 (MQ=12)

                                        gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                                  >  1:2117140/1‑36 (MQ=255)

                                            aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                  >  1:111947/1‑64 (MQ=255)

                                            aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat                                  >  1:1450022/1‑64 (MQ=255)

                                               gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                            >  1:1288974/1‑67 (MQ=255)

                                               gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                        <  1:43335/71‑1 (MQ=255)

                                               gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                        <  1:2648335/71‑1 (MQ=255)

                                               gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                        >  1:3004311/1‑71 (MQ=255)

                                               gCATGGGTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAttt                                                  >  1:538428/1‑45 (MQ=255)

                                                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                          <  1:1360880/50‑1 (MQ=255)

                                                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                               <  1:2636176/61‑1 (MQ=255)

                                                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                               <  1:4123172/61‑1 (MQ=255)

                                                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                               <  1:1040688/61‑1 (MQ=255)

                                                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg                        >  1:3760010/1‑68 (MQ=255)

                                                  gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                     >  1:1398399/3‑71 (MQ=255)

                                                   ggTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                     >  1:119737/1‑38 (MQ=255)

                                                    gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                   <  1:3498025/71‑1 (MQ=255)

                                                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGc                                                      <  1:408750/35‑1 (MQ=255)

                                                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                              <  1:3826659/43‑1 (MQ=255)

                                                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                      <  1:2918905/67‑1 (MQ=255)

                                                     ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                         <  1:3893971/48‑1 (MQ=25)

                                                       ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                    >  1:3826947/1‑35 (MQ=255)

                                                       ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                    >  1:3624840/1‑51 (MQ=255)

                                                       ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                    <  1:1990208/51‑1 (MQ=255)

                                                        tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                <  1:3395187/70‑1 (MQ=255)

                                                         cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca               <  1:11636/70‑1 (MQ=255)

                                                         cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca               <  1:1670833/70‑1 (MQ=255)

                                                         cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca               <  1:1199145/70‑1 (MQ=255)

                                                         cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca               <  1:1900468/70‑1 (MQ=255)

                                                            cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg           <  1:1472220/71‑1 (MQ=255)

                                                            cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGctt                                              <  1:2241829/36‑1 (MQ=16)

                                                             aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc          >  1:4119712/1‑71 (MQ=255)

                                                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg           >  1:4078426/1‑69 (MQ=255)

                                                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg           >  1:2729084/1‑69 (MQ=255)

                                                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg         <  1:2076321/71‑1 (MQ=255)

                                                              ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg         >  1:1221440/1‑71 (MQ=255)

                                                               gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGt        <  1:1555201/71‑1 (MQ=255)

                                                                tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa       <  1:232203/71‑1 (MQ=255)

                                                                 ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa              >  1:1117979/1‑63 (MQ=255)

                                                                 ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg      >  1:40416/1‑71 (MQ=255)

                                                                 ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg      >  1:593652/1‑71 (MQ=255)

                                                                  tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg      >  1:488209/1‑70 (MQ=255)

                                                                  tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt     >  1:509491/1‑71 (MQ=255)

                                                                  tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt     >  1:2004154/1‑71 (MQ=255)

                                                                  tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt     >  1:2525102/1‑71 (MQ=255)

                                                                  tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt     >  1:3324975/1‑71 (MQ=255)

                                                                    cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                           <  1:1916882/47‑1 (MQ=255)

                                                                     aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAt  <  1:3381657/71‑1 (MQ=255)

AACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAT  >  minE/608877‑609016

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