BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,423,602	T→C	R99R (CGT→CGC)	stfR →	prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,423,602	0	T	C	100.0%	57.3 / NA	17	R99R (CGT→CGC)	stfR	prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (8/9); total (8/9)

GGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTG  >  NZ_CP009273/1423523‑1423665

ggACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCAc                                                       <  2:217764/90‑1 (MQ=255)

     aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                  <  1:350440/90‑1 (MQ=255)

                 gAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTg                                      <  2:562894/90‑1 (MQ=255)

                    gATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtg                                   <  2:158681/90‑1 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                >  1:550715/1‑90 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                <  1:651563/90‑1 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                >  1:698433/1‑90 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                >  1:751906/1‑90 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                    <  1:385413/86‑1 (MQ=255)

                       tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                    >  2:385413/1‑86 (MQ=255)

                         tcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGa                              <  2:553217/90‑1 (MQ=255)

                            aCCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAaga                           <  1:92357/90‑1 (MQ=255)

                                       tGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTa                >  2:69775/1‑90 (MQ=255)

                                        gAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTaa               >  1:200188/1‑90 (MQ=255)

                                                ttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                  <  2:545454/47‑1 (MQ=255)

                                                ttCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                  >  1:545454/1‑47 (MQ=255)

                                                     ggTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTg  >  2:884793/1‑90 (MQ=255)

GGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTG  >  NZ_CP009273/1423523‑1423665

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