BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,816	T→G	9.7%	M480L (ATG→CTG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,816	0	T	G	9.7%	248.6 / 20.9	103	M480L (ATG→CTG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (52/41); new base G (4/6); total (56/47)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 105 total aligned reads displayed.

TCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATT  >  Exported/608746‑608886

tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                                                                        >  1:4930906/1‑71 (MQ=255)

 cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                                                       >  1:5374400/1‑71 (MQ=255)

 cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                                                       >  1:5019723/1‑71 (MQ=255)

 cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                                                       >  1:5019602/1‑71 (MQ=255)

     cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcgc                                                                    <  1:35490/70‑1 (MQ=255)

     cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                                                                   >  1:5178894/1‑71 (MQ=255)

     cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                                                                   <  1:3455623/71‑1 (MQ=255)

       tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                                                  >  1:738708/1‑70 (MQ=255)

       tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                                                  >  1:3015956/1‑70 (MQ=255)

       tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                                                  >  1:4003165/1‑70 (MQ=255)

         tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                <  1:4580983/70‑1 (MQ=255)

         tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt                                                               <  1:4527083/71‑1 (MQ=255)

           gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                             <  1:2696829/71‑1 (MQ=255)

           gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                             <  1:3008225/71‑1 (MQ=255)

           gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                             >  1:2323210/1‑71 (MQ=255)

            tgaatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAc                                                            <  1:159215/68‑1 (MQ=255)

            tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                                                            <  1:2240940/71‑1 (MQ=255)

            tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                                                            <  1:954640/71‑1 (MQ=255)

              gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                                                           <  1:2768294/70‑1 (MQ=255)

              aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtt                                                          >  1:1423752/2‑71 (MQ=255)

               atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                                                         >  1:4764583/1‑70 (MQ=255)

               atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                                                         <  1:4779975/71‑2 (MQ=255)

                tgatgTTCAGTGCAGTGATGTACTGTTTACCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtca                                                        <  1:2401335/71‑1 (MQ=255)

                tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                        <  1:2684964/71‑3 (MQ=255)

                tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                        >  1:3939805/1‑69 (MQ=255)

                tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                        <  1:1182116/71‑3 (MQ=255)

                   tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                     <  1:472241/71‑1 (MQ=255)

                   tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                     <  1:1855108/71‑1 (MQ=255)

                     ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                     >  1:4852286/1‑69 (MQ=255)

                     ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                     >  1:5059907/1‑69 (MQ=255)

                     ttCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                             <  1:5063728/61‑1 (MQ=255)

                      tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat                                                       >  1:1206162/1‑66 (MQ=255)

                      tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                  >  1:1383646/1‑71 (MQ=255)

                      tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                  >  1:4597266/1‑71 (MQ=255)

                       cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                  >  1:528916/1‑70 (MQ=255)

                       cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                  >  1:5662568/1‑70 (MQ=255)

                       cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                                                 >  1:87301/1‑71 (MQ=255)

                         gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                               <  1:5513822/71‑1 (MQ=255)

                         gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                               <  1:1142408/71‑1 (MQ=255)

                         gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                               <  1:3086584/71‑1 (MQ=255)

                            cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                             <  1:5444385/70‑1 (MQ=255)

                            cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                             <  1:3585521/70‑1 (MQ=255)

                             agtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAgg                                            <  1:3987384/70‑1 (MQ=255)

                              gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                          <  1:2126794/71‑1 (MQ=255)

                              gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                          <  1:5240645/71‑1 (MQ=255)

                               tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCa                                              <  1:2086920/66‑1 (MQ=255)

                               tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                          <  1:1278528/70‑1 (MQ=255)

                               tgATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                             <  1:5043685/51‑1 (MQ=255)

                                 aTGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                             >  1:3079151/1‑65 (MQ=255)

                                    tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                <  1:4835358/43‑1 (MQ=255)

                                      cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCtt                                  >  1:295020/1‑71 (MQ=255)

                                      cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCtt                                  >  1:1035317/1‑71 (MQ=255)

                                       tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                 >  1:687099/1‑71 (MQ=255)

                                        gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                >  1:4808103/1‑39 (MQ=255)

                                        gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                >  1:1775835/1‑39 (MQ=255)

                                        gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                >  1:524488/1‑39 (MQ=255)

                                        gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                 >  1:610017/1‑70 (MQ=255)

                                           tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                             >  1:2599519/1‑39 (MQ=255)

                                           tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                             >  1:1327243/1‑39 (MQ=255)

                                           tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc                             >  1:2914802/1‑71 (MQ=255)

                                           tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc                             >  1:2921449/1‑71 (MQ=255)

                                             gccagccaGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                           <  1:772482/71‑1 (MQ=255)

                                             gccagccaGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                           <  1:4783903/71‑1 (MQ=255)

                                              ccagccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGt                          >  1:3239341/1‑71 (MQ=255)

                                               cagccagTCCATGAAGTGATCCATGCGCTACATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         <  1:1554525/71‑1 (MQ=255)

                                                 gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:5283357/1‑69 (MQ=255)

                                                 gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:3053173/1‑69 (MQ=255)

                                                 gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCa                       <  1:4025166/71‑1 (MQ=255)

                                                  ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:1093147/1‑68 (MQ=255)

                                                  ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:658166/1‑68 (MQ=255)

                                                     gTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTg                                <  1:3393067/58‑1 (MQ=255)

                                                     gTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGAc                    >  1:1429563/1‑70 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                >  1:2315182/1‑71 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                >  1:3291293/1‑71 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                >  1:5658349/1‑71 (MQ=255)

                                                          tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGAt                              <  1:5092448/55‑1 (MQ=255)

                                                          tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGa              >  1:3450673/1‑71 (MQ=255)

                                                          tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGa              >  1:2971293/1‑71 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                          >  1:4064017/1‑42 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                           >  1:3950326/1‑57 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                >  1:5571038/1‑68 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                >  1:4070583/1‑68 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGAc             >  1:621988/1‑71 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                 >  1:258179/1‑65 (MQ=255)

                                                              gTGATCCATGGGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCaa          >  1:2444895/1‑71 (MQ=255)

                                                              gTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                             <  1:4032457/36‑1 (MQ=255)

                                                              gTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                             <  1:4025868/36‑1 (MQ=255)

                                                              gTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCaa          >  1:5097586/1‑71 (MQ=255)

                                                               tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTg                                <  1:1825032/48‑1 (MQ=255)

                                                               tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGa                     >  1:1382234/1‑59 (MQ=255)

                                                               tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc         >  1:30573/1‑71 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc         <  1:5052332/68‑1 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc         <  1:2257628/68‑1 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc         <  1:1831273/68‑1 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTg      <  1:1473586/71‑1 (MQ=255)

                                                                   ccATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                 >  1:5469228/1‑59 (MQ=255)

                                                                    cATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                 <  1:669781/58‑1 (MQ=255)

                                                                    cATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                 <  1:3597882/58‑1 (MQ=255)

                                                                     aTGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCa                       <  1:4017214/51‑1 (MQ=255)

                                                                      tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt  <  1:2256830/71‑1 (MQ=255)

TCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATT  >  Exported/608746‑608886

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