BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,190	A→G	7.9%	G688G (GGT→GGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,190	0	A	G	7.9%	279.6 / 11.2	101	G688G (GGT→GGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (58/35); new base G (4/4); total (62/39)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 104 total aligned reads displayed.

CATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGG  >  Exported/608121‑608260

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

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

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

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

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

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

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

    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                                                    <  1:1231279/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)

         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:2440528/1‑71 (MQ=255)

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

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

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

                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                     <  1:4439662/69‑1 (MQ=255)

                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                                    >  1:1493020/1‑70 (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:3098408/1‑66 (MQ=255)

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

                        accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                <  1:1572301/70‑1 (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:4045519/1‑50 (MQ=255)

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

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

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

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

                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                                             >  1:4920653/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)

                                  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:756104/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                         cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                              >  1:5692686/1‑71 (MQ=255)

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

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

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

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

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

                                                gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                       >  1:5161950/1‑71 (MQ=255)

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

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

                                                         aCCATCTATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa               <  1:5236490/70‑1 (MQ=255)

                                                         aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa                 >  1:1746875/1‑68 (MQ=255)

                                                          ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga             <  1:2638319/71‑1 (MQ=255)

                                                          ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga             <  1:304087/71‑1 (MQ=255)

                                                           catccatcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                          <  1:3559112/41‑1 (MQ=255)

                                                           catccatcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                          <  1:208932/41‑1 (MQ=255)

                                                           catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag            <  1:837849/71‑1 (MQ=255)

                                                           catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag            <  1:5529521/71‑1 (MQ=255)

                                                            atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga           <  1:3290948/71‑1 (MQ=255)

                                                             tccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                       >  1:1827767/1‑58 (MQ=255)

                                                              ccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                      >  1:3014024/1‑42 (MQ=255)

                                                               catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                >  1:2290821/1‑47 (MQ=255)

                                                                 tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      <  1:151381/71‑1 (MQ=255)

                                                                  aaGACCAACCTGGCTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgt       >  1:1873971/2‑69 (MQ=38)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGaa                    >  1:5671284/1‑55 (MQ=255)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga             >  1:676437/1‑62 (MQ=255)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga           >  1:1938276/1‑64 (MQ=255)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAg         >  1:362555/1‑66 (MQ=255)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      >  1:1957170/1‑69 (MQ=255)

                                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      >  1:106044/1‑69 (MQ=255)

                                                                     aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga             <  1:5558185/60‑1 (MQ=255)

                                                                     aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg  <  1:2844479/71‑1 (MQ=255)

CATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGG  >  Exported/608121‑608260

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