BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,465,502	G→A	E141K (GAG→AAG)	gspM →	Type II secretion system protein GspM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,502	0	G	A	100.0%	55.5 / NA	20	E141K (GAG→AAG)	gspM	Type II secretion system protein GspM
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (11/9); total (11/9)

CACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAGAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAGGTATTGGTTCATTTCAGGCAAAAAAAATATGCCTGGCATGATACAGTGCCCCTGATACTTTGTGTTGC  >  NC_000913/3465270‑3465727

cACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCAtt                                                                                                                                                                                                                 <  1:256139/251‑1 (MQ=255)

     aTGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGACTTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTcc                                                                                                                                                                                                            >  2:157384/1‑251 (MQ=255)

                                              aaTATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACTGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCGTGCCGGAAACGAAAAGATCACCCTCAGGCATGTGTCCTTCAGGGAGCAACAAGAATGACAATGCTACTGCCTCTCTTc                                                                                                                                                                   <  2:46766/251‑1 (MQ=255)

                                                           aGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCATCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTggttggtt                                                                                                                                                      >  2:292173/1‑251 (MQ=255)

                                                            ggAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTCACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGttt                                                                                                                                                     >  2:243533/1‑251 (MQ=255)

                                                              aaGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATAATCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATCCGGGATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGCAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGtttt                                                                                                                                                    >  1:323400/1‑250 (MQ=255)

                                                                 cAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGCAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTg                                                                                                                                                >  1:186962/1‑251 (MQ=255)

                                                                             atgtGGCTATCACTCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTAACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGACACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACACTGCTACTGCCGCTCTTCATTCTGGTTGGTTTTCTTGCAGATTATTTTg                                                                                                                                    >  2:316399/4‑251 (MQ=255)

                                                                                  gCTATCGCGCTGGAAAGCGGCCCGGAGAATTCTCTGACAATTCATCCTGTCACAGCACCCCGCGAAAACGTCTCCTGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAATATCACCGTCAGGCATTTGTCCTTCAGTGAGCAACAATTATGATAATGGTACTGCCGCTCTTCATTGTGGGGGGTTTTATTTCAGAATGTTTTGTTAAt                                                                                                                               <  2:31914/251‑1 (MQ=255)

                                                                                     aTCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGAAAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGcc                                                                                                                            >  2:75860/1‑251 (MQ=255)

                                                                                                             aaTACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTg                                                                                                    <  1:157384/251‑1 (MQ=255)

                                                                                                                                    cacaGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGtt                                                                             >  1:278219/1‑251 (MQ=255)

                                                                                                                                            cccTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAgg                                                                     <  1:75860/251‑1 (MQ=255)

                                                                                                                                              cTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTCATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAGGt                                                                    >  2:22973/1‑250 (MQ=255)

                                                                                                                                                                     cACTACGGCACAGGTAACATAAGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACGATTCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTAGTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTGAACGTTTCGCCAGGTATTGGTTCATTTCATGCaaaaaa                                             <  2:114824/250‑1 (MQ=255)

                                                                                                                                                                                aGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCa                                                                                                                                              <  2:123679/142‑1 (MQ=255)

                                                                                                                                                                                aGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCa                                                                                                                                              >  1:123679/1‑142 (MQ=255)

                                                                                                                                                                                           ggTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAGGTATTGGTTCATTTCAGGCAAAAAAAATATGCCTGGCATGATACAGTg                      <  1:243533/251‑1 (MQ=255)

                                                                                                                                                                                             gaTTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAAAGATCACCCTCATGCATTTTTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAGGTATTGGTTCATTTCAGGCAAAAAAAATATGCCTGGCATGATACAGTGcc                    <  1:219090/250‑1 (MQ=255)

                                                                                                                                                                                                               tCAGTTTACTCTTGCCGGACACGACACGCTCCACATCAGGCATTTGTCCTTCAGGGAGCCCCAATAATGACACTGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGCTTATTTTGTTAATGCCATCGCCTATCCCCTCTCGCCGTTGGCGGATAAAACGGCGTTAACGTTTCGCCAGGTATTGGTTCATTTCAGGCAAAAAAAATATGCCTGGCATGCTACAGTGCCCCTGATACTTTGTGTtgc  >  2:203474/1‑251 (MQ=255)

CACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAAAACGTCTCCCGCTGGCTCACTACGGCACAGGTAACATACGGTATTGTTATTGAAGACCTTCAGTTTACTCTTGCCGGAAACGAAGAGATCACCCTCAGGCATTTGTCCTTCAGGGAGCAACAATAATGACAATGCTACTGCCGCTCTTCATTCTGGTTGGTTTTATTGCAGATTATTTTGTTAATGCCATCGCCTATCACCTCTCGCCGTTGGAGGATAAAACGGCGTTAACGTTTCGCCAGGTATTGGTTCATTTCAGGCAAAAAAAATATGCCTGGCATGATACAGTGCCCCTGATACTTTGTGTTGC  >  NC_000913/3465270‑3465727

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