BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,175	G→A	18.5%	Y693Y (TAC→TAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,175	0	G	A	18.5%	274.9 / 46.0	119	Y693Y (TAC→TAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (63/34); new base A (13/9); total (76/43)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 122 total aligned reads displayed.

GCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGA  >  Exported/608105‑608245

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

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

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

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

    caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                                    <  1:3871174/71‑3 (MQ=9)

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

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

     aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                                    >  1:1072629/1‑68 (MQ=9)

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

      tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                                  >  1:2504342/1‑71 (MQ=255)

      tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                                   >  1:4274676/1‑70 (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:3082681/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                     <  1:348123/70‑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:5228483/1‑71 (MQ=255)

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

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

                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                                               >  1:3673025/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)

                                 ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                       >  1:2654879/1‑71 (MQ=255)

                                  gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                                                             <  1:5696337/48‑2 (MQ=255)

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

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

                                  gTGATGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                                                       >  1:4922195/1‑38 (MQ=255)

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

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

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

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

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

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

                                      tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                       <  1:4780395/66‑1 (MQ=255)

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

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

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

                                        accaccTACGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                                >  1:4563970/1‑71 (MQ=255)

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

                                            ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCAAGGTTGGTCTTACGa                                    >  1:3175673/1‑63 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                  atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                       >  1:1830649/1‑70 (MQ=255)

                                                   tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                      <  1:5088166/70‑1 (MQ=255)

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

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

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

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

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

                                                       tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                         >  1:3246150/1‑47 (MQ=255)

                                                       tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                         >  1:3721636/1‑47 (MQ=255)

                                                       tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                         >  1:4405462/1‑47 (MQ=255)

                                                       tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                         >  1:4465836/1‑47 (MQ=255)

                                                       tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                   >  1:64198/1‑69 (MQ=255)

                                                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                         >  1:774782/1‑62 (MQ=255)

                                                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                         >  1:1995157/1‑62 (MQ=255)

                                                         cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                            >  1:251112/1‑58 (MQ=255)

                                                          ttCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                      <  1:3583775/63‑1 (MQ=255)

                                                           tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg              >  1:4196184/1‑68 (MQ=255)

                                                            cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                     <  1:2983188/46‑1 (MQ=255)

                                                                 tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg       <  1:2945434/71‑1 (MQ=255)

                                                                      gTAACCATCCATCAGACCAGCCAGGTTGGTCTTTCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa  <  1:4596744/71‑1 (MQ=255)

GCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGA  >  Exported/608105‑608245

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