BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,955	G→A	13.0%	F433F (TTC→TTT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,955	0	G	A	13.0%	256.9 / 24.6	92	F433F (TTC→TTT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (41/39); new base A (4/8); total (45/47)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 108 total aligned reads displayed.

TTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTG  >  minE/608887‑609024

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

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

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

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

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

  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                                    <  1:1401060/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)

       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)

              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)

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

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

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

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

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

                         tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                             >  1:657558/1‑70 (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)

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

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

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

                                  aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat                                          >  1:1450022/1‑64 (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:1040688/61‑1 (MQ=255)

                                        tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                       <  1:4123172/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)

                                             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:1199145/70‑1 (MQ=255)

                                               cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                       <  1:1670833/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)

                                                       ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                    >  1:2022989/1‑65 (MQ=255)

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

                                                       ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg              >  1:40416/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)

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

                                                           aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                        <  1:3363694/57‑1 (MQ=255)

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

                                                             gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                      >  1:282008/1‑57 (MQ=255)

                                                             gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                      >  1:2363300/1‑57 (MQ=255)

                                                             gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg        <  1:2744558/71‑1 (MQ=255)

                                                              tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa               <  1:456201/62‑1 (MQ=255)

                                                              tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATAGCGTa               <  1:802870/62‑1 (MQ=255)

                                                               gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt      >  1:3364338/1‑71 (MQ=255)

                                                               gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt      >  1:1211478/1‑71 (MQ=255)

                                                               gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct      >  1:1385780/1‑69 (MQ=255)

                                                               gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct      >  1:808507/1‑69 (MQ=255)

                                                                  ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                        <  1:2944900/50‑1 (MQ=255)

                                                                  ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt    >  1:123820/1‑70 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                     <  1:1006285/52‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                 <  1:2124285/56‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt   >  1:2056076/1‑70 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt   >  1:2061036/1‑70 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg  <  1:1725687/71‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg  <  1:151874/71‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg  <  1:2539150/71‑1 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg     >  1:528020/1‑67 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg  >  1:3036046/1‑70 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGAAGCAAGCAATAGCGTAGTCATCgttg     >  1:2618494/1‑67 (MQ=255)

TTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTG  >  minE/608887‑609024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`