BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,892,993	G→A	E113K (GAA→AAA)	ygcN →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,892,993	0	G	A	90.9%	127.5 / 3.6	45	E113K (GAA→AAA)	ygcN	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (3/1); new base A (18/22); total (22/23)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGGCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGAT  >  NC_000913/2892898‑2893085

ccGATGGCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAaaaa                                                                                          <  2:582295/100‑1 (MQ=255)

 cGATGGCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCaaaaaa                                                                                         <  1:73967/100‑1 (MQ=255)

     ggCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAaagaag                                                                                     <  1:203188/100‑1 (MQ=255)

     ggCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAaagaag                                                                                     <  1:83498/100‑1 (MQ=255)

     ggCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAaagaag                                                                                     <  1:543652/100‑1 (MQ=255)

      gCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAgg                                                                                    <  1:802748/100‑1 (MQ=255)

           acgacgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCg                                                                               <  2:839415/100‑1 (MQ=255)

           acgacgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCg                                                                               <  2:1150020/100‑1 (MQ=255)

            cgacgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGa                                                                              >  2:79323/1‑100 (MQ=255)

              acgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGaa                                                                             <  2:39513/99‑1 (MQ=255)

               cgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATg                                                                           >  2:633094/1‑100 (MQ=255)

               cgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATg                                                                           >  2:632079/1‑100 (MQ=255)

               cgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATg                                                                           <  1:1005684/100‑1 (MQ=255)

               cgTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATg                                                                           <  1:764057/100‑1 (MQ=255)

                gTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGc                                                                          >  2:292029/1‑100 (MQ=255)

                          ttACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAg                                                                >  2:999159/1‑100 (MQ=255)

                           tACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGc                                                               <  1:478802/100‑1 (MQ=255)

                           tACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGc                                                               >  1:483864/1‑100 (MQ=255)

                             cAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGa                                                             >  2:1176884/1‑100 (MQ=255)

                              aGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGAc                                                            <  1:1244377/100‑1 (MQ=255)

                               gCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACg                                                           >  1:966804/1‑100 (MQ=255)

                                   ggcgtTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTgg                                                       <  2:346870/100‑1 (MQ=255)

                                       gTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGc                                                   <  2:305739/100‑1 (MQ=255)

                                       gTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGc                                                   <  2:483860/100‑1 (MQ=255)

                                           aTCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCCAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTg                                               >  1:132134/1‑100 (MQ=255)

                                             ccTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTa                                             >  2:1285414/1‑100 (MQ=255)

                                              cTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTAt                                            <  2:665882/100‑1 (MQ=255)

                                              cTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTAt                                            >  2:521412/1‑100 (MQ=255)

                                              cTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGGTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTAt                                            >  1:665156/1‑100 (MQ=255)

                                                  aGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATgaag                                        >  2:264885/1‑100 (MQ=255)

                                                  aGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATgaag                                        <  2:723595/100‑1 (MQ=255)

                                                     gtAATACGTGCACGTTTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATgaagaa                                      <  2:221222/99‑1 (MQ=255)

                                                         tACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACgg                                 >  1:343902/1‑100 (MQ=255)

                                                          aCGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGc                                >  1:844809/1‑100 (MQ=255)

                                                          aCGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGc                                <  2:217055/100‑1 (MQ=255)

                                                              gCACGATTCGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCagag                            <  1:958854/100‑1 (MQ=255)

                                                                  gATTCGTTCCGTGGCTGGTTGCCGAAGGCGAAAAAGAAGGTGGCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTc                          <  2:431171/98‑1 (MQ=255)

                                                                   aTTCGATCCGTGGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCtgt                       >  1:516606/1‑100 (MQ=255)

                                                                      cGATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCtgtgac                    >  1:382076/1‑98 (MQ=255)

                                                                       gATCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCg                   >  1:581359/1‑100 (MQ=255)

                                                                         tCCGTGGCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGtt                 <  1:959787/100‑1 (MQ=255)

                                                                              ggCTGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTtg            >  1:212951/1‑100 (MQ=255)

                                                                                 tGGTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTgtgg         <  2:599991/100‑1 (MQ=255)

                                                                                  ggTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTgtggt        >  1:124946/1‑100 (MQ=255)

                                                                                   gTTGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCG‑TATTTGTGGTGa      >  1:509966/1‑100 (MQ=255)

                                                                                     tGCCGAAGCCAAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGAc     >  2:1172033/1‑100 (MQ=255)

                                                                                      gccgaagccaaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACg    >  2:676958/1‑100 (MQ=255)

                                                                                       ccgaagccaaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGa   <  1:800615/100‑1 (MQ=255)

                                                                                       ccgaagccaaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGa   >  2:840486/1‑100 (MQ=255)

                                                                                        cgaagccaaaAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGat  <  1:292032/100‑1 (MQ=255)

CCGATGGCGTAACGACGTTTTCCAGCTTACAGCCCGGCGGTGAATCCTGGAGTGTATTACGTGCACGATTCGATCCGTGGCTGGTTGCCGAAGCCGAAAAAGAAGGTGTCGAATGCATCCCCGGAGCGACGGTGGATGCACTGTATGAAGAAAACGGCAGAGTCTGTGGCGTTATTTGTGGTGACGAT  >  NC_000913/2892898‑2893085

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