Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,825,251 T→C 9.1% intergenic (‑66/+782) cap8A_2 ← / ← icaR Capsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,825,2510TC 9.1% 122.0 / 3.4 44intergenic (‑66/+782)cap8A_2/icaRCapsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR
Reads supporting (aligned to +/- strand):  ref base T (23/17);  new base C (2/2);  total (25/19)
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.49e-01

CGGTAAAATAAGCAATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATATAATATTGTCC  >  USA300TCH1516_ALE/2825111‑2825390
                                                                                                                                            |                                                                                                                                           
cGGTAAAATAAGCAATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAAGTTTGCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCt                                                                                                                                             <  2:154326/141‑1 (MQ=255)
  gTAAAATAAGCCATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAg                                                                                                                                           <  1:253417/141‑1 (MQ=255)
         aaGCAATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACg                                                                                                                                    <  2:318690/141‑1 (MQ=255)
          aGCAATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCTGACAATGATCTTTAGTTTAAcca                                                                                                                                   >  1:261087/1‑139 (MQ=255)
                   ttAATATGCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACaaacc                                                                                                                          <  1:68585/141‑3 (MQ=255)
                      atatTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTATTTTTACGACCACAAAACATa                                                                                                                       >  2:160250/1‑141 (MQ=255)
                                          tAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTaa                                                                                                   <  1:48692/141‑1 (MQ=255)
                                               gTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAattatt                                                                                              <  2:160619/141‑1 (MQ=255)
                                                  ttACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCCAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAAttatttat                                                                                           >  2:72030/1‑141 (MQ=255)
                                                  ttACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCCAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAAttatttat                                                                                           >  2:285358/1‑141 (MQ=255)
                                                    aCTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTtatat                                                                                         >  1:333100/1‑141 (MQ=255)
                                                        aaTCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACTACAAAACATACACAACTACATTTTCTCTTATTATTTATATAAat                                                                                     >  2:248477/1‑141 (MQ=255)
                                                         aTCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAata                                                                                    >  1:340539/1‑141 (MQ=255)
                                                             aaCTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATAtttt                                                                                >  2:336809/1‑141 (MQ=255)
                                                              aCTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTa                                                                               >  2:185577/1‑141 (MQ=255)
                                                              aCTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTa                                                                               >  2:287732/1‑141 (MQ=255)
                                                              aCTTTTCCTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGCCAATTAGATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTa                                                                               >  2:417200/1‑141 (MQ=255)
                                                                 tttctttcATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCg                                                                            >  2:7395/1‑141 (MQ=255)
                                                                   tctttcATGGTTTCCTCCACATAATCAATCATTCTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATACTTTATCGtt                                                                          >  2:12771/1‑141 (MQ=255)
                                                                    ctttcATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGttt                                                                         >  2:333545/1‑141 (MQ=255)
                                                                          tGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGAAAATTATCTCTAGTTTACAGACCACATAACATACACAACTACATTTTCTTTACTTTTTTATATAAATATTTTATCGTTTAAAAtt                                                                   <  2:398546/141‑1 (MQ=255)
                                                                                  tccACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTGATTATTTATATAAATATTTTATCGTTTAAAATTATATCATg                                                           >  1:333523/1‑141 (MQ=255)
                                                                                    cacaTAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGAt                                                         >  1:139334/1‑141 (MQ=255)
                                                                                        tAATCAATCATTGTATTCATTATGTATGTTTTATAAATCTGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATtctc                                                     <  2:88911/141‑1 (MQ=255)
                                                                                        tAATCAATCATTGTATTCATTATGGATGTTTTATAAATCGGACAATTATATATAGTGTAACGACCACAAAACATACACACCTACATTTTGTCTTAGTACTTATATAAATATTTTATCGTTTAAAATTATATCATGATtctc                                                     <  2:418249/141‑1 (MQ=255)
                                                                                          atcaatcaTTGTATTCATTATGTATGTTATATAAATCTTACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTATCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCtata                                                   >  2:215466/1‑141 (MQ=255)
                                                                                            caatcaTTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTAcc                                                 >  1:63967/1‑141 (MQ=255)
                                                                                              atcaTTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCAt                                               <  1:27671/141‑1 (MQ=255)
                                                                                                      aTTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAATATTATATCATGATTCTCTACCATTATGtatc                                       >  2:376883/1‑140 (MQ=255)
                                                                                                      aTTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTCCCATTATTtata                                       >  1:322067/1‑141 (MQ=255)
                                                                                                      aTTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGtata                                       >  1:240577/1‑141 (MQ=255)
                                                                                                          aTTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACtt                                   >  2:419272/1‑141 (MQ=255)
                                                                                                          aTTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACtt                                   <  1:256811/141‑1 (MQ=255)
                                                                                                               ttatgtTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACttattta                              <  2:124817/140‑1 (MQ=255)
                                                                                                                    ttttATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATAttt                         <  1:243756/141‑1 (MQ=255)
                                                                                                                      ttATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATAttttt                       >  2:38479/1‑141 (MQ=255)
                                                                                                                          aaaTCGGACAATTTTATATAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGcac                   <  1:268844/141‑1 (MQ=255)
                                                                                                                              cGGACAATTATATCCAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGa               <  1:72030/141‑1 (MQ=255)
                                                                                                                              cGGACAATTATATCCAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGa               <  1:285358/141‑1 (MQ=255)
                                                                                                                                 aCAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGatat            <  1:56815/141‑1 (MQ=255)
                                                                                                                                    aTTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATataat         <  1:7395/141‑1 (MQ=255)
                                                                                                                                       atatCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATTTATAACTTATTTATATTTTTGCACAATATCTAATAtt      >  2:59100/1‑141 (MQ=255)
                                                                                                                                        tatCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATATAATATTg     <  2:115124/141‑1 (MQ=255)
                                                                                                                                         atCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATATAATATTGt    >  1:207374/1‑141 (MQ=255)
                                                                                                                                           cTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATATAATATTGTcc  <  1:236267/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CGGTAAAATAAGCAATAATTTAATATTCTTCTTTAGAATATTTAATAGTTTTACTAAATCAAACTTTTCTTTCATGGTTTCCTCCACATAATCAATCATTGTATTCATTATGTATGTTTTATAAATCGGACAATTATATCTAGTTTAACGACCACAAAACATACACAACTACATTTTCTCTAATTATTTATATAAATATTTTATCGTTTAAAATTATATCATGATTCTCTACCATTATGTATAACTTATTTATATTTTTGCACAAGATATAATATTGTCC  >  USA300TCH1516_ALE/2825111‑2825390

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: