BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,141,130	C→A	12.5%	A16E (GCG→GAG)	PP_2757 →	sugar‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,141,130	0	C	A	12.5%	87.4 / 3.3	32	A16E (GCG→GAG)	PP_2757	sugar‑binding protein
Reads supporting (aligned to +/- strand): ref base C (17/11); new base A (2/2); total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

AGCGCCGAACCGGCACTGCCCACTCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGCCTTGGC  >  NC_002947/3140992‑3141275

aGCGCCGAACCGGCACTGCCCACTCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTg                                                                                                                                         <  6:85707/149‑1 (MQ=255)

                    cACTCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCgagg                                                                                                                            >  7:190697/1‑142 (MQ=255)

                    cACTCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCgagg                                                                                                                            <  8:190697/142‑1 (MQ=255)

                       tCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGt                                                                                                                    >  1:115340/1‑147 (MQ=255)

                              cGACTGCGCGGCGACCTCTACCTCTACGGGGCCGCGTTGGCAGATGAGCTAAGGGCGCAGGCATGTGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTa                                                                                                           <  1:211601/149‑1 (MQ=255)

                                aCTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTAcc                                                                                                         >  3:264301/1‑149 (MQ=255)

                                  tGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGt                                                                                                       >  5:17697/1‑149 (MQ=255)

                                    cTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTc                                                                                                      <  6:17697/148‑1 (MQ=255)

                                      cGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTAc                                                                                                          <  3:84166/142‑1 (MQ=255)

                                      cGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTAc                                                                                                          >  4:84166/1‑142 (MQ=255)

                                       ggCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCt                                                                                                   <  5:224123/148‑1 (MQ=255)

                                          gACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTc                                                                                               <  2:115340/149‑1 (MQ=255)

                                          gACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGAGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTc                                                                                               >  8:228967/1‑149 (MQ=255)

                                          gACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGAGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTc                                                                                               >  8:232605/1‑149 (MQ=255)

                                             cACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCgg                                                                                             >  8:141724/1‑148 (MQ=255)

                                                                             gCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAg                                                            <  7:35292/149‑1 (MQ=255)

                                                                                 aGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGTCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCAt                                                        >  3:234634/1‑149 (MQ=255)

                                                                                         ggCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGa                                                >  4:61883/1‑149 (MQ=255)

                                                                                          ggttctgTCTAATCCTGTCGTTGCTCACCGTTCTTGTGCTGACTTTCGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGTCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGTCCGGCAGCAGTCATTACTTCGAc                                               <  5:211160/143‑1 (MQ=255)

                                                                                                 tCGAATCCTTACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGCAGGGTGAGGCCGTACCGTCGCTGGTCGGCAATCGCATTGCCGTGAGCATGAACGGAACCAGTCATTACTTCTacaacatcg                                        >  6:211160/1‑146 (MQ=255)

                                                                                                 tCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAgg                                        >  5:118255/1‑149 (MQ=255)

                                                                                                    aaTCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTc                                                  >  2:271401/1‑136 (MQ=255)

                                                                                                    aaTCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTc                                                  <  1:271401/136‑1 (MQ=255)

                                                                                                    aaTCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGt                                     >  3:144830/1‑149 (MQ=255)

                                                                                                                    aCCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCaa                     >  2:215700/1‑149 (MQ=255)

                                                                                                                     ccGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAAt                    <  5:13463/149‑1 (MQ=255)

                                                                                                                           tGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGAt               >  7:125944/1‑148 (MQ=255)

                                                                                                                                  gACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGcc       >  5:272521/1‑149 (MQ=255)

                                                                                                                                   aCATTGGAGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGcc       <  7:228967/148‑1 (MQ=255)

                                                                                                                                   aCATTGGAGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGcc       <  7:232605/148‑1 (MQ=255)

                                                                                                                                     aTTGGCGGGCGGCCTGGCCCACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGCCTTg    >  3:247393/1‑149 (MQ=255)

                                                                                                                                     aTTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGCCTTg    >  8:72830/1‑149 (MQ=255)

                                                                                                                                        ggcgggcgGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGCCTTggc  <  7:72830/148‑1 (MQ=255)

AGCGCCGAACCGGCACTGCCCACTCAACGGCGACTGCTCGGCGACCACTACCTCAACGGCGCCGCCTTGGCAGATGAGCAAAGGGCGCAGGCATGAGTCGAATCCTGACGTTGCTCACCGTTCTGGTGCTGACATTGGCGGGCGGCCTGGCACACGCCGAGGAGGAGGGTGAGGCCGTACCGTCGCTGGTCGGCAAGCGCATTGCCGTGAGCATGACCGGCACCAGTCATTACTTCGACATCAAGGCGTTCCAGGCGCAGGTCAATGAGATCAAGCGCCTTGGC  >  NC_002947/3140992‑3141275

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