BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,849,874	A→G	6.1%	G219G (GGT→GGC)	ygcI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,849,874	0	A	G	6.1%	292.3 / 7.8	98	G219G (GGT→GGC)	ygcI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (37/55); new base G (4/2); total (41/57)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGC  >  Exported/1849806‑1849941

ccTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATAc                                                                    <  1:3404345/70‑1 (MQ=255)

ccTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATAc                                                                    <  1:1660695/70‑1 (MQ=255)

  tGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATAcctc                                                                 >  1:601342/1‑71 (MQ=255)

    gTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATAcctcct                                                               >  1:2188352/1‑71 (MQ=255)

    gTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATAcctcct                                                               >  1:1791377/1‑71 (MQ=255)

     tGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCtt                                                              <  1:2047310/71‑1 (MQ=255)

         gTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAAt                                                          <  1:2861008/71‑1 (MQ=255)

            ggTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATca                                                        >  1:110693/1‑70 (MQ=255)

            ggTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATca                                                        >  1:2397685/1‑70 (MQ=255)

            ggTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATca                                                        >  1:3523342/1‑70 (MQ=255)

            ggTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATca                                                        >  1:3674264/1‑70 (MQ=255)

              tAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATc                                                         <  1:2117708/67‑1 (MQ=255)

              tAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATc                                                         <  1:3331357/67‑1 (MQ=255)

                aaGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATa                                                    <  1:918940/70‑1 (MQ=255)

                aaGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATa                                                    <  1:1170581/70‑1 (MQ=255)

                   aTCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCAt                                                >  1:2724/1‑71 (MQ=255)

                    tCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACAt                                                     <  1:2504705/65‑1 (MQ=255)

                    tCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACAt                                                     <  1:3137704/65‑1 (MQ=255)

                      ccTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtct                                             <  1:583101/71‑1 (MQ=255)

                      ccTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtct                                             >  1:3026626/1‑71 (MQ=255)

                      ccTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtc                                              <  1:1553071/70‑1 (MQ=255)

                        tGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCg                                           >  1:1316222/1‑71 (MQ=255)

                        tGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCg                                           >  1:3716063/1‑71 (MQ=255)

                          cTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATGCCTCCTTTAATCACATACCATTCTCg                                           <  1:3041990/69‑1 (MQ=255)

                          cTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATGCCTCCTTTAATCACATACCATTCTCg                                           <  1:2655425/69‑1 (MQ=255)

                          cTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAAt                                                          >  1:1762984/1‑54 (MQ=255)

                          cTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAAt                                                          >  1:584865/1‑54 (MQ=255)

                          cTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCg                                           <  1:1400969/69‑1 (MQ=255)

                             cAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGaa                                       <  1:3297206/70‑1 (MQ=255)

                             cAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGaa                                       <  1:203490/70‑1 (MQ=255)

                             cAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGaa                                       <  1:3611474/70‑1 (MQ=255)

                               ggCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATa                                                    <  1:1453889/55‑1 (MQ=255)

                                  ccTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAAt                                 <  1:274948/71‑1 (MQ=255)

                                   cTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAAtt                                <  1:559866/71‑1 (MQ=255)

                                   cTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAAtt                                <  1:2309163/71‑1 (MQ=255)

                                     ggCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtctc                                            <  1:1319710/57‑1 (MQ=255)

                                     ggCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtctc                                            <  1:594124/57‑1 (MQ=255)

                                      gCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGt                              <  1:1273185/70‑1 (MQ=255)

                                      gCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             <  1:1852512/71‑1 (MQ=255)

                                       cAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATtc                                              >  1:1281559/1‑53 (MQ=255)

                                        aaTGATGACTTTACTGAGATACATCCATGCCTCCtt                                                              >  1:2452396/1‑36 (MQ=255)

                                        aaTGATGACTTTACTGAGATACATCCATGCCTCCtt                                                              >  1:370585/1‑36 (MQ=255)

                                        aaTGATGACTTTACTGAGATACATCCATGCCTCCtt                                                              >  1:411784/1‑36 (MQ=255)

                                        aaTGATGACTTTACTGAGATACATCCATGCCTCCtt                                                              >  1:1685650/1‑36 (MQ=255)

                                        aaTGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             >  1:60544/1‑69 (MQ=255)

                                            atgaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGc                                     <  1:2259885/57‑1 (MQ=255)

                                            atgaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGc                                     <  1:3177677/57‑1 (MQ=255)

                                            atgaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGc                                     <  1:1767456/57‑1 (MQ=255)

                                            atgaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCa                       <  1:1419077/71‑1 (MQ=255)

                                             tgaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGc                        >  1:3001006/1‑69 (MQ=255)

                                              gaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCaa                                   >  1:721624/1‑57 (MQ=255)

                                              gaCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAg                     >  1:2601057/1‑71 (MQ=255)

                                               aCTTTACTGAGATACATCCATACCTCCTTTAATca                                                        <  1:2642201/35‑1 (MQ=255)

                                               aCTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAgg                    >  1:3398610/1‑71 (MQ=255)

                                                cTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCg                            >  1:2315454/1‑62 (MQ=255)

                                                cTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAg                          >  1:642467/1‑64 (MQ=255)

                                                cTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAg                          >  1:690084/1‑64 (MQ=255)

                                                cTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAg                          >  1:949913/1‑64 (MQ=255)

                                                 tttACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCTAGGCAAGGt                   >  1:267837/1‑70 (MQ=255)

                                                 tttACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGt                   >  1:197072/1‑70 (MQ=255)

                                                 tttACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGt                   >  1:3147092/1‑70 (MQ=255)

                                                   tACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAAtt                                >  1:2051730/1‑55 (MQ=255)

                                                   tACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGat                >  1:3139441/1‑71 (MQ=255)

                                                   tACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGa                 >  1:2578827/1‑70 (MQ=255)

                                                    aCTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGat                <  1:2855196/70‑1 (MQ=255)

                                                    aCTGAGATACATCCATACCTCCCTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGat                <  1:2501168/70‑1 (MQ=255)

                                                     cTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGatca              <  1:2720538/71‑1 (MQ=255)

                                                     cTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGatca              <  1:152177/71‑1 (MQ=255)

                                                     cTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGatca              <  1:796335/71‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:3543645/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:1907320/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:2026892/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:3010154/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:1282080/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:845081/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:3373898/46‑1 (MQ=255)

                                                      tGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAg                                      <  1:2296199/46‑1 (MQ=255)

                                                        agaTACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGatcat             <  1:228506/69‑1 (MQ=255)

                                                         gaTACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCgg          <  1:2117291/71‑1 (MQ=255)

                                                         gaTACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCgg          <  1:448673/71‑1 (MQ=255)

                                                            aCATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             >  1:1208928/1‑49 (MQ=255)

                                                            aCATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             >  1:2031734/1‑49 (MQ=255)

                                                            aCATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             >  1:1208521/1‑49 (MQ=255)

                                                            aCATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTc                             >  1:1848630/1‑49 (MQ=255)

                                                             cATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcg      <  1:1127674/71‑1 (MQ=255)

                                                             cATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcg      <  1:1594878/71‑1 (MQ=255)

                                                             cATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcg      <  1:1665002/71‑1 (MQ=255)

                                                             cATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcg      <  1:1901896/71‑1 (MQ=255)

                                                              aTCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcgt     >  1:605076/1‑71 (MQ=255)

                                                              aTCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcgt     <  1:1392107/71‑1 (MQ=255)

                                                              aTCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcgt     >  1:3584877/1‑71 (MQ=255)

                                                               tCCATACCTGCTTTAATCACATACCATTCTCGGGAAGCaaa                                  <  1:3538900/41‑1 (MQ=255)

                                                                ccATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcgtcg   >  1:2498556/1‑71 (MQ=255)

                                                                 cATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATTGGTTCGTcgc  <  1:2098250/71‑1 (MQ=255)

                                                                 cATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTtcgtcg   <  1:2967903/70‑1 (MQ=255)

                                                                 cATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTcgc  <  1:1389065/71‑1 (MQ=255)

                                                                 cATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTcgc  <  1:3561979/71‑1 (MQ=255)

                                                                 cATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTcgc  >  1:3386550/1‑71 (MQ=255)

CCTGGTGAAGTTGGTAAAGATCCCTGCTCCAGGCCCTGGCAATGATGACTTTACTGAGATACATCCATACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGC  >  Exported/1849806‑1849941

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