BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,306	A→T	18.6%	V316V (GTT→GTA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,306	0	A	T	18.6%	226.0 / 45.1	102	V316V (GTT→GTA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (55/28); new base T (9/10); total (64/38)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 108 total aligned reads displayed.

ATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCT  >  minE/609238‑609375

aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                                                     >  1:4952182/1‑71 (MQ=255)

 tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                                                       <  1:2988351/68‑1 (MQ=255)

  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                                    >  1:5282517/1‑70 (MQ=255)

  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                                                     >  1:4756260/1‑69 (MQ=255)

  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                                   >  1:399609/1‑71 (MQ=255)

   gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                                   <  1:4702818/70‑1 (MQ=255)

   acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                                                                   >  1:5649674/3‑70 (MQ=255)

       cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg                                                               >  1:4007674/1‑70 (MQ=255)

            gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                                                          >  1:1803539/1‑70 (MQ=255)

             cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                                        <  1:5186223/71‑1 (MQ=255)

             cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                                        <  1:3545860/71‑1 (MQ=255)

              ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                       >  1:3512965/1‑71 (MQ=255)

               cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                      <  1:2266684/71‑1 (MQ=255)

               cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                      <  1:3130415/71‑1 (MQ=255)

                aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                      >  1:3455217/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:3794498/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:2193704/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:3589419/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:4891136/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:123839/1‑70 (MQ=255)

                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                      >  1:4897799/1‑70 (MQ=255)

                 gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                    >  1:3407185/1‑71 (MQ=255)

                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                                                                       >  1:190192/1‑52 (MQ=255)

                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                                      >  1:4145257/1‑53 (MQ=255)

                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                                      >  1:5165173/1‑53 (MQ=255)

                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                                      >  1:2798119/1‑53 (MQ=255)

                   ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                    >  1:1024587/3‑69 (MQ=255)

                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                                                       >  1:1816350/1‑50 (MQ=255)

                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                                                 >  1:3601350/1‑56 (MQ=255)

                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                  <  1:3888221/71‑1 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                                    >  1:5177806/1‑52 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                                    >  1:3362506/1‑52 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                                    >  1:5637492/1‑52 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                                    >  1:2447193/1‑52 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                         >  1:722803/1‑63 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                       >  1:5626028/1‑65 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                  <  1:5527717/70‑1 (MQ=255)

                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                  <  1:2468755/70‑1 (MQ=255)

                      ggTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                                   <  1:2814762/51‑1 (MQ=255)

                         cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                                             >  1:4023666/1‑70 (MQ=255)

                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                                   >  1:4516847/1‑63 (MQ=255)

                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                           >  1:3271539/1‑71 (MQ=255)

                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                           >  1:5303577/1‑71 (MQ=255)

                            cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                         <  1:5031927/71‑1 (MQ=255)

                            cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                                                           <  1:228596/53‑1 (MQ=255)

                            cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                         <  1:278814/71‑1 (MQ=255)

                             agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                                                >  1:1972603/1‑63 (MQ=255)

                              gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                        >  1:878281/1‑69 (MQ=255)

                              gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                        >  1:2931471/1‑69 (MQ=255)

                                  aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                   >  1:1719359/1‑71 (MQ=255)

                                  aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                   >  1:4089161/1‑71 (MQ=255)

                                   cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                   >  1:2999037/1‑70 (MQ=255)

                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                         <  1:4537443/45‑1 (MQ=255)

                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                         <  1:4923374/45‑1 (MQ=255)

                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                         <  1:64744/45‑1 (MQ=255)

                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                         <  1:2220146/45‑1 (MQ=255)

                                       tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                                            >  1:4645613/1‑41 (MQ=255)

                                       tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                                            >  1:2543439/1‑41 (MQ=255)

                                        cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                       >  1:5633186/1‑45 (MQ=255)

                                        cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                       >  1:2829690/1‑45 (MQ=255)

                                        cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                             >  1:3684893/1‑69 (MQ=255)

                                        cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                             >  1:1981534/1‑71 (MQ=255)

                                        cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                             >  1:3693472/1‑71 (MQ=255)

                                        cATCGTATTCCGGAGTAAGCAGGAAGCGAACCATACGCAGTTTCa                                                       >  1:5174622/1‑45 (MQ=255)

                                           caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                           >  1:1560493/4‑70 (MQ=255)

                                            gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                         >  1:4107251/1‑71 (MQ=255)

                                            aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                          <  1:939224/68‑2 (MQ=255)

                                              aTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                    >  1:4940931/1‑42 (MQ=255)

                                               ttCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGatg                                                  >  1:4958715/1‑41 (MQ=255)

                                               ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTcta                                        >  1:958715/1‑51 (MQ=255)

                                               ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                      >  1:5283728/1‑71 (MQ=255)

                                                tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                   <  1:4064226/57‑1 (MQ=255)

                                                 ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                                        >  1:1607554/1‑51 (MQ=255)

                                                  cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                   <  1:4578514/71‑1 (MQ=255)

                                                  cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                   <  1:236672/71‑1 (MQ=255)

                                                    gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                              >  1:1729511/3‑58 (MQ=255)

                                                    gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                         <  1:3155037/47‑1 (MQ=255)

                                                    gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                  <  1:4530147/54‑1 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:4567928/60‑3 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:1205687/60‑3 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:1301603/60‑3 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:2205378/60‑3 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:4251922/60‑3 (MQ=255)

                                                       tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                         <  1:339874/60‑3 (MQ=255)

                                                         cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                 >  1:299515/1‑50 (MQ=255)

                                                         cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                 >  1:1787213/1‑50 (MQ=255)

                                                            aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                 >  1:4362599/1‑47 (MQ=255)

                                                            aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                 >  1:2102882/1‑47 (MQ=255)

                                                            aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                 >  1:1682521/1‑47 (MQ=255)

                                                             ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                   <  1:242077/44‑1 (MQ=255)

                                                             ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                   <  1:811417/44‑1 (MQ=255)

                                                             ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg        >  1:299446/1‑71 (MQ=255)

                                                             ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg        <  1:3169799/71‑1 (MQ=255)

                                                             ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg        >  1:1523225/1‑71 (MQ=255)

                                                              gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt          >  1:4128862/1‑68 (MQ=255)

                                                              gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc       >  1:5601422/1‑71 (MQ=255)

                                                                aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg        >  1:5283690/1‑68 (MQ=255)

                                                                aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc       >  1:2070553/1‑69 (MQ=255)

                                                                   gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc           <  1:4481714/62‑1 (MQ=255)

                                                                   gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt  <  1:1970934/71‑1 (MQ=255)

ATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCT  >  minE/609238‑609375

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