BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,688,071	C→G	12.2%	R349R (CGC→CGG)	PP_4148 →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,688,071	0	C	G	12.2%	77.4 / 4.3	33	R349R (CGC→CGG)	PP_4148	ABC transporter permease
Reads supporting (aligned to +/- strand): ref base C (15/14); new base G (2/2); total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

GCCTGATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAACTG  >  NC_002947/4687941‑4688216

gcctgATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCGTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCgg                                                                                                                                 <  1:104284/149‑1 (MQ=255)

gcctgATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCgg                                                                                                                                 <  5:22341/149‑1 (MQ=255)

    gatgaGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAg                                                                                                                              >  4:178278/1‑148 (MQ=255)

     atgaGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCa                                                                                                                            <  3:232654/149‑1 (MQ=255)

             cGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACAt                                                                                                                     <  3:178278/148‑1 (MQ=255)

                    gCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTg                                                                                                              >  8:141590/1‑148 (MQ=255)

                         cAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTccccc                                                                                                        >  6:322495/1‑149 (MQ=255)

                            ccGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTc                                                                                                     <  7:141590/149‑1 (MQ=255)

                                            gtCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTc                                                                                      <  5:322495/148‑1 (MQ=255)

                                                   gCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGACAGCCGGGAGCACTGACATGGCCTTGTCCACCCTCAACCGCCGGCGCTTCGAGCGtt                                                                               >  3:92126/1‑148 (MQ=255)

                                                   gCACCCGGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGttt                                                                              <  7:245329/149‑1 (MQ=255)

                                                                cTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCACCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAacgc                                                                 >  4:225139/1‑149 (MQ=255)

                                                                                 tggtggTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCa                                                                            <  1:161083/121‑1 (MQ=255)

                                                                                 tggtggTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCa                                                                            >  2:161083/1‑121 (MQ=255)

                                                                                         aaaCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtgg                                                        >  8:205065/1‑133 (MQ=255)

                                                                                         aaaCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtgg                                                        <  7:205065/133‑1 (MQ=255)

                                                                                           aCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTc                                      >  6:177396/1‑149 (MQ=255)

                                                                                                 cGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGAt                                 >  4:36226/1‑148 (MQ=255)

                                                                                                   gCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGAAATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCt                              >  6:162362/1‑149 (MQ=255)

                                                                                                     gacctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGt                            >  1:11040/1‑149 (MQ=255)

                                                                                                       cctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCAGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTACGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTc                          >  7:70863/1‑149 (MQ=255)

                                                                                                        ctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTTGTCGCTGTGGCTGTTCCTGATTCTGTTCg                         <  5:162362/149‑1 (MQ=255)

                                                                                                                   cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCtgt                                              >  4:77565/1‑117 (MQ=255)

                                                                                                                   cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCtgt                                              <  3:77565/117‑1 (MQ=255)

                                                                                                                   cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTg              >  4:141349/1‑149 (MQ=255)

                                                                                                                   cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTg              >  4:136357/1‑149 (MQ=255)

                                                                                                                     aagTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAc                                                                    >  7:28700/3‑93 (MQ=255)

                                                                                                                     aagTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAc                                                                    <  8:28700/91‑1 (MQ=255)

                                                                                                                      agTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtggtg                                                      <  2:247476/105‑1 (MQ=255)

                                                                                                                      agTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtggtg                                                      >  1:247476/2‑106 (MQ=255)

                                                                                                                            tCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGaa     <  3:141349/149‑1 (MQ=255)

                                                                                                                            tCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGaa     <  3:136357/149‑1 (MQ=255)

                                                                                                                              cgcgCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAAc    <  2:11040/148‑1 (MQ=255)

                                                                                                                               gcgcATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAACTg  <  3:36226/149‑1 (MQ=255)

GCCTGATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAACTG  >  NC_002947/4687941‑4688216

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