BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,043	(C)_6→7	coding (372/1035 nt)	yehA ←	putative fimbrial‑like adhesin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,188,037	1	.	C	100.0%	111.1 / NA	33	G126G (GGC→GGG)	yehA	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (19/14); total (19/14)

CACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAA  >  NC_000913/2187916‑2188158

cacaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACAt                                                                                                           >  1:1256160/1‑139 (MQ=255)

 acaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACAtt                                                                                                          <  1:72439/139‑1 (MQ=255)

                   gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAAccgccg                                                                                                                            >  1:741870/1‑103 (MQ=255)

                   gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAAccgccg                                                                                                                            <  2:741870/103‑1 (MQ=255)

                   gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGc                                                                                        >  1:1230958/1‑139 (MQ=255)

                   gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑cccc                                                                                                                       >  2:713272/1‑107 (MQ=255)

                   gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑cccc                                                                                                                       <  1:713272/107‑1 (MQ=255)

                                  caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACCCCATTAATCATCACTCAACCGCCGCCCCCCCGCAGAAACATTCAc                                                                                                       <  2:747916/109‑1 (MQ=255)

                                  caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCAc                                                                                                       >  1:747916/1‑109 (MQ=255)

                                  caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGt                                                                         >  1:845823/1‑139 (MQ=255)

                                   aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCgccg‑cc                                                                                                                         <  1:1162833/89‑1 (MQ=255)

                                   aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCgccg‑cc                                                                                                                         >  2:1162833/1‑89 (MQ=255)

                                     tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc                                                                            >  1:1123400/1‑133 (MQ=255)

                                     tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc                                                                            <  2:1123400/133‑1 (MQ=255)

                                             gTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa                                                                        >  1:364410/1‑129 (MQ=255)

                                             gTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa                                                                        <  2:364410/129‑1 (MQ=255)

                                               aCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCa                                                            >  2:241995/1‑139 (MQ=255)

                                               aCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAAGCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCa                                                            >  1:970811/1‑139 (MQ=255)

                                                   cGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCACAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTgg                                                        >  1:343381/1‑139 (MQ=255)

                                                        cTCCCGCGTGTCGATTATGAAGCACTGGATACTGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGaaa                                                   <  2:874352/139‑1 (MQ=255)

                                                        cTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGaaa                                                   >  1:872276/1‑139 (MQ=255)

                                                             gcgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGt                                              <  2:872276/139‑1 (MQ=255)

                                                             gcgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCACAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGt                                              <  2:343381/139‑1 (MQ=255)

                                                              cgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTc                                             >  2:1025787/1‑139 (MQ=255)

                                                                                  ggATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAt                         >  2:1294243/1‑139 (MQ=255)

                                                                                       aTGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAAc                                                                      <  2:993358/89‑1 (MQ=255)

                                                                                       aTGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAAc                                                                      >  1:993358/1‑89 (MQ=255)

                                                                                        tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCt                                            >  2:185278/1‑114 (MQ=255)

                                                                                        tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCt                                            <  1:185278/114‑1 (MQ=255)

                                                                                        tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACg                             >  2:701348/1‑129 (MQ=255)

                                                                                        tGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACg                             <  1:701348/129‑1 (MQ=255)

                                                                                           cGCAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGttttt                >  2:770357/1‑139 (MQ=255)

                                                                                             cAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgttgt                    <  1:726226/133‑1 (MQ=255)

                                                                                             cAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtgttgt                    >  2:726226/1‑133 (MQ=255)

                                                                                             cAACACCATTAATCATCACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTAt              >  2:474193/1‑139 (MQ=255)

                                                                                                        atcatcACTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCa   <  1:474193/139‑1 (MQ=255)

                                                                                                         tcatcaCTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCaa  <  1:1294243/139‑1 (MQ=255)

                                                                                                               cTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc                                                                            <  2:52770/59‑1 (MQ=39)

                                                                                                               cTGAACCGCCGCCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc                                                                            >  1:52770/1‑59 (MQ=39)

CACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAA  >  NC_000913/2187916‑2188158

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