BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,656,933	A→G	9.5%	C51R (TGC→CGC)	fdx ←	[2Fe‑2S] ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,656,933	0	A	G	9.5%	116.9 / 3.5	42	C51R (TGC→CGC)	fdx	[2Fe‑2S] ferredoxin
Reads supporting (aligned to +/- strand): ref base A (20/18); new base G (2/2); total (22/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 = 1.00e+00

CTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATAACAATCTTTGGC  >  NC_000913/2656792‑2657081

ctactaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAg                                                                                                                                                     <  1:7704/143‑1 (MQ=255)

ctactaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAg                                                                                                                                                     >  2:7704/1‑143 (MQ=255)

   ctaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTGACGAACGATGCAGTGGCAGGTGGTGCaa                                                                                                                                            <  6:116482/149‑1 (MQ=255)

            cgtcgGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGGCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGcaca                                                                                                                                       <  8:76878/145‑1 (MQ=255)

                ggTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt                                                                                                                               <  8:147391/149‑1 (MQ=255)

                 gTAACGCGCGCCTGGCAGCTGAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAACACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTc                                                                                                                              <  4:222835/149‑1 (MQ=255)

                   aaCGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGCGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTca                                                                                                                             >  5:144818/1‑148 (MQ=255)

                                  gCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGGCTTCCTGCTCTGAGCTTTCGGGCAGTGAGTCAATACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAAtct                                                                                                             <  5:200001/149‑1 (MQ=255)

                                   cTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAAtctc                                                                                                            <  5:190561/149‑1 (MQ=255)

                                    ttAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAAtctc                                                                                                            >  7:14587/1‑148 (MQ=255)

                                       aaCGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATa                                                                                                        >  8:65761/1‑149 (MQ=255)

                                          ggCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCAAGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCTATACCg                                                                                                     >  5:103413/1‑149 (MQ=255)

                                             tttCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGa                                                                                                          >  2:102138/1‑141 (MQ=255)

                                             tttCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGa                                                                                                          <  1:102138/141‑1 (MQ=255)

                                                   gCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCa                                                                                              >  3:34276/1‑147 (MQ=255)

                                                         gTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCAGGCAGTGAGTCAAAACATTCACGAAAGATGCATTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCg                                                                                      >  1:144665/1‑149 (MQ=255)

                                                         gTCCCCAGGCTTTGTCCAGCATGTCGGCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCg                                                                                      <  1:75087/149‑1 (MQ=255)

                                                            cccAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCg                                                                                   <  5:115785/149‑1 (MQ=255)

                                                               aGGCTTTGTCCAGCATTTCGTCTTCGTGCTCTGAGCTTTCCGGCAGTGAGGCAAAACCTTCACGAACGATGCACTGGCAGGGGGTGCTAGCACAGGATTTTTCACAGGCGTGGTCAATCTCGATACCGTTACGCAGAGCTGCGTCgaga                                                                                <  6:144818/149‑1 (MQ=255)

                                                               aGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCgag                                                                                 >  7:123189/1‑148 (MQ=255)

                                                                ggCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGaa                                                                               <  7:65761/149‑1 (MQ=255)

                                                                       tCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGttt                                                                         <  8:123189/148‑1 (MQ=255)

                                                                         cAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCAc                                                                      <  4:34276/149‑1 (MQ=255)

                                                                          aGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTCGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCAcc                                                                     >  3:152984/1‑149 (MQ=255)

                                                                                tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTAt                                                                <  8:14587/148‑1 (MQ=255)

                                                                                  gtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAg                                                             >  4:169460/1‑149 (MQ=255)

                                                                                              tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGc                                                                                       >  4:11907/1‑111 (MQ=255)

                                                                                              tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGc                                                                                       <  3:11907/111‑1 (MQ=255)

                                                                                              tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCGTCCAGAACAg                                                 >  4:46400/1‑149 (MQ=255)

                                                                                                   tttCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCGGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCa                                            >  2:28510/1‑149 (MQ=255)

                                                                                                   tttCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCGGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCa                                            >  2:24751/1‑149 (MQ=255)

                                                                                                        ggCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAgg                                       >  7:205918/1‑149 (MQ=255)

                                                                                                           aGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAggg                                      >  7:194618/1‑147 (MQ=255)

                                                                                                                gTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCa                                    >  5:202227/1‑144 (MQ=255)

                                                                                                                   aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATg                                                                             <  4:13142/100‑1 (MQ=255)

                                                                                                                   aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATg                                                                             >  3:13142/1‑100 (MQ=255)

                                                                                                                    aaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCt                            >  8:86485/1‑148 (MQ=255)

                                                                                                                           cACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAg                      <  3:169460/147‑1 (MQ=255)

                                                                                                                              gaacgaTGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCaa                  >  4:178894/1‑148 (MQ=255)

                                                                                                                                 cgaTGCAGTGGCGGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATa              <  1:28510/149‑1 (MQ=255)

                                                                                                                                 cgaTGCAGTGGCGGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATa              <  1:24751/149‑1 (MQ=255)

                                                                                                                                         tGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATAACAATCtt      >  6:54106/1‑149 (MQ=255)

                                                                                                                                             aGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATAACAATCTTTgga  <  8:205918/149‑2 (MQ=255)

CTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGATACCGTTACGCAGAGCTGCGTCGAGAATGGTTTCACCGCTATTAGCTTCCAGAACAGCGCCATCAGGGCAGAGATCCTGATGAGGCAAAATAACAATCTTTGGC  >  NC_000913/2656792‑2657081

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