BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,688	G→A	11.4%	T522T (ACC→ACT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,688	0	G	A	11.4%	258.7 / 17.7	97	T522T (ACC→ACT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (46/39); new base A (4/7); total (51/46)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 103 total aligned reads displayed.

TAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGT  >  Exported/608618‑608755

tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                                                      >  1:3085328/1‑71 (MQ=255)

tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                                                      >  1:115620/1‑71 (MQ=255)

tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                                                      >  1:5075723/1‑71 (MQ=255)

 aaCTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                                                      >  1:2842431/1‑70 (MQ=255)

   cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                                                   <  1:117589/71‑1 (MQ=255)

    tttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                                                  >  1:2713393/1‑71 (MQ=255)

     ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                                                   >  1:5359604/1‑69 (MQ=255)

     ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCg                                                                 >  1:4107728/1‑71 (MQ=255)

      tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                                                  <  1:2390300/69‑1 (MQ=255)

      tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                                                  <  1:2585078/69‑1 (MQ=255)

      gggCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCga                                                                <  1:373464/69‑2 (MQ=255)

       cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                                                     >  1:1881742/1‑65 (MQ=255)

       cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                                                     >  1:1824049/1‑65 (MQ=255)

        gCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGAt                                                              <  1:937516/71‑1 (MQ=255)

         cttatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                                                             <  1:2431491/69‑1 (MQ=255)

         cATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                                                             <  1:920176/71‑1 (MQ=255)

            atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACg                                                         >  1:1149481/1‑71 (MQ=255)

             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                                        >  1:4120457/1‑71 (MQ=255)

             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                                        >  1:3464094/1‑71 (MQ=255)

             tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                                        >  1:4624953/1‑71 (MQ=255)

              ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                                        >  1:4482926/1‑70 (MQ=255)

              ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGt                                                        >  1:3549713/1‑70 (MQ=255)

               tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTc                                                       <  1:4495706/70‑1 (MQ=255)

               tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACACCATGGTGCGGATA‑ACGTCa                                                      >  1:1826472/1‑71 (MQ=255)

                  ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACgg                                                   <  1:709857/71‑1 (MQ=255)

                   tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGt                                                  >  1:3881570/1‑71 (MQ=255)

                    gCAGACAGGGAGTCAGCAGCGACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                                                 >  1:1409918/1‑71 (MQ=255)

                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa‑a                                                           >  1:5666001/1‑61 (MQ=255)

                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                                                 >  1:5365467/1‑71 (MQ=255)

                    gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                                                 >  1:4431275/1‑71 (MQ=255)

                    gCAGACAGGGAGTCAGCAGCAACGGACAGAACAGCGATACCACACGCCATGGTGCGGATa‑a                                                           >  1:898897/1‑61 (MQ=255)

                     cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                                                 <  1:70206/70‑1 (MQ=255)

                     cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCg                                                >  1:3893114/1‑71 (MQ=255)

                       gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                              <  1:2664318/71‑1 (MQ=255)

                        aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                              >  1:4380746/1‑70 (MQ=255)

                        aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGc                                             <  1:259841/71‑1 (MQ=255)

                         cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGc                                             <  1:5063168/70‑1 (MQ=255)

                         cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                            <  1:2870297/71‑1 (MQ=255)

                          aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAg                                           >  1:4929491/1‑71 (MQ=255)

                            ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                            >  1:4325930/1‑68 (MQ=255)

                            ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                            >  1:767944/1‑68 (MQ=255)

                                tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:5488704/71‑1 (MQ=255)

                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAg                                           >  1:4966893/1‑63 (MQ=255)

                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     >  1:3859894/1‑69 (MQ=255)

                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     >  1:2339797/1‑69 (MQ=255)

                                  agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                   <  1:5218487/71‑1 (MQ=255)

                                       caACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑Ac                                                          <  1:1283199/43‑1 (MQ=255)

                                        aaCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑Ac                                                          <  1:1612146/42‑1 (MQ=255)

                                         aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagaga                               >  1:5597367/1‑68 (MQ=255)

                                         aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagaga                               >  1:606469/1‑68 (MQ=255)

                                         aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGc                            >  1:225641/1‑71 (MQ=255)

                                          ccgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                           >  1:2536424/3‑71 (MQ=255)

                                           ggACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                          <  1:203083/71‑1 (MQ=255)

                                             aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                        >  1:190803/1‑71 (MQ=255)

                                                   ccgCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTg                  >  1:2590246/3‑71 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                          >  1:4804471/1‑63 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                          >  1:789345/1‑63 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                          >  1:449662/1‑63 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                          >  1:4797998/1‑63 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                  <  1:4899071/71‑1 (MQ=255)

                                                   cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGGCGTGCAGCGCCAt                                     <  1:4551539/52‑1 (MQ=255)

                                                     gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:4739138/50‑1 (MQ=255)

                                                     gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:5505039/50‑1 (MQ=255)

                                                     gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:774347/50‑1 (MQ=255)

                                                     gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:1076884/50‑1 (MQ=255)

                                                     gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                >  1:1757715/1‑71 (MQ=255)

                                                      cGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                     <  1:3946373/49‑1 (MQ=255)

                                                       gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                              >  1:5614179/1‑39 (MQ=255)

                                                       gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                 >  1:2844041/1‑52 (MQ=255)

                                                       gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                   >  1:761863/1‑50 (MQ=255)

                                                       gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTAc               >  1:4218048/1‑70 (MQ=255)

                                                       gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACt              >  1:4981778/1‑71 (MQ=255)

                                                        aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                         >  1:2933757/1‑59 (MQ=255)

                                                         tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAgc                                          >  1:350219/1‑41 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                >  1:1411553/1‑49 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                >  1:2049538/1‑49 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                  <  1:4109552/63‑1 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                  <  1:1775790/63‑1 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt           <  1:4063157/70‑1 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt           <  1:1880320/70‑1 (MQ=255)

                                                           ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc          <  1:4524482/71‑1 (MQ=255)

                                                            cacaCGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt           <  1:3882643/69‑1 (MQ=255)

                                                             acacGCCATGGTGCTGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt        <  1:5588762/71‑1 (MQ=255)

                                                             acacGCCATGGTGCGGATG‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt        <  1:5425087/71‑1 (MQ=255)

                                                             acacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                >  1:5594075/1‑47 (MQ=255)

                                                             acacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                >  1:2945116/1‑47 (MQ=255)

                                                              catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg       <  1:3557553/68‑1 (MQ=255)

                                                              catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg       <  1:4083297/68‑1 (MQ=255)

                                                              catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg       <  1:3102444/68‑1 (MQ=255)

                                                              cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg       >  1:3314966/1‑71 (MQ=255)

                                                               acGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGc      <  1:3358513/71‑1 (MQ=255)

                                                                tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATc                                    <  1:965006/39‑1 (MQ=9)

                                                                tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATc                                    <  1:126045/39‑1 (MQ=9)

                                                                tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATc                                    <  1:1138975/39‑1 (MQ=9)

                                                                tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa     >  1:568748/2‑71 (MQ=255)

                                                                tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa     >  1:3838314/2‑71 (MQ=255)

                                                                  ccATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg   <  1:4580196/71‑1 (MQ=255)

                                                                   cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt  <  1:4603465/71‑1 (MQ=255)

                                                                    aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg       <  1:5161873/65‑1 (MQ=255)

                                                                    aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg       <  1:3682641/65‑1 (MQ=255)

TAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGT  >  Exported/608618‑608755

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