BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,477,418	T→C	10.5%	I24I (ATT→ATC)	PP_1291 →	PhoH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,477,418	0	T	C	10.5%	99.1 / 2.7	38	I24I (ATT→ATC)	PP_1291	PhoH family protein
Reads supporting (aligned to +/- strand): ref base T (13/21); new base C (2/2); total (15/23)
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

TACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAATGCCGCCAGGCCATCCGCCTGATAGAC  >  NC_002947/1477275‑1477556

tACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGa                                                                                                                                       <  3:32436/149‑1 (MQ=255)

         aGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAACCAATCCTGTATGTGCTCGATACCAACGTCCTGATCCACGACCCTAAcgc                                                                                                                              >  2:31145/1‑149 (MQ=255)

         aGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAACCAATCCTGTATGTGCTCGATACCAACGTCCTGATCCACGACCCTAAcgc                                                                                                                              >  2:36771/1‑149 (MQ=255)

            ggCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGtt                                                                                                                           <  7:131725/149‑1 (MQ=255)

             gCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTa                                                                                                                          <  5:265299/149‑1 (MQ=255)

             gCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAACAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGtt                                                                                                                           >  6:336073/1‑148 (MQ=255)

                     tctGCGGCGTAACAGACCACGCCGTGTGGCTTCGGCATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACtt                                                                                                                  >  1:67827/1‑149 (MQ=255)

                         ctgcgTAACAGACCACTCCGTGTGGCTTCGACATAAGGTGACCGAGTATGTATGACCAACGACGCAACCCTTCCTCCAGCAAGCCATTCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCga                                                                                                               <  1:243606/146‑1 (MQ=255)

                         cggcgTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTc                                                                                                                 <  3:261196/146‑1 (MQ=255)

                         cggcgTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCgag                                                                                                              <  3:49627/149‑1 (MQ=255)

                           gcgTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAAc                                                                                                                    >  8:85334/1‑141 (MQ=255)

                           gcgTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAAc                                                                                                                    <  7:85334/141‑1 (MQ=255)

                              tAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCa                                                                                                         <  2:34714/149‑1 (MQ=255)

                                aCAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCAcc                                                                                                       >  4:39348/1‑149 (MQ=255)

                                                  cTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGa                                                                                     >  8:157577/1‑149 (MQ=255)

                                                   ttCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACAAACGTCCTGATTCACGACCCTAACTCGTTACTGAACTTCGAGGAGAACAATGTCGCAATCCCGATGAc                                                                                    >  6:116998/1‑149 (MQ=255)

                                                      gACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAACAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACgg                                                                                  <  5:336073/148‑1 (MQ=255)

                                                            aGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAACCAATCCTGTATGTGCTCGATACCAACGTCCTGATCCACGACCCTAACGCGTTACTGAACTTCGAGGAGCNCCATGTCGCCATCCCGATGACGGTGCTgga                                                                           <  1:31145/149‑1 (MQ=255)

                                                            aGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAACCAATCCTGTATGTGCTCGATACCAACGTCCTGATCCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTgga                                                                           <  1:36771/149‑1 (MQ=255)

                                                                   cGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTc                                                                    <  6:322356/149‑1 (MQ=255)

                                                                         tGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACaa                                                               <  3:39348/148‑1 (MQ=255)

                                                                             tGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCa                                                          <  3:178559/149‑1 (MQ=255)

                                                                                  aaGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGACCTCGACAAGCTCAAGAcc                                                     >  8:245886/1‑149 (MQ=255)

                                                                                        ggaaCCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTg                                                                                >  3:65255/3‑116 (MQ=255)

                                                                                        gCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTg                                                                                <  4:65255/116‑1 (MQ=255)

                                                                                            cccTTCCTCCAGCAATCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAAgaga                                           <  5:116998/149‑1 (MQ=255)

                                                                                               ttCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTggag                                                                          <  7:169811/115‑1 (MQ=255)

                                                                                               ttCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTggag                                                                          >  8:169811/1‑115 (MQ=255)

                                                                                                   tccAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGCCCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATcgc                                    >  2:365334/1‑149 (MQ=255)

                                                                                                               aTCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACaga                                           >  5:367777/1‑130 (MQ=255)

                                                                                                               aTCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACaga                                           <  6:367777/130‑1 (MQ=255)

                                                                                                                   tGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAAt                            >  8:349836/1‑141 (MQ=255)

                                                                                                                   tGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAAt                            <  7:349836/141‑1 (MQ=255)

                                                                                                                                  ccAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAATGCCGCCAGGCCATCCGCCTGATa     <  7:245886/149‑1 (MQ=255)

                                                                                                                                     aCGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAATGCCGCCAGGCCATCCGCCTGATagac  <  4:353102/149‑1 (MQ=255)

                                                                                                                                       gTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAATGCCGCCAGGCCATCCGCCTGATagac  <  1:145796/147‑1 (MQ=255)

                                                                                                                                          cTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACaga                                           >  1:367658/1‑103 (MQ=255)

                                                                                                                                          cTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACaga                                           <  2:367658/103‑1 (MQ=255)

TACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGACCAAGGACGCAACCCTTCCTCCAGCAAGCCAATCCTGTATGTGCTCGATACCAACGTCCTGATTCACGACCCTAACGCGTTACTGAACTTCGAGGAGCACCATGTCGCCATCCCGATGACGGTGCTGGAGGAACTCGACAAGCTCAAGACCGGCAAACAGAGCATCGCCGCCGAATGCCGCCAGGCCATCCGCCTGATAGAC  >  NC_002947/1477275‑1477556

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