BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	800,934	A→G	13.1%	K59R (AAA→AGA)	ybhH →	putative isomerase YbhH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	800,934	0	A	G	13.1%	73.1 / 2.7	36	K59R (AAA→AGA)	ybhH	putative isomerase YbhH
Reads supporting (aligned to +/- strand): ref base A (16/14); new base G (3/2); total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTG  >  NC_000913/800796‑801074

ccTCGAGGGGCTCGTTCCTGTTAGCGGAACATATACACGAAGAGCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAGAGTc                                                                                                                                          >  2:191859/1‑143 (MQ=255)

  tCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGCTCTGGCAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGc                                                                                                                                        >  1:439364/1‑143 (MQ=255)

          cttgTTCCGGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCAttatta                                                                                                                                <  2:4812/140‑1 (MQ=255)

                    ttAGCGGAACGTTTACCCGAAGATCAACCGCTGTGTCGTTAGATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAg                                                                                                                      <  2:391452/143‑1 (MQ=255)

                      aGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCGGTAAAGTCGCCATTATTAGCCGGTCCAGCg                                                                                                                    >  1:384673/1‑143 (MQ=255)

                          gAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATACACTGACTGGAGGTTTGGCCATTATTAGCCGTTCCAGCGATcc                                                                                                                >  1:604767/1‑143 (MQ=255)

                            aCATTTACCCGAAGATCAAACGCATCGCGATAAAATATTGATGGCATTTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCcgc                                                                                                              <  2:374366/143‑1 (MQ=255)

                                        aGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCg                                                                                                  <  1:384179/143‑1 (MQ=255)

                                             aaaCGCAGCGCGATATAATAGTGATGGCAATTATGGGTTCCGGTAACGATCTGGAGATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTAt                                                                                             <  2:431103/143‑1 (MQ=255)

                                                      gcgATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAGAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcg                                                                                    >  2:393604/1‑143 (MQ=255)

                                                       cgATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTTACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTgcgc                                                                                   <  1:610519/143‑1 (MQ=255)

                                                             aaTATTGATGGCAAGTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGTCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTaa                                                                             <  2:120002/143‑1 (MQ=255)

                                                              atatTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAAt                                                                            >  1:152123/1‑143 (MQ=255)

                                                                         aagtctGGGTTCCGGTAACGATCTGGAAATTGACGGTATCGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGGCTATCTGCTCGCGCAGGTAATCGTTCATGAGc                                                                 >  2:858155/6‑143 (MQ=255)

                                                                               gggTTCCGGTAACGATCTGGAAATTGACTGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATCAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACgtg                                                           <  1:273066/143‑1 (MQ=255)

                                                                                 gTTTCGGTAACGTGCTGGAAATTGATGCTATTGGTGGCGGTAATCCACTCACCAGTATCGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCTTGTAATCGTTCATGAGCAACGTGTc                                                         <  2:312563/143‑1 (MQ=255)

                                                                                 gTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAGAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTc                                                         <  1:535066/143‑1 (MQ=255)

                                                                                  ttCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCg                                                        >  1:636490/1‑143 (MQ=255)

                                                                                             aTCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCACCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGCAATCGTTCATGAGCCACGTGTCGTTAGCACGACt                                             >  2:272978/1‑143 (MQ=255)

                                                                                                 ggAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACt                                         <  1:191859/143‑1 (MQ=255)

                                                                                                              aTTGGCGGCGGGAAGCCACTGACCAGTAAAGTCGCTATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCt                            <  2:333242/143‑1 (MQ=255)

                                                                                                               ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTg                           >  1:113092/1‑143 (MQ=255)

                                                                                                               ttGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGAGGTCGCTTATCTGTTTGCGCGGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAAGTGCGGCACTATGCTg                           >  2:447376/1‑143 (MQ=255)

                                                                                                                   cggcggTAATCCACTGACCAGTAGAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTg                       <  1:393604/143‑1 (MQ=255)

                                                                                                                    ggcggTAATCCACTGACCAGTAAAGTCGCCATCATTAGCCGCACCAGCGATCCGCGTGCTGATGTCGCGTATCTGCCGTCGCCGCGAAACGCGCACGAGCAACGACACGATACCACGCCTAACTGCGGCAATATGcgtacgga                      >  2:678730/1‑136 (MQ=255)

                                                                                                                           acccACCGACCAGTAAAGGCGCCATTATTAGCCGTTCCAGCGAGCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTggg               <  2:242106/141‑1 (MQ=255)

                                                                                                                               aCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCACCGTGTCGACACCACGTATAAGTGCGGCAACAGCCCGCCTGGTGTTGGGGcct           >  1:796943/1‑141 (MQ=255)

                                                                                                                               aCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATCTGCTGTCTGGCGTTGGGGCAt           >  1:817072/1‑143 (MQ=255)

                                                                                                                                 tGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTc         >  1:190457/1‑143 (MQ=255)

                                                                                                                                 tGACCAGTAAAGTCGCCATTATTAGCCGTTACGTTGATCCGCGTGCTGATGTCGATTATCTCTTCGCGCAGGTAATCGCTCAGTAGCAACGCGTCGATACCACGCCCAACCACGGCACTATGCTGTGGGGTGTGGGTGCATTc         >  2:16564/1‑143 (MQ=255)

                                                                                                                                   aCCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGc       >  1:639432/1‑143 (MQ=255)

                                                                                                                                   aCCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGc       >  1:568085/1‑143 (MQ=255)

                                                                                                                                   aCCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGCGTTGGGACATTCGc       >  1:526980/1‑143 (MQ=255)

                                                                                                                                    ccAGTAAAGTCGCCAATATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAAGCGGTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGcc      <  2:140504/143‑1 (MQ=255)

                                                                                                                                        tAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTg  <  1:376930/143‑1 (MQ=255)

                                                                                                                                        tAAAGGCGCCATTATTATCCGTCCAACCGCTCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAACCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTg  <  1:583364/143‑1 (MQ=255)

CCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTG  >  NC_000913/800796‑801074

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