BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,986,394	G→A	11.8%	G324S (GGC→AGC)	cpsG →	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,986,394	0	G	A	11.8%	76.4 / 3.9	34	G324S (GGC→AGC)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand): ref base G (14/16); new base A (2/2); total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

CCCAGGCAGCAAAATCATTCACGACCCGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAAAAGCCTGGGGCAACTGGTGGACGAACGTATTGCAGC  >  NC_002947/1986249‑1986530

cccAGGCAGCAAAATCATTCACGACCCGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGa                                                                                                                                      >  1:485220/1‑150 (MQ=255)

    ggCAGCAAAATCATTCACGACCCGCGCCTGCCCTCGCA‑ACCAACACACCAGGGTGCAGAAGCAAGTGGCACTTCAAGCCAGAGCAAGACCGGCCATGACTGCATAAAAAAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATg                                                                                                                                  >  1:233141/1‑150 (MQ=255)

       agcaAAATCATTCACGACCCGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGTGGGAATGTCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAgc                                                                                                                               <  1:844202/150‑1 (MQ=255)

               cATTCACGACCCGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCTAAGATGCTGTATATGGCGGCGAAATGAGCGCccacca                                                                                                                       <  2:389359/150‑1 (MQ=255)

                         ccGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGt                                                                                                             >  2:252944/1‑150 (MQ=255)

                         ccGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGt                                                                                                             >  1:773723/1‑150 (MQ=255)

                                tGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTg                                                                                                      <  1:215883/150‑1 (MQ=255)

                                             cATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACtt                                                                                                        <  1:136120/135‑1 (MQ=255)

                                             cATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACtt                                                                                                        >  2:136120/1‑135 (MQ=255)

                                              aTCA‑ACCAGCTCGCAGAAGCAGGGGGGATGGCAAGCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTGTGCGCGTCGAAGATGCTTTATATTGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACacc                                                                                        >  2:270621/1‑148 (MQ=255)

                                               acA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCg                                                                                       <  2:509290/149‑1 (MQ=255)

                                                           gaagaaGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCt                                                                            <  1:434004/150‑1 (MQ=255)

                                                           gaagaaGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGGATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCt                                                                            >  1:218940/1‑150 (MQ=255)

                                                              gaagCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTgg                                                                         <  2:656650/150‑1 (MQ=255)

                                                                         aaTGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCgg                                                                                      >  1:238029/1‑126 (MQ=255)

                                                                         aaTGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCgg                                                                                      <  2:238029/126‑1 (MQ=255)

                                                                                  ccAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTtgat                                                     <  1:252944/150‑1 (MQ=255)

                                                                                    agagCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTtgatga                                                   >  1:261281/1‑150 (MQ=255)

                                                                                      agCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCAGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTc                                                                      <  1:248401/129‑1 (MQ=255)

                                                                                      agCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCAGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTc                                                                      <  1:248314/129‑1 (MQ=255)

                                                                                      agCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCAGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTc                                                                      >  2:248314/1‑129 (MQ=255)

                                                                                      agCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCAGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTc                                                                      >  2:248401/1‑129 (MQ=255)

                                                                                        cAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGCCGATGTTGAAAATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGt                                               >  1:150267/1‑150 (MQ=255)

                                                                                              cGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCCCCTGGTCGCTCAGTTGATGAGCGTAACggg                                         <  2:191745/150‑1 (MQ=255)

                                                                                              cGGCCATATCTCCGTAAAAGAGCGTATGCGCGGCGACGAGCCTGTACATGGCGGCGAAATGAACGCCCACCATTACTTCTGCGACTTTGCATACTGCGACAGCGGTAACATCCCTAGACTCCTGGTCGCTAAGTTGATGAGCGTAACggg                                         >  2:564484/1‑150 (MQ=255)

                                                                                                gCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCa                                       >  2:257439/1‑150 (MQ=255)

                                                                                                   aTGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCaaaa                                    <  1:257439/150‑1 (MQ=255)

                                                                                                       cTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAACAGCCt                                <  2:673646/150‑1 (MQ=255)

                                                                                                       cTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCCCAGTTGATGAGCGTAACGGGCAAAAGCCt                                >  1:537089/1‑150 (MQ=255)

                                                                                                        ttCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAAAAGCCTg                               <  2:773723/150‑1 (MQ=255)

                                                                                                                      tATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAAAAGCCTGGGGCAACTGGTGGa                 <  2:261281/150‑1 (MQ=255)

                                                                                                                                     tGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGATAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAAAAGCCTGGGGCAACTGGTGGACGAACGTATTGCAGc  <  1:481484/150‑1 (MQ=255)

                                                                                                                                        tGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCa                                                           >  2:320652/1‑90 (MQ=255)

                                                                                                                                        tGTATATGGCGGCGAAATGAGCGCCCACCATTACGTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCa                                                           <  1:320652/90‑1 (MQ=255)

CCCAGGCAGCAAAATCATTCACGACCCGCGCCTGACCTGGAACACCATCA‑AACAGGTCGAAGAAGCAGGGGGAATGGCAATCCAGAGCAAGACCGGCCATGCCTTCATAAAAGAGCGTATGCGCGCCGAAGATGCTGTATATGGCGGCGAAATGAGCGCCCACCATTACTTCCGTGACTTTGCCTACTGCGACAGCGGTAACATCCCTTGGCTCCTGGTCGCTCAGTTGATGAGCGTAACGGGCAAAAGCCTGGGGCAACTGGTGGACGAACGTATTGCAGC  >  NC_002947/1986249‑1986530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`