BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,021	T→G	8.8%	K745Q (AAA→CAA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,021	0	T	G	8.8%	226.2 / 14.2	91	K745Q (AAA→CAA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (53/30); new base G (3/5); total (56/35)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGG  >  Exported/607953‑608087

gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAg                                                                  <  1:4346096/71‑1 (MQ=255)

  tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTc                                                                <  1:355191/71‑1 (MQ=255)

      aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGc                                                             <  1:2413541/70‑1 (MQ=255)

      aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGc                                                             <  1:2762980/70‑1 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        >  1:973397/1‑71 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        >  1:229397/1‑71 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        >  1:4352840/1‑71 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        >  1:2591745/1‑71 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        <  1:2151974/71‑1 (MQ=255)

          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                        <  1:645509/71‑1 (MQ=255)

             tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGa                                                           <  1:1309261/65‑1 (MQ=255)

              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                    >  1:4477486/1‑71 (MQ=255)

              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                    >  1:2170650/1‑71 (MQ=255)

              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                    >  1:982748/1‑71 (MQ=255)

                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                    <  1:4419698/68‑1 (MQ=255)

                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  <  1:1962975/70‑1 (MQ=255)

                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  >  1:2298041/1‑70 (MQ=255)

                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                 <  1:3026259/71‑1 (MQ=255)

                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                 >  1:4396549/1‑70 (MQ=255)

                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                >  1:2300988/1‑71 (MQ=255)

                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                 >  1:3494523/1‑70 (MQ=255)

                     cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACt                                              <  1:4815956/70‑1 (MQ=255)

                      aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTa                                                                   <  1:756690/48‑1 (MQ=255)

                      aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTa                                                                   <  1:3725039/48‑1 (MQ=255)

                      aTAGATTGAGTGAAGGTACGAGTAATAAAGTCCTGCTGCT‑GTTCTTTa                                                                   <  1:4201361/48‑1 (MQ=255)

                       tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGGCAGCGAGTTGAAACGTACTGc                                            <  1:3749529/70‑1 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:2338224/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:176854/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:71531/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:3197352/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:210154/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:2244981/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                         >  1:2561674/1‑56 (MQ=255)

                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                           >  1:3536577/1‑70 (MQ=255)

                           ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                       <  1:2617571/71‑1 (MQ=255)

                            tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  >  1:1801633/1‑59 (MQ=255)

                            tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  >  1:1953076/1‑59 (MQ=255)

                             gAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                                      <  1:83408/70‑1 (MQ=255)

                               gtaaaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                                   <  1:350587/68‑1 (MQ=255)

                                tGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                 <  1:1294776/56‑1 (MQ=255)

                                  aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                 <  1:404675/70‑3 (MQ=255)

                                  aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                 <  1:2969174/70‑3 (MQ=255)

                                  aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                 >  1:5315911/1‑70 (MQ=255)

                                  aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                 >  1:227011/1‑70 (MQ=255)

                                   aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  >  1:5010636/1‑52 (MQ=255)

                                   aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                  >  1:805418/1‑52 (MQ=255)

                                      tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                              >  1:3279435/1‑69 (MQ=255)

                                      tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                             >  1:398507/1‑70 (MQ=255)

                                          aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                        <  1:87811/71‑1 (MQ=255)

                                          aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                        <  1:3814349/71‑1 (MQ=255)

                                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                            <  1:3326983/65‑1 (MQ=255)

                                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                            <  1:5158368/65‑1 (MQ=255)

                                            taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                             <  1:366473/64‑1 (MQ=255)

                                             aataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                     <  1:569736/71‑1 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                          >  1:3593189/1‑65 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                          >  1:531529/1‑65 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                          >  1:5538430/1‑65 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                          >  1:2715948/1‑65 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                    <  1:472346/71‑1 (MQ=255)

                                              ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                    >  1:1064665/1‑71 (MQ=255)

                                              ataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                    >  1:4040379/1‑71 (MQ=255)

                                               taaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                     <  1:3450180/69‑1 (MQ=255)

                                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                >  1:2849456/1‑41 (MQ=255)

                                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                >  1:739180/1‑41 (MQ=255)

                                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                >  1:3235411/1‑41 (MQ=255)

                                                aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                        <  1:1961296/49‑1 (MQ=255)

                                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                          >  1:1115134/1‑46 (MQ=255)

                                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                            >  1:2975706/1‑60 (MQ=255)

                                                 aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                     >  1:2340301/1‑67 (MQ=255)

                                                   gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg               >  1:3563115/1‑71 (MQ=255)

                                                   gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg               >  1:5527442/1‑71 (MQ=255)

                                                    tCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                       <  1:5225139/62‑1 (MQ=255)

                                                     ccTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTAc                                               >  1:2986568/1‑37 (MQ=255)

                                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                     >  1:753132/1‑62 (MQ=255)

                                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa             >  1:3375138/1‑70 (MQ=255)

                                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa             >  1:1178530/1‑70 (MQ=255)

                                                      ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa             >  1:2966344/1‑70 (MQ=255)

                                                      ctgctgca‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa             >  1:562899/1‑70 (MQ=255)

                                                       tgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                  <  1:1232043/64‑1 (MQ=255)

                                                         ctgctgGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                            >  1:4867945/1‑37 (MQ=38)

                                                            ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                      <  1:1052531/39‑1 (MQ=255)

                                                            ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt      <  1:1119205/71‑1 (MQ=255)

                                                                ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc    <  1:1488244/70‑1 (MQ=255)

                                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg  >  1:2144533/1‑71 (MQ=255)

                                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg  >  1:5385737/1‑71 (MQ=255)

                                                                 tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg  >  1:3291109/1‑71 (MQ=255)

                                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                   >  1:5045808/1‑37 (MQ=255)

                                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                   >  1:4501310/1‑37 (MQ=255)

                                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                 >  1:2000940/1‑39 (MQ=255)

                                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                 >  1:2851405/1‑39 (MQ=255)

                                                                  cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                    >  1:4205704/1‑52 (MQ=255)

GATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGG  >  Exported/607953‑608087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`