BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,224,350	G→A	11.1%	R41C (CGC→TGC)	PP_2825 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,224,350	0	G	A	11.1%	92.2 / 2.6	36	R41C (CGC→TGC)	PP_2825	integrase
Reads supporting (aligned to +/- strand): ref base G (13/19); new base A (2/2); total (15/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTACCACGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGCGATTTGATCTCAGT  >  NC_002947/3224208‑3224494

gTTACCACGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAgg                                                                                                                                             <  5:81344/148‑1 (MQ=255)

gTTACCACGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAg                                                                                                                                              <  1:108620/147‑1 (MQ=255)

gTTACCACGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGt                                                                                                                                            >  6:266252/1‑149 (MQ=255)

       cGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGGGaag                                                                                                                                     >  3:94724/1‑147 (MQ=255)

       cGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCg                                                                                                                                        <  4:94724/146‑1 (MQ=255)

           gCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTAtct                                                                                                                                 <  1:134673/149‑1 (MQ=255)

             gAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTAtctc                                                                                                                                >  5:85070/1‑148 (MQ=255)

               aCCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCgat                                                                                                                             <  7:177269/149‑1 (MQ=255)

                                 ttGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGa                                                                                                            <  6:85070/148‑1 (MQ=255)

                                    cTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGAtctc                                                                                                        >  8:20383/1‑149 (MQ=255)

                                          gCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGAc                                                                                                  <  3:287439/149‑1 (MQ=255)

                                              cttagcgcTCCACCCAGCCGCTCCCATCGGCGAAACCTATCACGGCTCGAATGATGGTGCGTACGCTTTTGGACATGGCCAGTGTGAGTGAGAGTCGCAAGGTGGCGATTATCTCGATGATGTTATCCAGCTTTATCTCGCTGACgggg                                                                                              <  2:302808/145‑1 (MQ=255)

                                              cGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACggg                                                                                               <  5:266252/148‑1 (MQ=255)

                                               gTTGCGCTCCACCCAGCCCCGCCCCACGGCGAAACCGAGCACGGCACGAATGATGGTGCGTAAGCCTTTGGCCAAGGCCAGTGTGCCAGCGCGGAGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGAAAGAGCTGACGGGGa                                                                                             >  8:30453/1‑149 (MQ=255)

                                                     cTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTc                                                                                       <  5:61214/149‑1 (MQ=255)

                                                      tCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTcc                                                                                      >  8:6992/1‑149 (MQ=255)

                                                        cACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGa                                                                                    >  8:122997/1‑149 (MQ=255)

                                                                cccgccccTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTGGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACGTCCGAGCGTGGGc                                                                            >  5:223261/1‑149 (MQ=255)

                                                                 ccgccccTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCACAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTggg                                                                             >  3:253412/1‑147 (MQ=255)

                                                                     gccTCGGCGATACCGAGCACGGCACGATTGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGAt                                                                       <  6:35525/148‑1 (MQ=255)

                                                                      ccTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCATGGGCAGGAtg                                                                      >  3:245674/1‑149 (MQ=255)

                                                                            cGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGt                                                                                >  8:242015/1‑133 (MQ=255)

                                                                            cGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGt                                                                                <  7:242015/133‑1 (MQ=255)

                                                                                           cacgAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACgg                                                                                                <  8:189093/102‑1 (MQ=255)

                                                                                           cacgAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACgg                                                                                                >  7:189093/1‑102 (MQ=255)

                                                                                             cgAATGATGGGGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGaa                                               <  8:15882/149‑1 (MQ=255)

                                                                                               aaTGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGc                                             <  7:122997/149‑1 (MQ=255)

                                                                                                 tgatgGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGccc                                           >  2:114404/1‑149 (MQ=255)

                                                                                                             cgcgTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGAttt                                <  7:127504/148‑1 (MQ=255)

                                                                                                              gcgTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGc                              <  7:6992/149‑1 (MQ=255)

                                                                                                                     gccatggccaGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAg                       <  1:114404/149‑1 (MQ=255)

                                                                                                                            ccaGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGc                >  4:140852/1‑149 (MQ=255)

                                                                                                                             caGTGTGCGTGCGCGGCACAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAgg                                     <  6:85145/127‑1 (MQ=255)

                                                                                                                             caGTGTGCGTGCGCGGCACAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAgg                                     >  5:85145/1‑127 (MQ=255)

                                                                                                                                tgtgCGTGCGCGGCGCAAGGGGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGCGAtt            <  7:30453/149‑1 (MQ=255)

                                                                                                                                    cgtgcgCGGCACAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGCGATTTGa         <  4:253412/148‑1 (MQ=255)

                                                                                                                                         gcggcgCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGt                  <  3:62549/134‑1 (MQ=255)

                                                                                                                                         gcggcgCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGt                  >  4:62549/1‑134 (MQ=255)

                                                                                                                                          cggcgCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCATGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGCGATTTGATCTCAGt  <  4:245674/149‑1 (MQ=255)

GTTACCACGTGGCGAACCTTGAGTTCCTCAGTGTTGCTCAGGGCTACGTTGCGCTCCACCCAGCCCCGCCCCTCGGCGAAACCGAGCACGGCACGAATGATGGTGCGTACGCGTTTGGCCATGGCCAGTGTGCGTGCGCGGCGCAAGGTGGCGATTATCTCGATGATGTGATCCAGCTTGATCTCGCTGACGGGGATACTTCCGATCGTGGGCAGGATGTGATTGTTCAACGCGCCGCGGAAGCCCGCGAGGGATTTGCTCACGAGGTCGTGCGATTTGATCTCAGT  >  NC_002947/3224208‑3224494

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