Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 495,289 G→A 7.7% Q277* (CAG→TAG)  USA300HOU_0470 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730495,2890GA 7.7% 154.5 / 2.6 52Q277* (CAG→TAG) USA300HOU_0470hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (25/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 = 8.55e-01

TTTAAAAAGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTT  >  CP000730/495150‑495426
                                                                                                                                           |                                                                                                                                         
tttAAAAAGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTGCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGt                                                                                                                                          <  2:132466/141‑1 (MQ=255)
 tgaaaaaGATTAAAAACATAGGTATGGCACCGGATAAATATGGGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTa                                                                                                                                         <  2:671625/139‑1 (MQ=255)
  tAAAAAGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTaa                                                                                                                                        >  1:276932/1‑141 (MQ=255)
     aaaGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTATAAAAt                                                                                                                                     >  2:400105/1‑141 (MQ=255)
     aaaGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTATAAAAt                                                                                                                                     >  2:400128/1‑141 (MQ=255)
        gATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTc                                                                                                                                  <  1:555457/141‑1 (MQ=255)
         aTTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCa                                                                                                                                 <  2:697136/141‑1 (MQ=255)
         aTTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCa                                                                                                                                 <  1:326809/141‑1 (MQ=255)
                aCATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATAc                                                                                                                          <  2:62547/141‑1 (MQ=255)
                  aTAGGTGGGGCACCGGATAAATATGAGAATAGTAAGATGTGAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATCTCATGCATACCa                                                                                                                        <  2:275754/141‑1 (MQ=255)
                  aTAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCa                                                                                                                        >  2:672086/1‑141 (MQ=255)
                               cGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTTATGCTTACCACTAGCAATTGTat                                                                                                           >  1:706879/1‑140 (MQ=255)
                                  aaaaaTATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTg                                                                                                        <  1:339383/139‑1 (MQ=255)
                                    aaaTATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCa                                                                                                      <  1:610363/141‑1 (MQ=255)
                                      atatGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAAc                                                                                                    >  1:635457/1‑141 (MQ=255)
                                         tGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATc                                                                                                 <  1:4783/141‑1 (MQ=255)
                                         tGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATc                                                                                                 >  1:434324/1‑141 (MQ=255)
                                                 ataaGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCataataa                                                                                         >  2:490687/1‑140 (MQ=255)
                                                 ataaGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAg                                                                                         <  2:547871/141‑1 (MQ=255)
                                                   aaGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTATAAAATTTCCTGCATACCACTAGCAATTGTAATTGCAACATCCATAATAgcg                                                                                       <  1:400128/141‑1 (MQ=255)
                                                   aaGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTATAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAgcg                                                                                       <  1:400105/141‑1 (MQ=255)
                                                     gATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATCCCACTAGCAATTGTAATTGCAGCATCCATCATAGCGcc                                                                                     >  1:241824/1‑141 (MQ=255)
                                                     gATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGcc                                                                                     >  2:551614/1‑141 (MQ=255)
                                                           aGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGa                                                                               >  2:29933/1‑141 (MQ=255)
                                                           aGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCAAAGTGa                                                                               >  1:62455/1‑141 (MQ=255)
                                                                 tGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTAt                                                                         >  2:114218/1‑141 (MQ=255)
                                                                    tCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACTCTTAGCGCGTGTAATTGCAACGTCCATAAGAGCGCCAAGTGAACCtcgtaa                                                                      >  1:670470/1‑136 (MQ=255)
                                                                     cccATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTaaa                                                                     >  2:655845/1‑140 (MQ=255)
                                                                             atCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTATCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTTAGGc                                                             <  1:518070/141‑1 (MQ=255)
                                                                                      gATGTTTCGCCCAGCAAGCGCCCCTCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGat                                                    <  1:632089/141‑1 (MQ=255)
                                                                                         gtctagCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATc                                                 <  1:588114/136‑1 (MQ=255)
                                                                                          tttCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCGTGAATACCACTAGCAATTGTAATTGCAACGTCCATAATAGCGCCAAGTGAACCTATAAACCGTGAGTCTAGGAAGATATCt                                                >  1:634682/1‑141 (MQ=255)
                                                                                                     aaGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTaa                                     <  1:346749/141‑1 (MQ=255)
                                                                                                         gcCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGAt                                 >  2:22313/1‑141 (MQ=255)
                                                                                                         gcCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGCCTAGGAAGATATCTTTCGGTGGTAATGAt                                 >  1:598569/1‑141 (MQ=255)
                                                                                                          cccATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATAGCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATa                                >  2:558275/1‑141 (MQ=255)
                                                                                                            aatatcctcATTGAAATGTGTGGCCTACGCTGGAAAATTTCATGCATACCAGTAGCGAGTGTAATTGCAACATCCATAATAGCGACAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTCTCGGTGGTAATGATaaa                              <  2:195572/137‑1 (MQ=255)
                                                                                                                 cctcATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGtt                         <  1:138414/141‑1 (MQ=255)
                                                                                                                  ctcATTGAAATGTGTGGCGTATGCTGTAAAATGTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCAATGAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTc                        <  1:309277/141‑1 (MQ=255)
                                                                                                                  ctcATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTTATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTc                        <  1:651502/141‑1 (MQ=255)
                                                                                                                  ctcATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTc                        <  1:309304/141‑1 (MQ=255)
                                                                                                                      ttGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCg                    <  1:490687/141‑1 (MQ=255)
                                                                                                                      ttGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCg                    <  1:114218/141‑1 (MQ=255)
                                                                                                                      ttGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCg                    <  1:627411/141‑1 (MQ=255)
                                                                                                                         aaaTGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGttt                 >  1:659557/1‑141 (MQ=255)
                                                                                                                          aaTGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTc                >  2:411068/1‑141 (MQ=255)
                                                                                                                          aaTGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTc                >  1:228646/1‑141 (MQ=255)
                                                                                                                               gtgGCGTACGCTGTAAAATGTCATGCATACCACTAGAGAATGTGATTGCACATTCAATAATAGCGCCAAGTGAACCTAGTATCACTGAGGCGCGGAAGATAGATATCGGTGGTAAGGGTAAAAGGTTCATCGGTTcgcaat           >  2:124880/1‑136 (MQ=255)
                                                                                                                                tgGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACCGAGGCTAGGAAGATAACTTTCGGTGGTAATGAAAAAAAGTTCATCGTTTCATAttt          >  2:31190/1‑141 (MQ=255)
                                                                                                                                 ggCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTa         >  1:617135/1‑141 (MQ=255)
                                                                                                                                      aCGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGcc    <  1:452746/141‑1 (MQ=255)
                                                                                                                                        gCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCtt  <  2:635457/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TTTAAAAAGATTAAAAACATAGGTAGTGCACCGGATAAATATGAGAATAATAAGATGTTAGTCATTGTTCCCATAATATCTTGGCCGATGTTTCGCCCAGCAAGCGCCCATCTCCTCATTGAAATGTGTGGCGTACGCTGTAAAATTTCATGCATACCACTAGCAATTGTAATTGCAACATCCATAATAGCGCCAAGTGAACCTATTAACACTGAGGCTAGGAAGATATCTTTCGGTGGTAATGATAAAAAGTTCATCGTTTCATATTTAATGCCTT  >  CP000730/495150‑495426

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: