BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,297,843	G→A	11.7%	A760V (GCC→GTC)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,297,843	0	G	A	11.7%	85.0 / 3.3	34	A760V (GCC→GTC)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base G (15/15); new base A (2/2); total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGAT  >  NC_002947/2297706‑2297982

tctgcAGGCGCTGGGAAATCTGCTTGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCATTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTg                                                                                                                                   <  2:218821/147‑1 (MQ=255)

gctgcAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTg                                                                                                                                   <  2:222930/148‑1 (MQ=255)

   gcAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGc                                                                                                                               >  3:189629/1‑149 (MQ=255)

               cAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAgg                                                                                                                   >  8:275841/1‑149 (MQ=255)

                         gTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTggg                                                                                                         >  5:127765/1‑149 (MQ=255)

                            gacgccGAGGCAGGCAGTGGCCGGGCCGGCGTTGCGGGCCTTGAGCGCTTCCTCGGGTAGTACGCCGGCCAGGGCGTTCAGGCCTTGCGGCAACAGCTGCTCGTCGGTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCt                                                                                                      <  2:20445/144‑1 (MQ=255)

                             aGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCtg                                                                                                     <  6:127765/149‑1 (MQ=255)

                               ggCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCtgtt                                                                                                   >  7:51940/1‑149 (MQ=255)

                                gCGAGGAATGCATTGGGCGGGTCGTGGTTGAGGGCCTTGAGCGCTTCATCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGGTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCtgttg                                                                                                  <  4:229317/149‑1 (MQ=255)

                                  gAGGAAGGCAGTGGCCGGGTCGGCGTTGACGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTgttgtt                                                                                                <  8:51940/149‑1 (MQ=255)

                                    ggAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGttgttg                                                                                               <  4:189629/148‑1 (MQ=255)

                                    ggAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTgg                                                                                              <  7:275841/149‑1 (MQ=255)

                                       ctgctttGGCGGGGTCGTCTTTGAGGGGCTTGAGGGCTTCGTCTGGTTGTACGCGGGCGAGGTCGTTCAGGCCTTGCTGCAACAGGTGCTCGTCGTTGGGCAGGGCTTGGTGCTGTTGCTCGAGGTGGCGCTGGGCCTGTTGTTGGGc                                                                                            <  4:130130/143‑1 (MQ=255)

                                        ggCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGcgt                                                                                          >  6:195791/1‑149 (MQ=255)

                                          cATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGc                                                                                            >  3:130130/1‑145 (MQ=255)

                                                ccGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTggg                                                                                    >  6:11227/1‑147 (MQ=255)

                                                                ggCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGt                                                                  >  4:147762/1‑149 (MQ=255)

                                                                           cTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGcc                                                       >  5:196960/1‑149 (MQ=255)

                                                                            ttCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                             <  4:129894/142‑1 (MQ=255)

                                                                            ttCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAgg                                                             >  3:129894/1‑142 (MQ=255)

                                                                                tctGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTgg                                                                   <  4:212684/132‑1 (MQ=255)

                                                                                tctGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTgg                                                                   >  3:212684/1‑132 (MQ=255)

                                                                                                   ggTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAgg                               <  5:195791/149‑1 (MQ=255)

                                                                                                     tCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTCCTCCAGGTGGCGCTGGGCCTGTTGTTGCGCGTCGTGGGCAGCCTGCAATTGCTgg                                                                   >  2:286305/1‑111 (MQ=255)

                                                                                                             gCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGc                     <  2:133084/149‑1 (MQ=255)

                                                                                                               gatcttgcATCAGCTGGTCGTCGTTGACCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGt                   <  4:30023/144‑1 (MQ=255)

                                                                                                               cTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGt                   >  3:170590/1‑149 (MQ=255)

                                                                                                                 tgctgcAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCg                                               <  5:209176/119‑1 (MQ=255)

                                                                                                                 tgctgcAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCg                                               >  6:209176/1‑119 (MQ=255)

                                                                                                                   ctgcAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGTCGCTCTGGCGTTtc                <  5:11227/148‑1 (MQ=255)

                                                                                                                     gcAACAGCTGCTCGTCGTTGACCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGt                   <  4:33919/143‑1 (MQ=255)

                                                                                                                     gcAACAGCTGCTCGTCGTTGACCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGt                   >  3:33919/1‑143 (MQ=255)

                                                                                                                     gcAACAGATGCTCGTCGTTGACCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCAGTTGTTGGGCGTCGTGGCCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGAGt                   >  3:30023/1‑143 (MQ=255)

                                                                                                                                tcgtcgTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGAt  >  2:107050/1‑149 (MQ=255)

GCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGGTAGTACGCCGGCCAGGTCGTTCAGGCCTTGCTGCAACAGCTGCTCGTCGTTGGCCAGGGCTTGGTGCTGTTGCTCCAGGTGGCGCTGGGCCTGTTGTTGGGCGTCGTGGGCAGCCTGCAATTGCTGGTTGAGGCGGGCCGCGTCCCGTTGCAGGGCGAGCAGGGCGCTCTGGCGTTTCTCGTCCTGGCTGAT  >  NC_002947/2297706‑2297982

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: `‑`