BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,220,044	A→G	S236P (TCG→CCG)	yohF ←	SDR family oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,220,044	0	A	G	88.5%	60.2 / 2.9	26	S236P (TCG→CCG)	yohF	SDR family oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (11/12); total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

CAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT  >  NZ_CP009273/2219960‑2220110

cAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCG‑CATCCACTATCAACGACTGccc                                                              <  2:342943/90‑1 (MQ=255)

 aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  1:735416/90‑1 (MQ=255)

 aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGccc                                                              <  2:410327/90‑1 (MQ=255)

  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  2:825627/89‑1 (MQ=255)

  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              >  1:825627/1‑89 (MQ=255)

  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTCTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                             >  1:96775/1‑90 (MQ=255)

            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:658715/1‑90 (MQ=255)

            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:201853/1‑90 (MQ=255)

            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:278530/1‑90 (MQ=255)

              cGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAttt                                                 >  2:500151/1‑90 (MQ=255)

                 aCCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                              <  2:540699/90‑1 (MQ=255)

                 aCCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                              <  1:22006/90‑1 (MQ=255)

                  ccTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcgc                                             <  1:125013/90‑1 (MQ=255)

                      aCTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTc                                         <  1:419369/90‑1 (MQ=255)

                           gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                    <  2:182510/90‑1 (MQ=255)

                           gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                    <  2:240506/90‑1 (MQ=255)

                                   gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                            <  1:201853/90‑1 (MQ=255)

                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              >  1:275415/1‑51 (MQ=255)

                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  2:275415/51‑1 (MQ=255)

                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  2:780305/1‑90 (MQ=255)

                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  1:265085/1‑90 (MQ=255)

                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  1:391348/1‑90 (MQ=255)

                                             gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                  >  1:701541/1‑90 (MQ=255)

                                                  ttCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa             >  1:182375/1‑90 (MQ=255)

                                                   tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa            <  1:657703/90‑1 (MQ=255)

                                                             aaaGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGt  <  1:249057/90‑1 (MQ=255)

CAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT  >  NZ_CP009273/2219960‑2220110

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