Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 129,796 G→C intergenic (‑173/+212) spa ← / ← sarH1 immunoglobulin G binding protein A/staphylococcal accessory regulator H1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730129,7960GC97.1% 101.7 / NA 34intergenic (‑173/+212)spa/sarH1immunoglobulin G binding protein A/staphylococcal accessory regulator H1
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (13/20);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATATAGAAAATAGAAAGT  >  CP000730/129660‑129924
                                                                                                                                        |                                                                                                                                
tatTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACatat                                                                                                                              >  1:677401/1‑141 (MQ=255)
        aTTAAATTTAAAGGTGAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTtata                                                                                                                      >  1:884606/1‑141 (MQ=255)
        aTTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAAAGATACTTTTtcta                                                                                                                      >  1:895484/1‑141 (MQ=255)
        aTTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAAGATAGACATATTTTTtata                                                                                                                      >  1:185554/1‑141 (MQ=255)
          tAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAAt                                                                                                                    <  1:228032/141‑1 (MQ=255)
           aaaTTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATg                                                                                                                   <  1:715398/141‑1 (MQ=255)
           aaaTTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATg                                                                                                                   <  1:641258/141‑1 (MQ=255)
            aaTTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGa                                                                                                                  <  1:1058860/141‑1 (MQ=255)
             aTTTAAAGGGAAGTTTTACAACTTATAAAATACATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGaa                                                                                                                 <  1:1039322/141‑1 (MQ=255)
                tAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAgg                                                                                                              >  1:494032/1‑140 (MQ=255)
                      tAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATAAAGACATATTTTTTATAATGAAAGAGGTTTa                                                                                                        >  1:543408/1‑141 (MQ=255)
                         gTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCCCGCTCTCTCGAGTTATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATa                                                                                                     >  1:567652/1‑141 (MQ=255)
                                 aCTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAg                                                                                             <  1:994439/141‑1 (MQ=255)
                                                 ttttGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAAGATTTTTa                                                                             >  1:1001964/1‑141 (MQ=255)
                                                                                ggCTTAAACCCGCAAGATCACGCTATTTCGATTGATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAgg                                              <  1:747755/141‑1 (MQ=255)
                                                                                   ttAAAACCTCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGttt                                           <  1:369215/141‑1 (MQ=255)
                                                                                   ttAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGttt                                           >  1:205969/1‑141 (MQ=255)
                                                                                      aaaCCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAAt                                        >  1:1151117/1‑141 (MQ=255)
                                                                                       aaCCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAAtt                                       >  1:484558/1‑141 (MQ=255)
                                                                                          cacAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCg                                    <  1:60715/139‑1 (MQ=255)
                                                                                          cGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCg                                    >  1:1090527/1‑141 (MQ=255)
                                                                                            cAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGa                                  <  1:310491/141‑1 (MQ=255)
                                                                                               aaTCACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGaaaa                               >  1:132999/1‑140 (MQ=255)
                                                                                                  cACGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTCCTGAGGTTTAATTTCGGAAATGGt                            <  1:767267/141‑1 (MQ=255)
                                                                                                   aCGCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGtt                           <  1:984051/141‑1 (MQ=255)
                                                                                                     gCTATTTCGATTAATAAAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGtttt                         <  1:939317/141‑1 (MQ=255)
                                                                                                               ataataaAAACAGAAAAAAGATAGACAAATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATTGTTTTAAGTGAtata               <  1:598022/140‑1 (MQ=255)
                                                                                                                 aataaAAACAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATataga             <  1:304686/141‑1 (MQ=255)
                                                                                                                      aaaCAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTCAGTGATatagaaaata        <  1:711116/141‑1 (MQ=255)
                                                                                                                       aaCAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATatagaaaatag       <  1:715048/141‑1 (MQ=255)
                                                                                                                       aaCAGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATCAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATatagaaattag       <  1:857981/141‑1 (MQ=255)
                                                                                                                       aaCAGATAAAAAATAAGCATATTTTTTTAAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATatagaaaatag       <  1:436587/141‑1 (MQ=255)
                                                                                                                          aGATAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATatagaaaatagaaa    >  1:1884/1‑141 (MQ=255)
                                                                                                                            aTAAAATATAGACATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATATAGAAAATAGAAAgt  <  1:495824/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
TATTTATAATTAAATTTAAAGGTAAGTTTTACAACTTATAAAATAAATATTTTGCTGAAAGATTTATTCAGGAAGTAAAAGGCTTAAAACCGCAAAATCACGCTATTTCGATTAATAAAAACAGATAAAATATAGAGATATTTTTTATAATGAAAGCGGTTTAATACATATGAGTTAAAAAATATTTTTAAGAATAAAGTGGGGCTTTGAATGTGCTGAGGTTTAATTTCGGAAATGGTTTTAAGTGATATAGAAAATAGAAAGT  >  CP000730/129660‑129924

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