BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,332,792	G→A	11.7%	Q105* (CAA→TAA)	PP_2930 ←	L‑serine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,332,792	0	G	A	11.7%	88.0 / 2.8	35	Q105* (CAA→TAA)	PP_2930	L‑serine dehydratase
Reads supporting (aligned to +/- strand): ref base G (14/16); new base A (2/2); total (16/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 = 6.05e-01

CCGGCTTCAGGCCGGCCTCGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCTGAC  >  NC_002947/3332645‑3332928

ccGTCTTCAGGCCGGCCTCGGCCACTTCTTGCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGc                                                                                                                                         <  5:78267/149‑1 (MQ=255)

         ggccggccTCGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCaccaccacc                                                                                                                                >  4:132715/1‑149 (MQ=255)

              gccTCGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGc                                                                                                                           <  4:171784/149‑1 (MQ=255)

                      cACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTc                                                                                                                     >  6:244082/1‑147 (MQ=255)

                        cTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGTGCGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGTACTGCCTTGCACCACCACCTTGGCGTCTTCCAg                                                                                                                  >  1:254397/1‑149 (MQ=255)

                        cTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGTACTGCCTTGCACCACCACCTTGGCGTCTTCCa                                                                                                                   <  2:254397/147‑1 (MQ=255)

                           cATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTg                                                                                                               >  2:238053/1‑148 (MQ=255)

                               caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTgcagca                                                                                                          <  3:132715/149‑1 (MQ=255)

                               caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTgcagc                                                                                                           >  3:105568/1‑148 (MQ=255)

                               caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGACTTGCACCACCAACTTGGAGTCTTCAAGGTgcagc                                                                                                           >  4:45690/1‑148 (MQ=255)

                                accaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCaa                                                                                                         <  3:17927/149‑1 (MQ=255)

                                 ccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGGCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAAc                                                                                                        <  3:45690/149‑1 (MQ=255)

                                         gCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTg                                                                                                  <  7:306173/147‑1 (MQ=255)

                                         gCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGCGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTg                                                                                                  <  1:296189/147‑1 (MQ=255)

                                             gCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGt                                                                                            >  6:182946/1‑149 (MQ=255)

                                              cccGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGtt                                                                                           <  3:157027/149‑1 (MQ=255)

                                                cGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCg                                                                                         >  1:237706/1‑149 (MQ=255)

                                                             tGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCTTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCgg                                                                            <  3:66022/149‑1 (MQ=255)

                                                             tGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCgg                                                                            <  1:238053/149‑1 (MQ=255)

                                                             tGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCgg                                                                            <  4:105568/149‑1 (MQ=255)

                                                             tGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCgg                                                                            <  5:244082/149‑1 (MQ=255)

                                                             tGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCgg                                                                            <  2:237706/149‑1 (MQ=255)

                                                                gtcgt‑cGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCACCTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTa                                                                         >  7:144450/1‑149 (MQ=255)

                                                                                    aaGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCAc                                                      <  2:89531/149‑1 (MQ=255)

                                                                                                                   ccAGGGCGTTGGACTCTGTCCAGGAACTGCCTTTCACCACCACCTTCTCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGc                        <  8:144450/148‑1 (MQ=255)

                                                                                                                              gCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCAtt            <  5:182946/149‑1 (MQ=255)

                                                                                                                              gCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGCTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCAtt            >  2:242098/1‑149 (MQ=255)

                                                                                                                                    tGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCggt                                                                                       <  7:17181/68‑1 (MQ=255)

                                                                                                                                    tGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCggt                                                                                       >  8:17181/1‑68 (MQ=255)

                                                                                                                                     gCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCt     >  3:161633/1‑149 (MQ=255)

                                                                                                                                      ccAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCt     >  5:3369/1‑148 (MQ=255)

                                                                                                                                      ccAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCt     >  6:99331/1‑148 (MQ=255)

                                                                                                                                       cAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCTGa   >  4:215570/1‑149 (MQ=255)

                                                                                                                                        aGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCAGAc  >  2:294555/1‑149 (MQ=255)

                                                                                                                                                ccTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTg                         <  6:313926/118‑1 (MQ=255)

CCGGCTTCAGGCCGGCCTCGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCGGTGGTGGTTTCGGGTATCACCACGGTCACCGGCACCCCCAGCTTGCGTCCGGCATAGGCGACTGCCAAACCGGCATTGCCACCTGAC  >  NC_002947/3332645‑3332928

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