D3A :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,589,939	A→G	S23P (TCA→CCG)	SAUSA300_RS07845 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,589,939	0	A	G	98.2%	118.3 / ‑5.6	54	S23S (TCA→TCG)	SAUSA300_RS07845	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (44/9); minor base C (1/0); total (45/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAAT  >  NC_007793/1589848‑1590023

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:4542906/7‑100 (MQ=255)

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

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

                       atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                                       >  1:5030204/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)

                                                tccgatctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTa                              <  1:1631534/94‑1 (MQ=255)

                                                  gctactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat                            >  1:2743385/2‑98 (MQ=255)

                                                  agtactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat                            >  1:5917729/3‑98 (MQ=255)

                                                     actCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATatt                         >  1:2700871/1‑100 (MQ=255)

                                                      ctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATattc                        <  1:1716201/100‑1 (MQ=255)

                                                                      aaaCAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATaaa        >  1:3159360/1‑100 (MQ=255)

                                                                            aaGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCaag  <  1:1511918/100‑2 (MQ=255)

CGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAAT  >  NC_007793/1589848‑1590023

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