BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	457,689	C→A	22.3%	H307N (CAC→AAC)	ybgK →	predicted enzyme subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	457,689	0	C	A	22.3%	150.5 / 26.3	67	H307N (CAC→AAC)	ybgK	predicted enzyme subunit
Reads supporting (aligned to +/- strand): ref base C (30/22); new base A (12/3); total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGCTGG  >  Exported/457625‑457759

ggAAGAGGCACTGAAAGCACGGCAAGATCAACAACGTTATTTCGAACAATTAGCGTGGCGGCTGCACAAt                                                                   >  1:793842/1‑70 (MQ=255)

ggAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAAt                                                                   >  1:2958328/1‑70 (MQ=255)

 gagcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 >  1:915889‑M1/56‑71 (MQ=255)

 gAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTT‑AACAATTAGCGTGGCGGCTGCACAATGaa                                                                <  1:1114765/71‑1 (MQ=255)

  agcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 <  1:2236549‑M1/16‑1 (MQ=255)

   gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                               >  1:3710344‑M1/54‑71 (MQ=255)

   gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                               >  1:3717041‑M1/54‑71 (MQ=255)

   gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaa                                                               >  1:2977298‑M1/54‑71 (MQ=255)

   gcttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 >  1:240379‑M1/54‑69 (MQ=255)

    gagGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaaa                                                              <  1:3948132/71‑1 (MQ=255)

    cttcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGaaaa                                                              >  1:297011‑M1/53‑71 (MQ=255)

      tcagtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                             <  1:5264732‑M1/20‑1 (MQ=255)

       gCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg                                                           <  1:1744701/71‑1 (MQ=255)

       gCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg                                                           <  1:3112133/71‑1 (MQ=255)

         gtatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAatc                                                            >  1:2557310‑M1/48‑67 (MQ=255)

          tatgatcccaaagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGaaca                                                                     >  1:2224193‑M1/47‑58 (MQ=255)

          cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACaa                                                                    <  1:5033258/59‑1 (MQ=255)

          cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                                <  1:205400/63‑1 (MQ=255)

          cTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                                <  1:159483/63‑1 (MQ=255)

           tAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAAt                                                                   >  1:5669189/1‑59 (MQ=255)

           tAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGa                                                          <  1:2727195/68‑1 (MQ=255)

            aaaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaa                                                                <  1:4995752/61‑1 (MQ=255)

             aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                         >  1:1015479/1‑67 (MQ=255)

             aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCt                                                       <  1:3095273/69‑1 (MQ=255)

             aaaGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTg                                                      >  1:856055/1‑70 (MQ=255)

               aGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAAc                                                   <  1:4796028/71‑1 (MQ=255)

                gCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaaa                                                              >  1:2377816/1‑59 (MQ=255)

                   cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGcc                                                >  1:1941208/1‑70 (MQ=255)

                   cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGcc                                                >  1:3988515/1‑70 (MQ=255)

                   cGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCg                                               <  1:471770/71‑1 (MQ=255)

                    aagtttatgatgaatatcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                             >  1:1064855‑M1/37‑56 (MQ=255)

                      cAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                         >  1:4074634/1‑58 (MQ=255)

                       aaGATCAGCAACGTTATTTTGAAGAATTAGCGTGGCGGCTGCACAATGAAAATTGACCt                                                       <  1:1424377/59‑1 (MQ=255)

                       aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                        >  1:2400825/1‑58 (MQ=255)

                       aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                        >  1:2918116/1‑58 (MQ=255)

                       aaGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc                                                        >  1:170246/1‑58 (MQ=255)

                        aGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTg                                          >  1:4701182/1‑71 (MQ=255)

                         gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                                 <  1:5634265/47‑1 (MQ=255)

                         gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                                 <  1:1830136/47‑1 (MQ=255)

                         gATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGa                                                                 <  1:3226732/47‑1 (MQ=255)

                          aTCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                         >  1:1799424/1‑54 (MQ=255)

                           tCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGc                                       >  1:971858/1‑71 (MQ=255)

                            cAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGc                                       >  1:739114/1‑70 (MQ=255)

                                aaCGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                                  >  1:2989288/1‑71 (MQ=255)

                                   gTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTgg                                         <  1:3048413/61‑1 (MQ=255)

                                   atcgtaagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                             >  1:2239276‑M1/22‑41 (MQ=255)

                                    ttATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgc                              <  1:1072977/71‑1 (MQ=255)

                                    tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 >  1:5139085‑M1/21‑36 (MQ=255)

                                    tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 >  1:2228061‑M1/21‑36 (MQ=255)

                                    tcgtaagacgaagcgggcggGGCGGTTGAACAATGa                                                                 >  1:3376426‑M1/21‑36 (MQ=255)

                                      aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                               >  1:3144296/1‑36 (MQ=255)

                                      aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                               >  1:2802601/1‑36 (MQ=255)

                                      aTTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                                                               >  1:4838057/1‑36 (MQ=255)

                                      aTTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgcgc                            >  1:698919/1‑71 (MQ=255)

                                      aTTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTgcgc                            >  1:2433623/1‑71 (MQ=255)

                                       ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGa                                                          >  1:4298451/1‑40 (MQ=255)

                                       ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                         >  1:5306963/1‑41 (MQ=255)

                                       ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAc                                                         >  1:5107746/1‑41 (MQ=255)

                                       ttttGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc                                            >  1:2187942/1‑54 (MQ=255)

                                        aagacgaagcgggcggGGCGGTTGAACAATGAAAAt                                                             <  1:4967640‑M1/20‑1 (MQ=255)

                                          tGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                                  >  1:1472115/1‑61 (MQ=255)

                                          tGAACAATTAGCGTGGCGGCTCCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAgg                                  >  1:3549911/1‑61 (MQ=255)

                                               aaTTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCt                                <  1:1749936/58‑1 (MQ=255)

                                               aaTTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCt                                <  1:793860/58‑1 (MQ=255)

                                                  tAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCa                  >  1:3027193/1‑69 (MQ=255)

                                                   aGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCagag               >  1:3856260/1‑71 (MQ=255)

                                                           ggcTGCACAATGAAAATTGACCTGAACGCCGATCTGGGGGAAGGCTGCGCCAGCGACGCAGAGCTATTAAc       <  1:1462941/71‑1 (MQ=255)

                                                             cTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGc     <  1:107209/71‑1 (MQ=255)

                                                                cacaATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGCTgg  >  1:2965180/1‑71 (MQ=255)

                                                                cacaATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGCTgg  >  1:3414952/1‑71 (MQ=255)

GGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATCTGGGCGAAGGCTGCGCCAGCGACGCAGAGCTATTAACGCTGG  >  Exported/457625‑457759

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