BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,943	T→A	17.0%	A437A (GCA→GCT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,943	0	T	A	17.0%	237.1 / 39.2	94	A437A (GCA→GCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (38/40); new base A (5/11); total (43/51)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 104 total aligned reads displayed.

ACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGT  >  minE/608874‑609013

aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                                       <  1:1038274/71‑1 (MQ=255)

aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                                       <  1:2792016/71‑1 (MQ=255)

aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                                       <  1:2491041/71‑1 (MQ=255)

aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                                       <  1:1971130/71‑1 (MQ=255)

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

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

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

   aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                    <  1:3040962/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:2381497/1‑69 (MQ=255)

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

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

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

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

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

               cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                                           <  1:1916967/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:3852655/1‑71 (MQ=255)

                   ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                    >  1:3491701/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:725150/69‑1 (MQ=255)

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

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

                              gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                     >  1:1648793/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: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:2648335/71‑1 (MQ=255)

                                                  gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                     <  1:43335/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:1040688/61‑1 (MQ=255)

                                                     tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                            <  1:4123172/61‑1 (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:1990208/51‑1 (MQ=255)

                                                          ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                 >  1:3624840/1‑51 (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:2729084/1‑69 (MQ=255)

                                                                 ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg        >  1:4078426/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)

                                                                    ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa         >  1:2022989/1‑65 (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)

ACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGT  >  minE/608874‑609013

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