BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,895	T→C	8.4%	E453E (GAA→GAG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,895	0	T	C	8.4%	297.5 / 13.1	107	E453E (GAA→GAG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (47/51); new base C (3/6); total (50/57)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 124 total aligned reads displayed.

TCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTG  >  Exported/608825‑608963

tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                      <  1:235755/71‑1 (MQ=255)

  aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                     <  1:1096068/70‑1 (MQ=255)

  aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                     <  1:2272679/70‑1 (MQ=255)

  aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                    >  1:4312569/1‑71 (MQ=255)

   cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                   <  1:1614455/71‑1 (MQ=255)

   cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                   <  1:912948/71‑1 (MQ=255)

   cTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                   <  1:569113/71‑1 (MQ=255)

    ttCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                      <  1:859462/67‑1 (MQ=255)

      catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                 <  1:948639/70‑1 (MQ=255)

      catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                 >  1:4100724/1‑70 (MQ=255)

      catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                 >  1:3276520/1‑70 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                               >  1:4876915/1‑71 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                               >  1:5548118/1‑71 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                               >  1:1828819/1‑71 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                <  1:4223156/70‑1 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                <  1:2193148/70‑1 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                <  1:3371711/70‑1 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                <  1:1591369/70‑1 (MQ=255)

       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                <  1:4369909/70‑1 (MQ=255)

        tcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                              <  1:1901824/71‑1 (MQ=255)

        tcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                              <  1:1910813/71‑1 (MQ=255)

        tcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                              <  1:706396/71‑1 (MQ=255)

         catAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                     >  1:2376842/1‑63 (MQ=255)

         catAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                 <  1:3962185/67‑1 (MQ=255)

          atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                            >  1:4099931/1‑71 (MQ=255)

                 aGGATATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                     >  1:4215052/1‑71 (MQ=255)

                 aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCg                                                                   <  1:4789398/57‑1 (MQ=255)

                 aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                     <  1:2648773/71‑1 (MQ=255)

                 aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                     >  1:4941894/1‑71 (MQ=255)

                 aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                     >  1:808107/1‑71 (MQ=255)

                  ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                    <  1:3891587/71‑1 (MQ=255)

                  ggACATAGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                     >  1:2870897/1‑70 (MQ=255)

                    aCATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                  >  1:1881430/1‑71 (MQ=255)

                      aTCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                   <  1:2356034/68‑1 (MQ=255)

                       tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                               <  1:2574888/71‑1 (MQ=255)

                       tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                               <  1:4092283/71‑1 (MQ=255)

                        cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                                  >  1:3859678/1‑51 (MQ=255)

                        cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                              >  1:3047628/1‑71 (MQ=255)

                         gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                               >  1:1206350/1‑53 (MQ=255)

                          ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                  >  1:4518497/1‑65 (MQ=255)

                          ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                            <  1:4843710/71‑1 (MQ=255)

                          ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                            <  1:4510550/71‑1 (MQ=255)

                          ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                            <  1:2868979/71‑1 (MQ=255)

                            tttGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                            >  1:1517222/1‑69 (MQ=255)

                              tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                 >  1:4690974/1‑46 (MQ=255)

                               gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                   >  1:2231155/1‑59 (MQ=255)

                                   ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAAGCCGCCGTTGATTGCGTACAGCATGGt                                   <  1:3284944/71‑1 (MQ=255)

                                     ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                 <  1:4326348/71‑1 (MQ=255)

                                       cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                <  1:2553028/70‑1 (MQ=255)

                                       cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                >  1:5268579/1‑70 (MQ=255)

                                         gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                             >  1:1264034/1‑71 (MQ=255)

                                          aCTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                            <  1:1821597/71‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                            >  1:545681/1‑37 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                       <  1:1190732/58‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                       <  1:4165177/58‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                       <  1:2527513/58‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                       <  1:3575453/58‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                     <  1:1342268/60‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                     <  1:4109264/60‑1 (MQ=255)

                                            ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                  >  1:2264402/1‑63 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                     >  1:3448639/1‑58 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                     >  1:4370299/1‑58 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                             >  1:2451882/1‑66 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                         <  1:1832260/70‑1 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                         <  1:3010404/70‑1 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                        >  1:4920491/1‑71 (MQ=255)

                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                        >  1:3961435/1‑71 (MQ=255)

                                                   cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcg                   >  1:5607572/3‑70 (MQ=255)

                                                   cAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                   >  1:2086234/1‑71 (MQ=255)

                                                    gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                  >  1:4024147/2‑69 (MQ=255)

                                                     aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                            >  1:2503781/1‑60 (MQ=255)

                                                     aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                 <  1:1813025/71‑1 (MQ=255)

                                                     aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                 <  1:2942887/71‑1 (MQ=255)

                                                     aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                 <  1:4610772/71‑1 (MQ=255)

                                                     aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                 <  1:62444/71‑1 (MQ=255)

                                                       cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg               >  1:2418760/1‑71 (MQ=255)

                                                       cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg               <  1:604767/71‑3 (MQ=255)

                                                         gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                            >  1:2643916/1‑40 (MQ=255)

                                                         gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa             >  1:2833316/1‑71 (MQ=255)

                                                         gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa             >  1:763935/1‑71 (MQ=255)

                                                          cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc            >  1:664257/1‑71 (MQ=255)

                                                            ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc           >  1:3300909/1‑70 (MQ=255)

                                                              ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                         <  1:3972192/38‑1 (MQ=255)

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

                                                               tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                         <  1:5693198/37‑1 (MQ=255)

                                                               tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg          >  1:4919861/1‑68 (MQ=255)

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

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

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

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

                                                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa     <  1:2485728/71‑1 (MQ=255)

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

                                                                   tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt   <  1:3820837/71‑1 (MQ=255)

                                                                    ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                >  1:899933/1‑41 (MQ=255)

                                                                     tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                          >  1:42513/1‑46 (MQ=255)

                                                                     tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                       >  1:3815029/1‑49 (MQ=255)

                                                                     tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                       >  1:2630020/1‑49 (MQ=255)

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

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

                                                                     tttCGTCAACGCCACCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                       >  1:1673985/1‑49 (MQ=255)

TCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTG  >  Exported/608825‑608963

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