BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,187,198	G→T	10.6%	L167L (CTG→CTT)	PP_3683 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,187,198	0	G	T	10.6%	101.4 / 2.7	38	L167L (CTG→CTT)	PP_3683	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (17/17); new base T (2/2); total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CTCCTTGAGTCAGCTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAGCTTCCTTAAA  >  NC_002947/4187049‑4187345

cTCCTTGAGTCAGCTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTg                                                                                                                                                      >  1:609173/1‑150 (MQ=255)

 tCCTTGGGTCAGCTTGGACGCAAGCCATCGAGTCGGATGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGGGCTTGGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTGTCGCCGAGAACATCTTGGATAGGCACGATCTgg                                                                                                                                                     <  2:473230/150‑1 (MQ=255)

    ttGAGTCAGCTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCt                                                                                                                                                  >  2:952560/1‑150 (MQ=255)

    ttGAGTCAGCTTGGACGCAAGCCATCCAGTAGGACGGCAATCATGCCTGGACGGTCAATTCTGCTAGTGTGGACGAGCTTCCCGCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCt                                                                                                                                                  >  2:357975/1‑150 (MQ=255)

    ttGAGTCAGCTTGGACGCAAGCCATACAGTCGGACAGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCACAGGTCGTGCAGAATTTGAGTTTTCGCCTCTCGCCGAGAACATGTTGTATAGGCACGTTCTGGGCt                                                                                                                                                  >  2:520036/1‑150 (MQ=255)

       ggTCGGCTTGGACGCAAGCCACCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTAt                                                                                                                                               <  1:158064/149‑1 (MQ=255)

           aGCTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGa                                                                                                                                           <  2:882706/150‑1 (MQ=255)

             cTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTCTGCTAGTGTGGACTTGCTTCGCTCTGGTCGTTCAGATTTTGAGTTTTCGACTTGCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGt                                                                                                                                         <  1:570810/150‑1 (MQ=255)

                             ccAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCGCGATCTGGGCTTATCCGAGTTCAAGGCCTTCCAGCg                                                                                                                         >  1:485578/1‑150 (MQ=255)

                                   ggACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCacac                                                                                                                   >  2:790890/1‑150 (MQ=255)

                                      cGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGtt                                                                                                                <  1:790890/150‑1 (MQ=255)

                                       ggCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGGTAGGAAAGCTTTGGCCTAGGTCTAGCTCAAGGACGTCCAGCGTCACACGTTg                                                                                                               >  2:398046/1‑150 (MQ=255)

                                         cAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGAt                                                                                                             >  1:547865/1‑150 (MQ=255)

                                         cAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGAt                                                                                                             >  1:208606/1‑150 (MQ=255)

                                              aTGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGt                                                                                                        >  1:1059552/1‑150 (MQ=255)

                                               tGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTa                                                                                                       <  2:940935/150‑1 (MQ=255)

                                                gCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTAc                                                                                                      <  2:430954/150‑1 (MQ=255)

                                                 ccTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTg                                                                                                                                                      >  2:265137/1‑101 (MQ=255)

                                                 ccTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTg                                                                                                                                                      <  1:265137/101‑1 (MQ=255)

                                                  cTGNACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTTGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGa                                                                                                    >  2:698889/1‑150 (MQ=255)

                                                  cTGNACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTTGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGa                                                                                                    >  2:695694/1‑150 (MQ=255)

                                                          tCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGAt                                                                                                             >  2:622072/1‑133 (MQ=255)

                                                          tCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGAt                                                                                                             <  1:622072/133‑1 (MQ=255)

                                                           tctttttGCTAGTGTGGACGAGCTTCGCTCAGTTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCNAGTTNAAGGACTTCCAGCGNCACACGTTGATCACGTACGAGGCTTTGga                                                                                           <  1:154566/147‑1 (MQ=255)

                                                               tttGCTAGTGTGGACGAGCTTCGCTCAGTTCGTTCAGATTTTGAGGTTTCGGCTATCGCCGAGAACATCTTGGATAGGCACGAGCTGGGCTTATCCGAGTTCAAAGACTTCCAGCGTCACAGGTTGCTTACGTCCGCGGTATTGgacaga                                                                                       >  2:154566/1‑147 (MQ=255)

                                                                           gACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATtgt                                                                           >  1:672986/1‑150 (MQ=255)

                                                                           gACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATtgt                                                                           <  2:1059552/150‑1 (MQ=255)

                                                                           gACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATtgt                                                                           >  1:211149/1‑150 (MQ=255)

                                                                             cGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTTGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGa                                                                         <  1:695694/150‑1 (MQ=255)

                                                                             cGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTTGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGa                                                                         <  1:698889/150‑1 (MQ=255)

                                                                                           gTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑Ct                                                          <  2:728663/150‑1 (MQ=255)

                                                                                             cGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑Cttt                                                        <  2:547865/150‑1 (MQ=255)

                                                                                                                 gACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGc                                    >  2:925031/1‑150 (MQ=255)

                                                                                                                 gACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGc                                    >  2:526929/1‑150 (MQ=255)

                                                                                                                        gCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGAAGACTCACGAAATAATGATAGTGC‑TCTTTACTTTGACGATCGCCTAGTGGAGGCTcaaaaa                             >  2:560906/1‑146 (MQ=255)

                                                                                                                           gagaACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAgg                          <  1:925031/150‑1 (MQ=255)

                                                                                                                                  tCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGt                      <  1:934396/147‑1 (MQ=255)

                                                                                                                                  tCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGt                      >  2:934396/1‑147 (MQ=255)

                                                                                                                                         tAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGCCTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAg            <  2:872024/150‑1 (MQ=255)

                                                                                                                                           ggCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCCCACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAGCt          <  1:398046/150‑1 (MQ=255)

                                                                                                                                              aCGATCTGCGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAGCTTcc       <  2:609173/150‑1 (MQ=255)

                                                                                                                                                   cTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAGCTTCCTTaaa  <  2:208606/150‑1 (MQ=255)

CTCCTTGAGTCAGCTTGGACGCAAGCCATCCAGTCGGACGGCAATCATGCCTGGACGGTCAATTTTGCTAGTGTGGACGAGCTTCGCTCAGGTCGTTCAGATTTTGAGTTTTCGACTATCGCCGAGAACATCTTGGATAGGCACGATCTGGGCTTATCCGAGTTCAAGGACTTCCAGCGTCACACGTTGATCACGTACGAGGCTTTGGACGACTCACGAAATTGTGACGCAGCGTCGGT‑CTTTCACGATCGCCTAGTGGACGCTCACACGAGGCGGTCTGATAAAAGCTTCCTTAAA  >  NC_002947/4187049‑4187345

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