D3A :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,589,898	A→G	T10A (ACT→GCT)	SAUSA300_RS07845 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,589,898	0	A	G	96.3%	124.2 / ‑4.7	55	T10A (ACT→GCT)	SAUSA300_RS07845	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (46/6); total (49/6)
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

AAAAAAATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGA  >  NC_007793/1589817‑1589992

ttactaaTCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCaa                                                                              >  1:2422635/6‑99 (MQ=21)

     aaTCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAc                                                                         >  1:3584357/1‑100 (MQ=25)

                   tacaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGaaaa                                                           >  1:2046978/1‑100 (MQ=38)

                    acaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACCAAAGTTTTGAAAAg                                                          >  1:4948690/4‑100 (MQ=38)

                     caaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgc                                                         >  1:3996684/3‑99 (MQ=38)

                     aaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAAAAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAgc                                                         >  1:298778/3‑99 (MQ=38)

                         tgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAAAAAAGTTTTGCAAAgccggg                                                     >  1:3375461/1‑95 (MQ=255)

                         cgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgccggg                                                     >  1:4670977/2‑95 (MQ=255)

                               tccgatctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAg                                               <  1:4220980/93‑1 (MQ=255)

                                 ctGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCg                                             >  1:3347942/2‑100 (MQ=255)

                                 ccgatctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTGGAAAAGCCGGGTGCCAGCg                                             <  1:3142452/93‑1 (MQ=255)

                                  tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGa                                            >  1:1680375/1‑100 (MQ=255)

                                  tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGTTTTGAAAAGCCGGGCGCCAGCGa                                            >  1:1104808/1‑100 (MQ=255)

                                   gAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaa                                           >  1:5777351/1‑100 (MQ=255)

                                    tctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                          >  1:3452206/3‑100 (MQ=255)

                                    taTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                          >  1:5444817/2‑100 (MQ=255)

                                    gctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                                          >  1:3312694/3‑100 (MQ=255)

                                    gaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTATTTGTAGCAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaag                                          >  1:953029/2‑99 (MQ=255)

                                    aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                                          >  1:828730/1‑100 (MQ=255)

                                    aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                          >  1:3914006/1‑100 (MQ=255)

                                     ttggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAGAGCCGGGTGCCAGCGaaaa                                         >  1:2617275/2‑100 (MQ=255)

                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                                         >  1:4588470/1‑100 (MQ=255)

                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaa                                         >  1:4740688/1‑100 (MQ=255)

                                     aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAAAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                                         >  1:626814/1‑100 (MQ=255)

                                      gtgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                                        >  1:2295226/3‑99 (MQ=255)

                                      gcctggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                                        >  1:2287713/4‑99 (MQ=255)

                                      ctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAACTACTCGTAATAACCCAATAAACAAAAGCTTTGAAAAACCGGGGGGCAGCGaagat                                        >  1:445744/3‑97 (MQ=255)

                                      acatggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAACTACTTGTAGTAAACCAATAAACAAAAGGGTTGAAAAGCCGGGGGCCAGCGgaaac                                        >  1:6321703/4‑99 (MQ=255)

                                       ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTTAAAAGCCGGGTGCCAGCGaaaact                                       >  1:5815377/1‑98 (MQ=255)

                                       ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGaaaact                                       >  1:5568292/1‑98 (MQ=255)

                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAGCCACTCGTAGTAATCCAATAAACAAAAGCTTTGAAAAGCCCGGCGCCAgatcggaaga                                      >  1:496662/1‑91 (MQ=255)

                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACtt                                      >  1:2350372/1‑100 (MQ=255)

                                        gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGATACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACtt                                      >  1:4786197/1‑100 (MQ=255)

                                          ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGAAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTaa                                    >  1:5107228/1‑100 (MQ=255)

                                          ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAACCTCCTCGTTGTAACCCCATTAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACTTaa                                    >  1:3111504/1‑100 (MQ=255)

                                                    ttATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTAAAAAGCACttt                          >  1:3428327/1‑100 (MQ=255)

                                                    ctATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACttt                          >  1:3021768/2‑100 (MQ=255)

                                                     tgtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                         >  1:1009260/8‑100 (MQ=255)

                                                     tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                         >  1:4147905/1‑100 (MQ=255)

                                                     tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                         >  1:4847659/1‑100 (MQ=255)

                                                      gtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:1486363/7‑100 (MQ=255)

                                                      ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:3590381/7‑100 (MQ=255)

                                                      ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:4542906/7‑100 (MQ=255)

                                                      atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:5030204/7‑100 (MQ=255)

                                                      atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:847219/7‑100 (MQ=255)

                                                      atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                        >  1:6167256/7‑100 (MQ=255)

                                                       tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                       <  1:2546694/95‑1 (MQ=255)

                                                       tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                       >  1:6307574/6‑100 (MQ=255)

                                                       tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGATTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                       >  1:5653315/6‑100 (MQ=255)

                                                       tGGCAGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                       <  1:1395034/100‑1 (MQ=255)

                                                        ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                      >  1:227284/5‑100 (MQ=255)

                                                        ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                      >  1:141162/5‑100 (MQ=255)

                                                                  aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTa            <  1:6159713/100‑1 (MQ=255)

                                                                  aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTCTATCAGAAAAAGCTa            <  1:4160601/100‑1 (MQ=255)

                                                                       aaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAGAGCTAAGatc       >  1:5303032/1‑98 (MQ=255)

                                                                            aagaagCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGa  >  1:365721/1‑100 (MQ=255)

AAAAAAATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGA  >  NC_007793/1589817‑1589992

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`