BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,159	A→C	21.5%	S699A (TCC→GCC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,159	0	A	C	21.5%	212.0 / 71.2	121	S699A (TCC→GCC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (61/34); new base C (15/11); total (76/45)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 132 total aligned reads displayed.

TTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  Exported/608089‑608229

ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                                        >  1:4004942/1‑71 (MQ=255)

    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:1071379/1‑71 (MQ=255)

    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:1309634/1‑71 (MQ=255)

    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:2694238/1‑71 (MQ=255)

         gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                                                               >  1:5125089/1‑71 (MQ=255)

         gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgt                                                                <  1:4659684/70‑1 (MQ=255)

             cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:1158104/1‑70 (MQ=255)

             cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:229330/1‑70 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:1019380/70‑1 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4259562/70‑1 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:3282078/70‑1 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2385207/70‑1 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2369404/70‑1 (MQ=255)

              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2140849/70‑1 (MQ=255)

              cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           <  1:4463507/69‑1 (MQ=255)

                gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2948280/1‑70 (MQ=12)

                gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2978923/1‑70 (MQ=12)

                gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                                                        >  1:5649424/1‑71 (MQ=255)

                 cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                                                 <  1:3371798/61‑1 (MQ=255)

                 cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                                                 <  1:1564587/61‑1 (MQ=255)

                  agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATa                                                      >  1:4622590/1‑71 (MQ=12)

                   gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATac                                                     >  1:5666752/1‑70 (MQ=12)

                    caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                    <  1:3379660/71‑1 (MQ=255)

                    caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                    >  1:1334170/1‑71 (MQ=255)

                      tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4915494/1‑71 (MQ=255)

                      tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   >  1:4274676/1‑70 (MQ=255)

                       ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2831214/1‑63 (MQ=255)

                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2085760/60‑1 (MQ=255)

                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                                                >  1:701406/1‑71 (MQ=255)

                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                                                 <  1:700995/70‑1 (MQ=255)

                          aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                                               <  1:2767935/70‑1 (MQ=255)

                           cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:4829242/71‑1 (MQ=255)

                           cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:3110125/71‑1 (MQ=255)

                           cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:5077221/71‑1 (MQ=255)

                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:3961175/1‑70 (MQ=255)

                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5401046/1‑70 (MQ=255)

                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5186034/1‑70 (MQ=255)

                             attCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:4609794/70‑1 (MQ=12)

                                cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:2194135/1‑69 (MQ=255)

                                cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:1950331/1‑69 (MQ=255)

                                catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                    <  1:1334738/59‑1 (MQ=255)

                                catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:4958019/1‑71 (MQ=255)

                                 attacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGcc                                       <  1:4386120/68‑1 (MQ=9)

                                  ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:2189330/64‑1 (MQ=9)

                                  ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:5129726/64‑1 (MQ=9)

                                  ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                      <  1:5048742/69‑2 (MQ=255)

                                  tcaCGTTAACGTTCAGGTCCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAg                                          <  1:244548/67‑1 (MQ=255)

                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                                         >  1:764585/1‑50 (MQ=255)

                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                                         >  1:3474597/1‑50 (MQ=255)

                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:3737946/1‑61 (MQ=255)

                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:1697490/1‑61 (MQ=255)

                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                     <  1:348123/70‑1 (MQ=255)

                                     cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4196643/47‑1 (MQ=255)

                                     cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                                    >  1:1497237/1‑70 (MQ=255)

                                     cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccaga                                    >  1:3046544/1‑69 (MQ=255)

                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagccag                                 >  1:5154923/1‑71 (MQ=255)

                                         gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:5356049/2‑59 (MQ=9)

                                         gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:847201/2‑59 (MQ=9)

                                         gaCGTTGAGGTGTTGACAGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:2651620/2‑59 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATcagcccagcc                                   >  1:1906388/1‑67 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcc                                       >  1:1661064/1‑63 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg                                <  1:2450950/70‑1 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:4374389/1‑71 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:2440528/1‑71 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3487062/1‑71 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3673025/1‑71 (MQ=255)

                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:5228483/1‑71 (MQ=255)

                                           cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGATCAGCCAGGTTg                             >  1:1063503/1‑71 (MQ=255)

                                            gTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTgg                            >  1:1924003/1‑71 (MQ=255)

                                            gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                            <  1:4275958/71‑1 (MQ=255)

                                              tgaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtt                          >  1:3347066/3‑70 (MQ=255)

                                               cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:4177110/1‑71 (MQ=255)

                                                  gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                       >  1:2324261/1‑70 (MQ=255)

                                                  gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      <  1:1212158/71‑1 (MQ=255)

                                                   tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4328242/1‑42 (MQ=255)

                                                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:5183963/1‑70 (MQ=255)

                                                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:2530585/1‑70 (MQ=255)

                                                     cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                   <  1:12246/71‑1 (MQ=255)

                                                      tGACCGCCTTCGACATACGCTTCGTGGTGGAAATACACATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                  >  1:1881386/1‑71 (MQ=255)

                                                      tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:1664053/1‑49 (MQ=255)

                                                      tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATAAGAc                                        >  1:3797227/1‑49 (MQ=255)

                                                       gACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     <  1:4104449/67‑1 (MQ=255)

                                                        accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:935576/1‑66 (MQ=255)

                                                        accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:3098408/1‑66 (MQ=255)

                                                            ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                    >  1:1552807/1‑63 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:4920653/1‑50 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:2331445/1‑50 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:4045519/1‑50 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5352547/1‑50 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5594752/1‑50 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:4803447/69‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:3116944/69‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt        >  1:330922/1‑71 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt        >  1:4551955/1‑71 (MQ=255)

                                                                  atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt      <  1:1041323/71‑1 (MQ=255)

                                                                     gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                           <  1:571495/47‑1 (MQ=255)

                                                                     gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   <  1:3715514/71‑1 (MQ=255)

                                                                     gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   >  1:4544237/1‑71 (MQ=255)

                                                                      atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:3626318/71‑1 (MQ=255)

                                                                      atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:260989/71‑1 (MQ=255)

TTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  Exported/608089‑608229

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