BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,055,557	C→T	10.5%	R23H (CGT→CAT)	PP_2667 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,055,557	0	C	T	10.5%	99.3 / 2.7	38	R23H (CGT→CAT)	PP_2667	ABC transporter permease
Reads supporting (aligned to +/- strand): ref base C (16/18); new base T (2/2); total (18/20)
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.39e-01

CAGGCCGGGGATGATGTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCGAGGTCACCGCCGTGCTCAAGGC  >  NC_002947/3055415‑3055697

caggccGGGGATGATGTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTg                                                                                                                                           >  2:250673/1‑147 (MQ=255)

 aggccGGGGATGATGTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTgcag                                                                                                                                        >  6:243137/1‑149 (MQ=255)

     cGGGGATGATGTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCAc                                                                                                                                     >  3:104728/1‑148 (MQ=255)

            gatgTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAgcg                                                                                                                               >  3:251510/1‑147 (MQ=255)

             atgTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAg                                                                                                                            <  5:243137/149‑1 (MQ=255)

                    ccTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTcg                                                                                                                     >  8:1746/1‑149 (MQ=255)

                                    ggTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGt                                                                                                     <  4:251510/149‑1 (MQ=255)

                                          gTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGaaa                                                                                                <  4:386292/148‑1 (MQ=255)

                                           tAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACa                                                                                              >  6:319760/1‑149 (MQ=255)

                                             gggCTCGATGATGGCAATGCACAGCGATGCGCGAAAC‑CTTGCAGAGAACACCAGCAACCACAGCAGCGGTCGCACCATCGCGCTGAGCAGGAGTGTACGCTGCAGCACTAAGCGCAGCCATTCGCGGATCATGATGCCGTTGAAACAt                                                                                             <  6:327456/148‑1 (MQ=255)

                                                   gatgatGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCa                                                                                        <  5:241800/147‑1 (MQ=255)

                                                      gatgGCAATGCCCAGCGCTGCGCGAAACACCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTATGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGt                                                                                      >  1:24894/1‑147 (MQ=255)

                                                           cAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGtt                                                                               >  7:182123/1‑148 (MQ=255)

                                                              tGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATc                                                                          <  5:319760/149‑1 (MQ=255)

                                                                          cgcgAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGc                                                                <  7:1746/147‑1 (MQ=255)

                                                                               aaC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAg                                                         <  8:72918/149‑1 (MQ=255)

                                                                               aaC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAg                                                         >  2:72375/1‑149 (MQ=255)

                                                                                       cAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGt                                                  <  4:104728/149‑1 (MQ=255)

                                                                                         tcGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTg                                                <  3:230565/148‑1 (MQ=255)

                                                                                          cGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTATGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ACCGGGCGTTGGCTCCTGAGGGGGCGTTgg                                               <  2:24894/149‑1 (MQ=255)

                                                                                          cGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTgg                                               >  1:248934/1‑149 (MQ=255)

                                                                                           gAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTg                                                >  4:230565/1‑147 (MQ=255)

                                                                                            aaCACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTc                                             <  7:99034/149‑1 (MQ=255)

                                                                                                 cAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTATGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCg                                        >  8:220198/1‑149 (MQ=255)

                                                                                                 cAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGt                                              >  2:349128/1‑143 (MQ=255)

                                                                                                 cAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGt                                              <  1:349128/143‑1 (MQ=255)

                                                                                                 cAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCg                                        >  4:211726/1‑149 (MQ=255)

                                                                                                  aGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCgg                                       >  2:185786/1‑149 (MQ=255)

                                                                                                  aGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCAATCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCgg                                       <  1:250673/149‑1 (MQ=255)

                                                                                                   gCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTCAATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCgg                                       <  1:72375/149‑1 (MQ=255)

                                                                                                    cAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAGCATTGCCAGGAAGCGTTC‑ATCGGGCGTTGGCTGGTGAGGGGGCGTTGGTCAGGCGGGc                                     >  4:117483/1‑149 (MQ=255)

                                                                                                        cacaGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCgg                                                                        <  8:129465/110‑1 (MQ=255)

                                                                                                          cagcagCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGc                               <  1:185786/149‑1 (MQ=255)

                                                                                                           agcagcGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCg                              >  6:275451/1‑149 (MQ=255)

                                                                                                                 gggCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCg                        <  3:211726/149‑1 (MQ=255)

                                                                                                                  ggCGCACCAGCGCGCTGAGCAGGCGTGTATGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCGa                       <  7:220198/149‑1 (MQ=255)

                                                                                                                                              aCGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCGAGGTCACCGCCGTGCTCAAGGc  <  7:229630/142‑1 (MQ=255)

                                                                                                                                              aCGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCGAGGTCACCGCCGTGCTCAAGGc  >  8:229630/1‑142 (MQ=255)

CAGGCCGGGGATGATGTACACCTCGTAGGGGATGTAGGTGTCGTAGGGCTCGATGATGGCAATGCCCAGCGCTGCGCGAAAC‑CCGGCAGCGAACACCAGCAACCACAGCAGCGGGCGCACCAGCGCGCTGAGCAGGCGTGTACGCTGCAGCACGAAGCGCAGCCATTCGCGCAGCATGATGCCGTTGAAACATTGCCAGTAAGCGTTC‑ATCGGGCGTTGGCTCCTGAGGGGGCGTTGGTCAGGCGGGCAAAGGCGTGATCGAGGTCACCGCCGTGCTCAAGGC  >  NC_002947/3055415‑3055697

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