BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,990,626	A→T	15.4%	V169D (GTC→GAC)	PP_4399 ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,990,626	0	A	T	15.4%	83.0 / 9.2	39	V169D (GTC→GAC)	PP_4399	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base A (16/17); new base T (4/2); total (20/19)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CGGGCGTACGCGCCAGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCATCGCT  >  NC_002947/4990480‑4990765

cGGGCGTACGCGCCAGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGaca                                                                                                                                           <  4:374294/149‑1 (MQ=255)

 gggCGTACGCGCCAGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGCTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGacac                                                                                                                                          <  2:344515/149‑1 (MQ=255)

          cgcCAGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGt                                                                                                                                 <  1:81163/149‑1 (MQ=255)

           gcCAGCGGCTTGGCCGTGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTg                                                                                                                                >  1:264581/1‑149 (MQ=255)

              aGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCAt                                                                                                                             >  2:262870/1‑149 (MQ=255)

               gCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCAt                                                                                                                             >  5:163530/1‑148 (MQ=255)

                 tggttGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGTGCAGGTGCTTGAAGAGGATGGTCTGGTCGGCCAACAGCAGCGGGTGGCCGTGCTGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGaa                                                                                                                          <  3:235883/146‑1 (MQ=255)

                  gCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAg                                                                                                                         <  6:163530/149‑1 (MQ=255)

                  gCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAg                                                                                                                         <  2:312924/149‑1 (MQ=255)

                      ggCCATGCAGACCACGGTCAGTCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAAGAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGgt                                                                                                                       <  8:172528/147‑1 (MQ=255)

                      ggCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGgt                                                                                                                       <  1:19243/147‑1 (MQ=255)

                        ccATGCAGACCACGGTCAGCCCGTGCTTGAGCGACACTTCCTTGAGCAGGTGCTTTAACAGGACGGTCTGGTCGGCCAACAGAAGCGGGTCGCCTTGCAGCAGGTTGATCTCGACCTGTCTGACACCCATTTCGATCATGAAGGTGTCa                                                                                                                   >  8:298682/1‑149 (MQ=255)

                             cAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTgg                                                                                                              <  3:343479/149‑1 (MQ=255)

                                 ccACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAg                                                                                                           <  2:201469/148‑1 (MQ=255)

                                  cACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGc                                                                                                         >  3:205736/1‑149 (MQ=255)

                                       tCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAgc                                                                                                    <  2:264581/149‑1 (MQ=255)

                                          gCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGcgcg                                                                                                 >  6:188928/1‑149 (MQ=255)

                                                 cTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATg                                                                                           <  5:188928/148‑1 (MQ=255)

                                                        gcCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGTCATGCAGTGGTa                                                                                   <  1:262870/149‑1 (MQ=255)

                                                                 ttGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCtt                                                                          >  8:32213/1‑149 (MQ=255)

                                                                    aGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTga                                                                        <  7:236314/148‑1 (MQ=255)

                                                                    aGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCATGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTg                                                                         >  1:295482/1‑147 (MQ=255)

                                                                         gTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGTTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACg                                                                  <  7:298682/149‑1 (MQ=255)

                                                                                     gAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCgt                                                      >  6:190078/1‑149 (MQ=255)

                                                                                        ggTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCgttgt                                                   <  4:168758/149‑1 (MQ=255)

                                                                                         gTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCgttgtt                                                  <  7:32213/149‑1 (MQ=255)

                                                                                            tGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACGTCCTTGAAGAACGGGCGCAGGCCGTTGTTgg                                                >  4:367941/1‑148 (MQ=255)

                                                                                             ggTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGaa                                              >  8:356236/1‑149 (MQ=255)

                                                                                             ggTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGaa                                              >  6:159875/1‑149 (MQ=255)

                                                                                              gTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGaa                                              <  4:205736/148‑1 (MQ=255)

                                                                                               tCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACt                                            >  3:279489/1‑149 (MQ=255)

                                                                                                                  cGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGACGGTGGCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATgg                         <  2:184071/149‑1 (MQ=255)

                                                                                                                               gTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCgcggc            >  5:6817/1‑149 (MQ=255)

                                                                                                                               gTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCgcggc            >  5:1385/1‑149 (MQ=255)

                                                                                                                                  atgcacGAAATGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACt                                            <  2:200090/109‑1 (MQ=38)

                                                                                                                                  atgcacGAAATGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACt                                            >  1:200090/6‑114 (MQ=38)

                                                                                                                                    gcacGAAATGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCa       >  5:286264/4‑149 (MQ=255)

                                                                                                                                     cacGAAATGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCAt      >  8:216686/3‑149 (MQ=255)

                                                                                                                                      acGAAATGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCATc     >  1:307004/2‑149 (MQ=255)

                                                                                                                                        gAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCATCGc   >  7:177622/1‑149 (MQ=255)

                                                                                                                                         aaatGGGTGTCAGCCAGTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCATCGCt  <  6:286264/146‑1 (MQ=255)

CGGGCGTACGCGCCAGCGGCTTGGCCATGCAGACCACGGTCAGCCCGTGCTTGAGCGCCACTTCCTTGAGCAGGTGCTTGAACAGGAAGGTCTGGTCGGCCAACAGCAGCGGGTCGCCGTGCAGCAGGTTGATCTCGAACTGGCTGACACCCATTTCGTGCATGAAGGTGTCACGTGGCAGGCCGAGCGCGGCCATGCAGTGGTACACCTCCTTGAAGAACGGGCGCAGGCCGTTGTTGGAACTGACACTGAACGCCGAATGGCCCAGCTCGCGGCGGCCATCGCT  >  NC_002947/4990480‑4990765

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