BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,892	T→C	8.4%	K454K (AAA→AAG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,892	0	T	C	8.4%	297.1 / 13.1	107	K454K (AAA→AAG)	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 = 9.19e-01

Only 100 of 118 total aligned reads displayed.

CATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACT  >  Exported/608823‑608962

cATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtt                                                                       >  1:3878317/1‑71 (MQ=255)

 aTCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                      <  1:281192/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)

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

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

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

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

         atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttt                                                                      <  1:4199609/63‑1 (MQ=255)

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

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

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

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

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

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

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

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

          tcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                <  1:4471931/68‑1 (MQ=255)

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

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

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

           catAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGttttt                                                                    >  1:2376842/1‑63 (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:4941894/1‑71 (MQ=255)

                   aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                    <  1:2648773/71‑1 (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:2707379/1‑71 (MQ=255)

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

                            ccTTTGATCGGTTCAGACTTCTGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                           <  1:2881538/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                      gTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                 <  1:65963/71‑1 (MQ=255)

                                       ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                <  1:3358694/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)

                                         cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                              >  1:4130640/1‑71 (MQ=255)

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

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

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

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

                                              ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                  <  1:4029474/46‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                         <  1:808652/69‑1 (MQ=255)

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

                                                cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                       >  1:4920491/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:2942887/71‑1 (MQ=255)

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

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

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

                                                         cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg              <  1:3899536/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)

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

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

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

                                                              ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc          >  1:3640368/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)

                                                                  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg              >  1:5664726/1‑60 (MQ=255)

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

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

CATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACT  >  Exported/608823‑608962

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