BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	481,213	A→G	11.7%	L383P (CTG→CCG)	PP_0396 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	481,213	0	A	G	11.7%	85.0 / 3.1	34	L383P (CTG→CCG)	PP_0396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (14/16); new base G (2/2); total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGA  >  NC_002947/481075‑481351

ctGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGgtg                                                                                                                                  >  6:31260/1‑149 (MQ=255)

             gTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTg                                                                                                                     <  8:272303/149‑1 (MQ=255)

                      cTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGACGATTCCAACTGCTGCTGGGAGGAAGTGTCGGGGAAAGCATCTCCGATACGCACGTATTAATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACg                                                                                                            >  2:29998/1‑149 (MQ=255)

                              cgcTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTgta                                                                                                    <  4:288049/149‑1 (MQ=255)

                               gcTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTgtag                                                                                                   <  4:152216/149‑1 (MQ=255)

                                   gAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACt                                                                                               >  3:20157/1‑149 (MQ=255)

                                   gAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACt                                                                                               >  3:213804/1‑149 (MQ=255)

                                     agagTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTg                                                                                              >  3:138365/1‑148 (MQ=255)

                                     agacgtCACCGAAGACCGCGTAGATACCGCCGGCCGATACCAACTGCGGCTGGGCGAACGTGGCGCGGAAAGCAACTCCGAAACGCTCGTATTCATACCACACCGCCGGGTGTTCACGGGGGCTGATCTCAACGTAAGTGTAGTACTGc                                                                                             <  4:44879/144‑1 (MQ=255)

                                      gagTTCACGGAAGACCGTGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGca                                                                                            <  7:239360/149‑1 (MQ=255)

                                        gTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGcaca                                                                                          <  5:31260/149‑1 (MQ=255)

                                                           gATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTgg                                                                       <  5:163067/149‑1 (MQ=255)

                                                           gATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTg                                                                        >  6:163067/1‑148 (MQ=255)

                                                              atCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCgg                                                                                     <  2:112015/132‑1 (MQ=255)

                                                              atCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCgg                                                                                     >  1:112015/1‑132 (MQ=255)

                                                               tCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGgaga                                                                    <  4:138365/148‑1 (MQ=255)

                                                                 gccggccgATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGa                                                                 <  2:220067/149‑1 (MQ=255)

                                                                             cAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCTTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGAGGCATGATCTGTTTGGAGAGGATGTCGCgg                                                         >  7:234258/1‑145 (MQ=255)

                                                                                ctgctgctgTGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATgg                                                  <  4:20157/149‑1 (MQ=255)

                                                                                                   tCGGTGAAAGCATCTGCGATACGCTCGTATTACTACCACAGGGCCTGGTGTTCAAGGGGGGTGATCTCAAGGTAAGTGTAGTACTTCACTTGCGTCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGtt                               <  7:279776/149‑1 (MQ=255)

                                                                                                              aTCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCAt                       >  5:213966/1‑146 (MQ=255)

                                                                                                              aTCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCa                     <  6:213966/148‑1 (MQ=255)

                                                                                                                           tCGTATTCATACCACGGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCt              <  4:110938/142‑1 (MQ=255)

                                                                                                                           tCGTATTCATACCACGGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCt              >  3:110938/1‑142 (MQ=255)

                                                                                                                           tCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTcg                                                            <  7:68927/96‑1 (MQ=255)

                                                                                                                           tCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTcg                                                            >  8:68927/1‑96 (MQ=255)

                                                                                                                             gTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCTGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCa      <  4:213804/148‑1 (MQ=255)

                                                                                                                              tATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATa    >  3:130168/1‑149 (MQ=255)

                                                                                                                               aTTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATATGTTTGGAGAGTATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATAAGAAGCCTGCGCGGCa      >  5:76512/1‑146 (MQ=255)

                                                                                                                               aTTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGCGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCa      >  5:237230/1‑146 (MQ=255)

                                                                                                                                 tCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGa  >  4:160520/1‑148 (MQ=255)

                                                                                                                                 tCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGa  <  3:160520/148‑1 (MQ=255)

                                                                                                                                   aTACCACGGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGa  >  4:143940/1‑146 (MQ=255)

                                                                                                                                   aTACCACGGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGa  <  3:143940/146‑1 (MQ=255)

CTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTAGATATCGCCGGCCGATACCAACTGCTGCTGGGCGAACGTGTCGGGGAAAGCATCTCCGATACGCTCGTATTCATACCACAGCGCCTGGTGTTCACGGGGGGTGATCTCAACGTAAGTGTAGTACTGCACATGCGGCATGATCTGTTTGGAGAGGATGTCGCGGCAGATGGGTGAGTCGTCGTTCCAGTTGTCGCCATCAGAAGCCTGCGCGGCATAGA  >  NC_002947/481075‑481351

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