D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,453,525	A→G	V281A (GTT→GCT)	SAUSA300_RS07200 ←	virulence factor C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,453,525	0	A	G	98.2%	141.8 / ‑4.8	56	V281A (GTT→GCT)	SAUSA300_RS07200	virulence factor C
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (4/51); total (4/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

GGATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAACAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTAACCGGTA  >  NC_007793/1453441‑1453621

ggATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCAc                                                                                   <  1:2421291/100‑1 (MQ=255)

ggATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCAc                                                                                   <  1:2235152/100‑1 (MQ=255)

 gATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCATGACTTTTATCACg                                                                                  <  1:5422372/100‑1 (MQ=255)

 gATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACg                                                                                  <  1:5346579/100‑1 (MQ=255)

 gATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACg                                                                                  <  1:3567837/100‑1 (MQ=255)

 gATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACg                                                                                  <  1:3007367/100‑1 (MQ=255)

 gATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACg                                                                                  <  1:4925919/100‑1 (MQ=255)

  atcatcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCaagg                                                                                 <  1:2895143/100‑4 (MQ=255)

  atcatcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACgg                                                                                 <  1:3142545/100‑2 (MQ=255)

  atcatcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACgg                                                                                 <  1:224014/100‑2 (MQ=255)

  atcatcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGa                                                                                 <  1:1180120/100‑1 (MQ=255)

   tcatcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGaa                                                                                <  1:2049435/100‑1 (MQ=255)

    catcTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGaat                                                                               <  1:805532/100‑2 (MQ=255)

          aTAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCtt                                                                         <  1:3983869/100‑1 (MQ=255)

          aTAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCta                                                                         <  1:2189961/100‑2 (MQ=255)

            agtagCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGGCTTTTATCACGAAGCCCCtttt                                                                       <  1:190639/100‑3 (MQ=255)

            agtagCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTtat                                                                       <  1:348468/100‑1 (MQ=255)

              tagCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTtatat                                                                     <  1:4625725/100‑1 (MQ=255)

               agCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTtatata                                                                    <  1:4311601/100‑1 (MQ=255)

                 cACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTtataagca                                                                  <  1:4076540/100‑5 (MQ=255)

                  aCCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTtatatacac                                                                 <  1:4178195/100‑4 (MQ=255)

                    cATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATaagcac                                                               <  1:997226/100‑5 (MQ=255)

                       ttCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGACCCTTATATAAATACGGt                                                            <  1:1259367/100‑1 (MQ=255)

                       ttCTGGCAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTGCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGt                                                            <  1:3740442/100‑1 (MQ=255)

                        tCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGtt                                                           <  1:3300106/100‑2 (MQ=255)

                        tCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTGCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTa                                                           <  1:4696876/100‑1 (MQ=255)

                            gTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTaactg                                                       <  1:1267101/100‑4 (MQ=255)

                                 gCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCtt                                                  <  1:820459/100‑1 (MQ=255)

                                  ccTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCttt                                                 <  1:5447972/100‑1 (MQ=255)

                                       ggATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACttt                                            >  1:3578871/1‑100 (MQ=255)

                                         ttgcccTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCa                                          <  1:1940281/97‑1 (MQ=255)

                                              cTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCAt                                     <  1:4369812/100‑1 (MQ=255)

                                                aaaTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATAc                                   >  1:4830602/1‑100 (MQ=255)

                                                 aaTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATAcc                                  <  1:2221058/100‑1 (MQ=255)

                                                 aaTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATAcc                                  <  1:939426/100‑1 (MQ=255)

                                                  aTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCt                                 <  1:3613251/100‑1 (MQ=255)

                                                  aTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCt                                 <  1:1198801/100‑1 (MQ=255)

                                                  aTCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCt                                 <  1:1591948/100‑1 (MQ=255)

                                                   tCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATAccca                                <  1:2853356/100‑3 (MQ=255)

                                                   tCGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCta                                <  1:58800/100‑1 (MQ=255)

                                                    cGCTTATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCtaa                               <  1:4106574/100‑1 (MQ=255)

                                                       ctATGCAATCCCCTGCTGTTCTTATTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCtaataa                            >  1:4126635/2‑100 (MQ=255)

                                                        tATGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTaataat                           >  1:3352859/1‑100 (MQ=255)

                                                          tGCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATAc                         <  1:3942984/100‑1 (MQ=255)

                                                           gCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACa                        <  1:5213329/100‑1 (MQ=255)

                                                           gCAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACa                        <  1:64654/100‑1 (MQ=255)

                                                            cAATCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACac                       <  1:5319574/100‑2 (MQ=255)

                                                             aaTCCCCTGCTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACaac                      <  1:1953718/100‑2 (MQ=255)

                                                             aaTCCCCTACTGTTCTTCTTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACaag                      <  1:5277174/100‑2 (MQ=255)

                                                              aTCCCCTGCTGTTCTTCTTACAACAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATc                     <  1:1637224/100‑1 (MQ=255)

                                                                          tcttctTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTa         <  1:634021/100‑1 (MQ=255)

                                                                          tcttctTACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTa         <  1:4369047/100‑1 (MQ=255)

                                                                           cttcttACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTag        <  1:348515/100‑2 (MQ=255)

                                                                            ttcttACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTagc       <  1:3136785/100‑3 (MQ=255)

                                                                            ttcttACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTagc       <  1:4539994/100‑3 (MQ=255)

                                                                                tACAGCAGGACTTTTATCACGAAGCCCCTTATATAAATGCGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTAACCGGt   <  1:3236596/100‑1 (MQ=255)

                                                                                 acagcaGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTagatcgga  >  1:3336576/1‑93 (MQ=255)

                                                                                 acagcaGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTAACCggcg  <  1:2411643/100‑3 (MQ=255)

GGATCATCTAATAGTAGCACCATTTCTGGTAGTGCCTCTGGATACCCTAAATCGCTTATGCAATCCCCTGCTGTTCTTCTTACAACAGGACTTTTATCACGAAGCCCCTTATATAAATACGGTAAAATTTCTTTACTTTCAATCATACCTAATAATACAATCGCTTGACGTCTAACCGGTA  >  NC_007793/1453441‑1453621

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