BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,759	A→C	31.8%	V165V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,759	0	A	C	31.8%	112.1 / 119.1	123	V165V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (57/27); new base C (18/21); total (75/48)
Fisher's exact test for biased strand distribution p-value = 2.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 127 total aligned reads displayed.

GTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTA  >  Exported/609690‑609828

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

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

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

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

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

   cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                                   <  1:2323848/71‑1 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                        aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                              >  1:830147/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)

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

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

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

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

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

                               atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                        <  1:1230781/70‑1 (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)

                                     ccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                 <  1:4890959/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:5222851/69‑1 (MQ=255)

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

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

                                        cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                              >  1:2896784/1‑71 (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:4223194/4‑70 (MQ=255)

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

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

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

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

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

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

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

                                                    tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                 <  1:5121823/56‑1 (MQ=255)

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

                                                      tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                                                 >  1:2553216/1‑38 (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:4702682/1‑60 (MQ=255)

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

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

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

                                                          aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt             >  1:1024750/3‑70 (MQ=255)

                                                           gACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                               >  1:1265131/1‑35 (MQ=255)

                                                           gACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                               <  1:1134571/35‑1 (MQ=255)

                                                            aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt          <  1:183171/71‑1 (MQ=255)

                                                              cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGcc                  <  1:870396/61‑1 (MQ=16)

                                                              cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatgg        >  1:3710304/1‑70 (MQ=255)

                                                               ggTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                 <  1:3069347/45‑1 (MQ=255)

                                                                gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                 >  1:4581481/1‑60 (MQ=12)

                                                                gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                 >  1:4575174/1‑60 (MQ=12)

                                                                gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                      >  1:1149188/1‑55 (MQ=255)

                                                                gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAac                       >  1:1333273/1‑54 (MQ=255)

                                                                 tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                               >  1:3094220/1‑45 (MQ=255)

                                                                 tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                               >  1:2025101/1‑45 (MQ=255)

                                                                  cAGAACACCAGATTTACGGCAACGCATGATGTCCGGAGTGTAAACGTCg                          <  1:417847/49‑1 (MQ=255)

                                                                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt    >  1:4540961/1‑70 (MQ=255)

                                                                    gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTTAACGTCGaa                        >  1:4909086/1‑49 (MQ=255)

                                                                    gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                        >  1:2810451/1‑49 (MQ=255)

                                                                    gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                        >  1:1354427/1‑49 (MQ=255)

                                                                     aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                      >  1:3661589/1‑50 (MQ=255)

GTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTA  >  Exported/609690‑609828

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