BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,690,814	G→A	12.9%	L321L (CTG→CTA)	ydgA →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,690,814	0	G	A	12.9%	69.4 / 4.7	31	L321L (CTG→CTA)	ydgA	DUF945 family protein
Reads supporting (aligned to +/- strand): ref base G (11/16); new base A (2/2); total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

TAGCCAACTGGATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAGGGTGAAA  >  NC_000913/1690670‑1690959

tAGCCAACTGGATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTgg                                                                                                                                               <  2:600753/149‑1 (MQ=255)

      aCTGGATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGc                                                                                                                                         <  2:584702/149‑1 (MQ=255)

          gATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCCGa                                                                                                                                     >  1:168627/1‑149 (MQ=255)

          gATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAAAGCGCAAACTCAGGCGCTGCGGGCCCAGCCAGa                                                                                                                                     >  1:222006/1‑149 (MQ=255)

                 cGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGcc                                                                                                                              <  2:203336/149‑1 (MQ=255)

                              gAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGCCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCCGCTGGCCCAGCCAGAAACTGCCAACAACCCCGGAc                                                                                                                 >  2:1503/1‑149 (MQ=255)

                                    aCAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATc                                                                                                           >  2:540652/1‑149 (MQ=255)

                                     cAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCa                                                                                                          <  1:107393/149‑1 (MQ=255)

                                      aGAATCATGATGTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAg                                                                                                         <  1:102256/149‑1 (MQ=255)

                                           cAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATTGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGGCAACAACCCCGAACTTTATCAGGAGaa                                                                                                    <  2:301815/149‑1 (MQ=255)

                                            aGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGaaa                                                                                                   <  2:198886/149‑1 (MQ=255)

                                             ggATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGTCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAg                                                                                                  <  1:98976/149‑1 (MQ=255)

                                                      cAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGCGACGGAAg                                                                                         >  2:487139/1‑149 (MQ=255)

                                                             aaGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCcttctt                                                                                  <  1:114437/149‑1 (MQ=255)

                                                             aaGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATAAGGAGAAAGTGACGGAAGCcttctt                                                                                  <  1:347857/149‑1 (MQ=255)

                                                                   aCTTTAAAAGTCGGCCAGAATGATGATGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAgcgc                                                                            <  2:86678/149‑1 (MQ=255)

                                                                    cTTTAAAAGGCGGTCAGATTTATGGGGAAGCCTGGCGTCAGTTTAGCCAGCAATAGCACGCGCAACCTCAGGCGCTGCTGGCACAGCCAGAAATTGCCTACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGcc                                                                           <  2:222767/149‑1 (MQ=255)

                                                                     tttAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCATTCTTTAGCGccc                                                                          >  1:240892/1‑149 (MQ=255)

                                                                               ggCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGAt                                                                >  2:231715/1‑149 (MQ=255)

                                                                                 ccAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCATCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGGTGATGc                                                              >  2:139817/1‑149 (MQ=255)

                                                                                    gATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTACTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGa                                                           >  2:170942/1‑149 (MQ=255)

                                                                                    gATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTACTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGa                                                           >  2:514160/1‑149 (MQ=255)

                                                                                            gTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGAt                                                   <  1:374952/149‑1 (MQ=255)

                                                                                                aGCCCGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTCTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCgg                                               <  2:418321/149‑1 (MQ=255)

                                                                                                            gTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTACTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGTAAGGCGATCCGGTGATTACTATcg                                   <  1:514160/149‑1 (MQ=255)

                                                                                                            gTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTACTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATcg                                   <  1:170942/149‑1 (MQ=255)

                                                                                                                      cAATATAACGCGCAAACGCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGTAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAg                         <  1:424956/149‑1 (MQ=255)

                                                                                                                                       tCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTAGTATCGCGCCGCTAAGCTGGAAAAACACGCAgg        >  1:601667/1‑149 (MQ=255)

                                                                                                                                       tCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAgg        <  2:219432/149‑1 (MQ=255)

                                                                                                                                       tCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTAATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAgg        >  1:601708/1‑149 (MQ=255)

                                                                                                                                         aGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAGGGt      >  1:378950/1‑149 (MQ=255)

                                                                                                                                             gctgctGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAGGGTGaaa  <  1:540652/149‑1 (MQ=255)

TAGCCAACTGGATTACTCGCTAAACAGCCTGAAGGTACAGAATCAGGATCTGGGCAGCGGCAAGCTGACTTTAAAAGTCGGCCAGATTGATGGTGAAGCCTGGCATCAGTTTAGCCAGCAATATAACGCGCAAACTCAGGCGCTGCTGGCACAGCCAGAAATTGCCAACAACCCCGAACTTTATCAGGAGAAAGTGACGGAAGCCTTCTTTAGCGCCCTGCCGCTGATGCTGAAAGGCGATCCGGTGATTACTATCGCGCCGCTAAGCTGGAAAAACAGTCAGGGTGAAA  >  NC_000913/1690670‑1690959

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