BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,849,930	T→C	V365A (GTG→GCG)	sppA →	protease IV, a signal peptide peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,849,930	0	T	C	100.0%	110.1 / NA	32	V365A (GTG→GCG)	sppA	protease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (20/12); total (20/12)

ACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGG  >  NC_000913/1849794‑1850065

aCGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCAGCGACGCTCGCCTTGACCCGAAAGCGaa                                                                                                                                      <  1:481476/140‑1 (MQ=255)

           tACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTcc                                                                                                                           <  2:583098/140‑1 (MQ=255)

                gTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCg                                                                                                                       <  1:122197/139‑1 (MQ=255)

                gTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCg                                                                                                                       >  2:122197/1‑139 (MQ=255)

                          cGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTaa                                                                                                                 >  2:762103/1‑135 (MQ=255)

                          cGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTaa                                                                                                                 <  1:762103/135‑1 (MQ=255)

                                 gtcTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAggcgg                                                                                                     <  2:543062/140‑1 (MQ=255)

                                          aaTGGCGCAATTATGGATGGCGTGGAAACTCGGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTAc                                                                                            <  2:289256/140‑1 (MQ=255)

                                          aaTGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTAc                                                                                            <  2:86819/140‑1 (MQ=255)

                                           aTGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAACGCGAAAGCGATTGTACGCCGTGTTAGTAGCAAAGGCGGCCGCGTTAcc                                                                                           >  2:421464/1‑140 (MQ=255)

                                           aTGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTAcc                                                                                           >  2:72363/1‑140 (MQ=255)

                                                       tggatggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTg                                                                               >  2:27687/1‑140 (MQ=255)

                                                            ggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAATCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcg                                                                          >  1:145154/1‑140 (MQ=255)

                                                            ggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcg                                                                          >  1:662488/1‑140 (MQ=255)

                                                              cGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGGGTGTTAATAGCCAAGGTGGCAGCGTAACCGGGTCTGAAGTGATTcgcg                                                                        >  2:275960/1‑140 (MQ=255)

                                                              cGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcgcg                                                                        >  2:14306/1‑140 (MQ=255)

                                                                aGGAAACTTAGGGGAATGTTGGCGGTGATACTACTTCGGCACAAATCCGCGATTCTCGCCTTGACCCGAAAACGACAGCGAATGTTCTGCGTGTTAATAGCCCACGCGCCAGCGTTAACGCGTCTGAAGTGATTCgcgcc                                                                      >  2:326622/1‑139 (MQ=255)

                                                                aGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGGACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGCCCTGCGTGTTAATAGTCCAAGCGCCAGCGTTAACGCGTCTGAAGTGATTCGCGCt                                                                      >  2:318122/1‑140 (MQ=255)

                                                                aGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCt                                                                      >  1:865997/1‑140 (MQ=255)

                                                                         aGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgca                                                             >  1:537395/1‑140 (MQ=255)

                                                                                tGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGccc                                                      >  1:573165/1‑140 (MQ=255)

                                                                                tGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACAGCGTCTGAAGTAATTCGCGCTAAACTGGCAGCAGccc                                                      >  1:287447/1‑140 (MQ=255)

                                                                                   tGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCggg                                                   <  2:865997/140‑1 (MQ=255)

                                                                                           aTACCACTGCGGCACAACTCCGCGACGCTCGCCTTGACCAGAAAGCGACAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGt                                           <  2:287447/140‑1 (MQ=255)

                                                                                                  tGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                    >  2:669759/1‑140 (MQ=255)

                                                                                                   gCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCACTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgt                                   >  2:440188/1‑140 (MQ=255)

                                                                                                    cGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgtg                                  >  1:742908/1‑140 (MQ=255)

                                                                                                     ggCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTgg                                 >  1:566323/1‑140 (MQ=255)

                                                                                                     ggCACAAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTgg                                 >  1:499900/1‑140 (MQ=255)

                                                                                                       cacaAATCCGCGACGCTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGCGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGtt                               >  1:448745/1‑140 (MQ=255)

                                                                                                                      cTCGCCTTGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggc                >  1:795013/1‑140 (MQ=255)

                                                                                                                            ttGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcg          <  2:742908/140‑1 (MQ=255)

                                                                                                                            ttGACCCGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcg          <  2:795013/140‑1 (MQ=255)

                                                                                                                                cccGAAAGCGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcgggatggcggcat      <  2:448745/140‑1 (MQ=255)

                                                                                                                                    aaagcgaaagcgATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtgg  >  1:498744/1‑140 (MQ=255)

ACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGG  >  NC_000913/1849794‑1850065

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