BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,747	C→T	31.1%	L169L (CTG→CTA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,747	0	C	T	31.1%	156.8 / 87.3	116	L169L (CTG→CTA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (51/28); new base T (14/22); total (66/50)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTG  >  Exported/609681‑609815

  tgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCg                                                                >  1:4832129/2‑71 (MQ=255)

   gATACCTTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                               >  1:2364034/1‑71 (MQ=255)

   gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                               >  1:4304771/1‑71 (MQ=255)

     tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                                                              >  1:1882342/1‑70 (MQ=255)

     tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                                                              >  1:652340/1‑70 (MQ=255)

       ccatACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATGTGGTAAACCGGTCAg                                                           <  1:3584754/68‑1 (MQ=255)

       ccatACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAg                                                           <  1:1652425/68‑1 (MQ=255)

       ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                                                           <  1:4061678/71‑1 (MQ=255)

       ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                                                           <  1:4882749/71‑1 (MQ=255)

       ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                                                           <  1:819147/71‑1 (MQ=255)

         gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                         >  1:1245798/1‑71 (MQ=255)

          tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                         >  1:1790933/1‑70 (MQ=255)

           aCAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCacg                                                       >  1:2452247/1‑70 (MQ=255)

            cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                      >  1:3290189/1‑71 (MQ=255)

            cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                      >  1:4111320/1‑71 (MQ=255)

            cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                      >  1:2925148/1‑71 (MQ=255)

               cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                      >  1:5112888/1‑68 (MQ=255)

               cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                      >  1:1550717/1‑68 (MQ=255)

                gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAg                                                   <  1:2540851/67‑1 (MQ=255)

                gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa                                                  <  1:3725503/68‑1 (MQ=255)

                  aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                <  1:2322149/71‑1 (MQ=255)

                   aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                >  1:343396/1‑70 (MQ=255)

                    cgcgGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACGGGTCAGCACGCCAGAtt                                                <  1:2976377/69‑1 (MQ=255)

                    cgcgGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCa                                                            >  1:2939260/1‑57 (MQ=12)

                    cgcgGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTa                                              >  1:1043154/1‑71 (MQ=255)

                    cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                              >  1:1977826/1‑71 (MQ=255)

                    cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                              <  1:3262171/71‑1 (MQ=255)

                    cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                              <  1:4904901/71‑1 (MQ=255)

                     gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                              <  1:3497576/70‑1 (MQ=255)

                     gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGCTCAGAACACCAGATTTAc                                             <  1:2408572/71‑1 (MQ=255)

                        aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                          <  1:1593060/71‑1 (MQ=255)

                        aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                          <  1:1586519/71‑1 (MQ=255)

                            tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                      >  1:3230450/1‑71 (MQ=255)

                             aGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAac                                                        >  1:4743364/1‑52 (MQ=255)

                             aGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAac                                                        >  1:4555955/1‑52 (MQ=255)

                             aGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAac                                                        >  1:2650020/1‑52 (MQ=255)

                              gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAttt                                               <  1:267034/60‑1 (MQ=255)

                              gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAttt                                               <  1:1618357/60‑1 (MQ=255)

                               tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCa                                                    <  1:5566384/54‑1 (MQ=255)

                               tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                   >  1:5050522/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACg                                            >  1:4480444/1‑60 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACg                                            >  1:5302453/1‑60 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:3452840/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:1634714/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:5506485/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:830147/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:3952179/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:5385352/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:2137121/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:1910296/1‑71 (MQ=255)

                                 aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAgg                                  <  1:438047/70‑1 (MQ=255)

                                   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                               <  1:269486/71‑1 (MQ=255)

                                    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                              >  1:4793356/1‑55 (MQ=255)

                                    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACg                                            >  1:3314651/1‑57 (MQ=255)

                                    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACg                                            >  1:2412064/1‑57 (MQ=255)

                                    gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                              >  1:840490/1‑69 (MQ=255)

                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                             >  1:3351202/1‑71 (MQ=255)

                                     atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                             >  1:417548/1‑71 (MQ=255)

                                        atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                           >  1:5603563/1‑70 (MQ=255)

                                           cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                >  1:2696991/4‑62 (MQ=255)

                                           cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                >  1:291098/4‑62 (MQ=255)

                                              ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                    <  1:119365/71‑1 (MQ=255)

                                               cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                   <  1:609488/71‑1 (MQ=255)

                                               cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                   <  1:5234997/71‑1 (MQ=255)

                                                acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                  >  1:1256356/1‑71 (MQ=255)

                                                acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                    <  1:826443/69‑1 (MQ=255)

                                                acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                    <  1:5222851/69‑1 (MQ=255)

                                                 cggccaTATGCATCTGGCAGACCTTTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                  >  1:2340349/1‑70 (MQ=255)

                                                 cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                  >  1:1699264/1‑70 (MQ=255)

                                                 cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                 >  1:5312540/1‑71 (MQ=255)

                                                  ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGAt                                                 >  1:391485/1‑38 (MQ=255)

                                                  ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                     >  1:709788/1‑50 (MQ=255)

                                                  ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                <  1:114108/71‑1 (MQ=255)

                                                  ggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAca                <  1:5066423/71‑2 (MQ=255)

                                                   gccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAgc                                  >  1:2408130/1‑51 (MQ=16)

                                                    ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATc              <  1:4805222/71‑1 (MQ=255)

                                                     caTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg             <  1:1392433/71‑1 (MQ=255)

                                                       tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                <  1:1389149/50‑1 (MQ=255)

                                                       tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                    <  1:2781513/62‑1 (MQ=255)

                                                       tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa           <  1:1112119/71‑1 (MQ=255)

                                                       tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa           <  1:4442888/71‑1 (MQ=255)

                                                        agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa           >  1:2467907/4‑70 (MQ=255)

                                                        agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa           >  1:4223194/4‑70 (MQ=255)

                                                        agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa           >  1:4424970/4‑70 (MQ=255)

                                                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                 >  1:1136120/1‑47 (MQ=255)

                                                         tGCATCTGGCAGACCGGGCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca         >  1:3822322/1‑71 (MQ=255)

                                                         gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                  >  1:35780/3‑62 (MQ=255)

                                                             tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                            <  1:2465098/48‑1 (MQ=16)

                                                             tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                            <  1:4288415/48‑1 (MQ=16)

                                                             tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggtg                        <  1:3678607/52‑3 (MQ=16)

                                                             tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                    <  1:145809/56‑1 (MQ=255)

                                                              cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc    >  1:4247352/1‑71 (MQ=255)

                                                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                    <  1:1846544/54‑1 (MQ=255)

                                                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg             >  1:3039087/1‑61 (MQ=255)

                                                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt   >  1:5075076/1‑71 (MQ=255)

                                                               tGGCAGACCGGTCAGAACACCAGATTTACGG‑AACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg  >  1:2475998/1‑71 (MQ=255)

                                                                  cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa           >  1:1376146/1‑60 (MQ=255)

                                                                  cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa           >  1:4702682/1‑60 (MQ=255)

                                                                  cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa           >  1:3868989/1‑60 (MQ=255)

                                                                  cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CGGAGTGTAAACGTCGaa           >  1:4552518/1‑59 (MQ=255)

GTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTG  >  Exported/609681‑609815

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