Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,065,589 T→C 8.2% K129E (AAG→GAG)  USA300HOU_1919 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,065,5890TC 8.2% 132.3 / 2.9 49K129E (AAG→GAG) USA300HOU_1919hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTACTTGAATGTCTAATTGTATTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTGTGTCCATTATA  >  CP000730/2065451‑2065718
                                                                                                                                          |                                                                                                                                 
gTACTTGAATGTCTAATTGTATTGCTTCAGAATAAAAGCGTTAGCCACTATAGACTATCCATCCTTCAAATGACCAGCTTTTAATATAGAGTCCGTTCATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTaa                                                                                                                                 <  1:269773/141‑1 (MQ=255)
     tGAATGTCTAATTGTATTGCTTCAGAATAAAAGCGTGCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTCATATAGAGTCCGTAAATCAAACAAGCATCAATTGCATGACCACTTCCTACATTCTTAATAGTc                                                                                                                            >  2:316678/1‑141 (MQ=255)
         tGTCTAATTGTATTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCacac                                                                                                                        <  1:265573/141‑1 (MQ=255)
            cTAATTGTATTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTa                                                                                                                     <  1:170959/141‑1 (MQ=255)
             tAATTGTATTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTAt                                                                                                                    >  2:111061/1‑141 (MQ=255)
                    aTTGCTTCAGAATAGAAGCGTTCGCCACTAGAGTCTATCAATCCTTTTGAGGGTTAGCTTTTAAGATATAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGa                                                                                                             <  1:177754/141‑1 (MQ=255)
                    aTTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGa                                                                                                             <  1:110722/141‑1 (MQ=255)
                       gCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAata                                                                                                          >  1:253830/1‑141 (MQ=255)
                          tCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCt                                                                                                       <  2:322638/141‑1 (MQ=255)
                            aGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTc                                                                                                     >  2:314537/1‑141 (MQ=255)
                               aTAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGc                                                                                                  <  1:258273/141‑1 (MQ=255)
                                tAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCGTCAATTGCATGACTACTCACTACATGCTTAATAGTCACACCTATTAATTGAATATCTtcatca                                                                                                 <  2:171961/141‑1 (MQ=255)
                                tAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCa                                                                                                 >  1:9675/1‑141 (MQ=255)
                                 aaaaGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCAt                                                                                                >  1:296265/1‑141 (MQ=255)
                                 aaaaGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCAt                                                                                                >  1:139806/1‑141 (MQ=255)
                                 aaaaGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCAt                                                                                                <  1:190118/141‑1 (MQ=255)
                                  aaaGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATg                                                                                               >  1:330227/1‑141 (MQ=255)
                                         tCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATg                                                                                        <  1:313527/141‑1 (MQ=255)
                                                 ataGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGc                                                                                <  2:200033/141‑1 (MQ=255)
                                                 ataGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGc                                                                                >  2:95022/1‑141 (MQ=255)
                                                     tCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTCAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTa                                                                            >  1:375753/1‑141 (MQ=255)
                                                     tCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTCAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTa                                                                            >  1:375752/1‑141 (MQ=255)
                                                      cTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTAt                                                                           <  1:430605/141‑1 (MQ=255)
                                                          gAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGtt                                                                       <  2:208294/141‑1 (MQ=255)
                                                              ccTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATaa                                                                   >  1:190415/1‑141 (MQ=255)
                                                               cTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAg                                                                  <  2:12022/141‑1 (MQ=255)
                                                                ttCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAgg                                                                 >  1:454544/1‑141 (MQ=255)
                                                                 tCAAATGAACAGCTTTTAATATAGAGTACGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGga                                                                <  2:61157/141‑1 (MQ=255)
                                                                        aaCAGCTTTTAATATACAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGgatacgata                                                         <  1:243121/141‑1 (MQ=255)
                                                                          cAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATaca                                                       <  1:33220/141‑1 (MQ=255)
                                                                               tttAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACA‑AGCTGTATTGTTATAAGGATACGATACACCATTc                                                 >  1:27268/1‑141 (MQ=255)
                                                                                       agagTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCaa                                          >  1:335657/1‑141 (MQ=255)
                                                                                         agTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAAc                                        >  2:153970/1‑141 (MQ=255)
                                                                                          gTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACt                                       <  1:173125/141‑1 (MQ=255)
                                                                                               ttAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTgccg                                  <  1:303049/141‑1 (MQ=255)
                                                                                                 aaTCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTgccgcc                                >  1:118906/1‑141 (MQ=255)
                                                                                                 aaTCCCACAAGCATCAATTGCATGACCACTCACTACATTCTCAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTgccgcc                                <  2:375752/141‑1 (MQ=255)
                                                                                                 aaTCCCACAAGCATCAATTGCATGACCACTCACTACATTCTCAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTgccgcc                                <  2:375753/141‑1 (MQ=255)
                                                                                                   tCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCtt                              <  1:100209/141‑1 (MQ=255)
                                                                                                    cccTCAAGCATCAATTGCATGACCACTCATTACGTTCTTAATGGTCACACCTATTAATTGAGTAGCTTCTGCATGCCCAATGCACAGAGCCGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCttt                             <  2:422236/141‑1 (MQ=255)
                                                                                                     ccACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTa                            <  1:210398/141‑1 (MQ=255)
                                                                                                           gCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAACCTTGCCGCCTTTaatatta                      >  2:314132/1‑139 (MQ=255)
                                                                                                              tCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTaatatgaata                   >  2:1987/1‑141 (MQ=255)
                                                                                                                   tGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACt              <  2:235857/141‑1 (MQ=255)
                                                                                                                       tGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTgtgt          >  2:107246/1‑141 (MQ=255)
                                                                                                                       tGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTgtgt          >  1:110067/1‑141 (MQ=255)
                                                                                                                        gACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTGTGTc         <  2:407594/141‑1 (MQ=255)
                                                                                                                          ccACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTGTGTCCa       <  2:356430/141‑1 (MQ=255)
                                                                                                                               cactACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTGTGTCCATtata  >  1:73212/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                 
GTACTTGAATGTCTAATTGTATTGCTTCAGAATAAAAGCGTTCGCCACTATAGTCTATGAATCCTTCAAATGAACAGCTTTTAATATAGAGTCCGTTAATCCCACAAGCATCAATTGCATGACCACTCACTACATTCTTAATGGTCACACCTATTAATTGAATATCTTCAGCATGCCCAATGCACATAGCTGTATTGTTATAAGGATACGATACACCATTCATATCAAACTTGCCGCCTTTAATATGAATATGACTGTGTCCATTATA  >  CP000730/2065451‑2065718

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