BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,938,656	Δ1 bp	coding (1008/1443 nt)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,938,656	0	C	.	97.6%	169.8 / ‑3.2	41	coding (1008/1443 nt)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (20/20); minor base G (1/0); total (21/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 = 9.95e-01

ATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGAT  >  NC_000913/1938519‑1938793

aTTACTAACCCGATGCCGACGCGTGCAGACGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCCGCCGAAACTGCCGCTGGGCAGTATCCGTCCGCAACCGTCGCAGCCATAGCGCGGGttgc                                                                                                                                      >  2:194684/1‑141 (MQ=255)

 ttcgtAATCCGAGGCCGAGCGGTGGAGAAGGCATGGAGGGAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACCGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGcc                                                                                                                                    <  2:73276/139‑1 (MQ=255)

     tAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCCTGGCGCG‑GTTTGCCTggc                                                                                                                                >  2:1010456/1‑142 (MQ=255)

         ccGATGCCGACGCGGGCAGAAGTCATGGACGTGGCACACGCCGTGCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCg                                                                                                                            <  1:202949/143‑1 (MQ=255)

                cGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAGACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGaaaaaa                                                                                                                     <  1:334274/143‑1 (MQ=255)

                 gACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCGGTCTGTAGTCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAAt                                                                                                                    >  1:325655/1‑143 (MQ=255)

                       tGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAg                                                                                                              >  2:691005/1‑143 (MQ=255)

                              gTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTTATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAAc                                                                                                       <  2:587649/143‑1 (MQ=255)

                                    gACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAGCTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCt                                                                                                 >  1:527689/1‑143 (MQ=255)

                                      cGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGCGCGGAACAAATCCCGAGCATCAACGTTTCtga                                                                                               >  2:1006775/1‑141 (MQ=255)

                                          gCAAACGCCCTTCTGGCTGGTACTGACGCTGTGATGCTGGCTGCAGAAACTGCCGCTGGGCAGTATCCGTTAGAAAGCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAacac                                                                                           <  1:141564/143‑1 (MQ=255)

                                            aaaCGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCGGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCg                                                                                         >  1:528132/1‑143 (MQ=255)

                                               cGCCGTTCTGGATGGTACTGGCGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCt                                                                                      <  2:940286/143‑1 (MQ=255)

                                                gCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGCGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGAGGAAAAAATCCCGAGCATCAACGTTTCTAAACAGCGtgtg                                                                                     >  2:410595/1‑143 (MQ=255)

                                                    ttCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACg                                                                                 >  1:568233/1‑143 (MQ=255)

                                                        ggATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCa                                                                             <  1:795464/143‑1 (MQ=255)

                                                           tcgTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGtt                                                                          <  1:152867/141‑1 (MQ=255)

                                                              tACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCAGCAACGTTACTAAACACCGTCTGGACGTTCAGGTCGa                                                                       >  1:303350/1‑143 (MQ=255)

                                                              tACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GGTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTCAACCCCGTCAGGACGTTCAGTTCGa                                                                       >  2:271053/1‑143 (MQ=255)

                                                                   aCGCGGTGATGCTGGCTGCAGGAACTGCCGCGGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAAtg                                                                  <  2:635062/143‑1 (MQ=255)

                                                                                  cggCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGCAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTg                                                   <  1:351041/141‑1 (MQ=255)

                                                                                  cTGCAGAAGCGGCCGCTGGGCAGTAGCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCGTGGGTGCGGAACAAACTCGGATCATCAACGTTTCGAAACACCGTCTGGACGGTCAGTTCGACAATGTGGAAGAAGCTATTg                                                   <  2:167937/143‑1 (MQ=255)

                                                                                   tGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTCTTGc                                                  >  2:799947/1‑143 (MQ=255)

                                                                                   tGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGc                                                  <  2:528132/143‑1 (MQ=255)

                                                                                       gAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCaag                                              >  1:789544/1‑141 (MQ=255)

                                                                                              ccGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAACAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGGCAGGTCAGCaa                                       >  2:195664/1‑143 (MQ=255)

                                                                                               cGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAAt                                      >  1:871610/1‑143 (MQ=255)

                                                                                                 cTGGGCAGTATCCGACAGAAACCCTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAACAATCCCGAGCATCAACGTTTCTGAGCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCCATTGCCATGTCAGCagcgt                                    >  2:537868/1‑139 (MQ=255)

                                                                                                   gggCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTAc                                  >  1:858059/1‑143 (MQ=255)

                                                                                                   gggCAGTATCCGTCAGAAACCGTGGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTAc                                  <  1:526064/143‑1 (MQ=255)

                                                                                                    ggCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACg                                 <  1:127821/143‑1 (MQ=255)

                                                                                                      cAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGCAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCACTGTACGCa                               >  2:883137/1‑143 (MQ=255)

                                                                                                           tCCGTCAGAAACCGTTGCAGCCATCGCGCG‑GTTTGCCTTGTTTCTGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTaa                          <  1:764697/143‑1 (MQ=255)

                                                                                                              gtgaGAAAGCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                       <  1:1005969/140‑1 (MQ=255)

                                                                                                              gTCAGAAACCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                       >  1:946473/1‑143 (MQ=255)

                                                                                                                   aaaCCGTGGCAGCCATGGCGCG‑GTTGGCCGGGGTGCGGAAGAAATCCCGAGCATCACCGTTTCTAAACACCGTCTGGCCGTTCAGTCCGACAATGTGGAAGCAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGa                  <  2:916732/143‑1 (MQ=255)

                                                                                                                     aCCGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAAACCCGAGCCTCAATGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGaaa                <  1:847843/143‑1 (MQ=255)

                                                                                                                       cggtggaGCCGTGGCGCG‑GTTGACGTGGGGGCGGAAGAGATCCCGAGCATCAACGTTTCTAGACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGTTAACCACCTGAAAgg              <  2:856987/137‑1 (MQ=255)

                                                                                                                       cGTTGCAGCCATGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAgg              >  1:941141/1‑143 (MQ=255)

                                                                                                                                 aTGGCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCGTGAACGATTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGCAAGGCATTACGGCg    >  2:517796/1‑143 (MQ=255)

                                                                                                                                   ggCGCG‑GTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGat  <  1:664445/143‑1 (MQ=255)

ATTACTAACCCGATGCCGACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGAT  >  NC_000913/1938519‑1938793

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