BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,237,300	T→G	10.0%	L265R (CTG→CGG)	PP_2835 →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,237,300	0	T	G	10.0%	111.0 / 4.2	40	L265R (CTG→CGG)	PP_2835	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (18/18); new base G (2/2); total (20/20)
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.62e-01

CAGCATCCTGCCGCTCGGCGGCCAGGCCTTTGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGA  >  NC_002947/3237160‑3237445

cagcaTCCTGCCGCTCGGCGGCCAGGCCTTTGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTg                                                                                                                                           <  6:74778/149‑1 (MQ=255)

        tGCCGCTCGGCGGCCAGGCCTTTGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCgcgc                                                                                                                                   >  8:5903/1‑149 (MQ=255)

                   ggcCAGGCCTTTGGCCGCAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTcggcgg                                                                                                                        <  7:218763/149‑1 (MQ=255)

                         gCCTTTGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTg                                                                                                                    >  4:140457/1‑147 (MQ=255)

                              tGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCg                                                                                                              <  4:116145/148‑1 (MQ=255)

                                gCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTc                                                                                                           >  7:167794/1‑149 (MQ=255)

                                   acaAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTc                                                                                                        >  1:93912/1‑149 (MQ=255)

                                     aaaGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTcga                                                                                                      <  7:45609/149‑1 (MQ=255)

                                       aGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTcgac                                                                                                     >  2:8872/1‑148 (MQ=255)

                                         gTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCa                                                                                                                                  >  3:205497/1‑117 (MQ=255)

                                         gTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCa                                                                                                                                  <  4:205497/117‑1 (MQ=255)

                                              gCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCa                                                                                              <  7:104431/148‑1 (MQ=255)

                                              gCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAg                                                                                             >  4:88114/1‑149 (MQ=255)

                                                          ttGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTggggga                                                                                 >  2:59676/1‑148 (MQ=255)

                                                               cGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTgg                                                                            <  2:93912/149‑1 (MQ=255)

                                                                gAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGaaa                                                                                                 <  2:262461/127‑1 (MQ=255)

                                                                gAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGaaa                                                                                                 >  1:262461/1‑127 (MQ=255)

                                                                  aGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTc                                                                         <  4:192278/149‑1 (MQ=255)

                                                                         aCCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCaa                                                                   <  3:140457/148‑1 (MQ=255)

                                                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGc                                                                                                               >  1:204454/1‑102 (MQ=255)

                                                                           cTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGc                                                                                                               <  2:204454/102‑1 (MQ=255)

                                                                                ggCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCT‑GGGGGATTGGTTCTGGCAAGccgcccgc                                                          <  1:59676/149‑1 (MQ=255)

                                                                                          ccgTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCg                                                 <  5:34250/149‑1 (MQ=255)

                                                                                             tGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAAc                                              <  8:167794/149‑1 (MQ=255)

                                                                                                gcgcAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTCCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGcc                                           >  5:177256/1‑149 (MQ=255)

                                                                                                  gcAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCg                                         >  7:155290/1‑149 (MQ=255)

                                                                                                    aaGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGt                                       <  3:9800/149‑1 (MQ=255)

                                                                                                        aTGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGc                                   >  4:229655/1‑149 (MQ=255)

                                                                                                                 aTCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCggg                          <  3:229655/149‑1 (MQ=255)

                                                                                                                       ggggCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCGTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGcc                    >  4:207587/1‑148 (MQ=255)

                                                                                                                        gggCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCa                    <  1:8872/148‑1 (MQ=255)

                                                                                                                               tCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGc            >  5:244321/1‑149 (MQ=255)

                                                                                                                               tCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGc            >  2:202481/1‑149 (MQ=255)

                                                                                                                               tCGACCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGTGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGc            >  8:87368/1‑149 (MQ=255)

                                                                                                                                     catGCGCGGGTCAATCACCAGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCgg                           <  2:175673/126‑1 (MQ=255)

                                                                                                                                     catGCGCGGGTCAATCACCAGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCgg                           <  4:3584/126‑1 (MQ=255)

                                                                                                                                     catGCGCGGGTCAATCACCAGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCgg                           >  1:175673/3‑128 (MQ=255)

                                                                                                                                     catGCGCGGGTCAATCACCAGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCgg                           >  3:3584/3‑128 (MQ=255)

                                                                                                                                       tGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTc    >  2:190412/1‑149 (MQ=255)

                                                                                                                                         cgcTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGTGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGa  <  7:87368/149‑1 (MQ=255)

                                                                                                                                            tGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGc            <  6:89957/136‑1 (MQ=255)

                                                                                                                                            tGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGc            >  5:89957/1‑136 (MQ=255)

CAGCATCCTGCCGCTCGGCGGCCAGGCCTTTGGCCACAAAGGTTTTGCCCTTGCGCTGTTGGTCGAAGCCCTCACCTCGGGGCTGGCCGGCCGTGGGCGCAAGGATGCACCGGATCAGTGGGGCGCCTCGGCCACTGCGCTGGTGATTGACCCGCGCATGTTCGGCGGCCTGGAGGCGTTCCTCGACGAAACCAGCTTCCTGGGGGGATTGGTTCTGGCAAGCCGCCCGCTAGACCCCGAACGCCCGGTACGCCTGCCCGGGCAAGCAGGCCTTGCCCTGCGTCGA  >  NC_002947/3237160‑3237445

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