Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,799 C→A 48.5% intergenic (‑882/‑286) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,7990CA48.5% 14.1 / 43.0 43intergenic (‑882/‑286)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base A (9/12);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATAATAAAATCGCCTACTTTCTTAAATTAACAATATC‑TATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCT  >  CP000730/1929661‑1929932
                                                                                                                                           |                                                                                                                                         
tCATAATAAAATCGCCTACTTACTCAAACTAACAATATCTTCTTATCAT‑GAATTTTGCCAATTAAGTGTCGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCg                                                                                                                                          <  1:231249/140‑1 (MQ=255)
       aaaaTCGCCTACTTTCTTAAATTAACAATATC‑TATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGatat                                                                                                                                   <  1:129154/140‑1 (MQ=255)
                tACTTGTTTAAATTAACAAGATC‑TATTCTCATAGAATTTGTCCATTTAAGTGTAGACGATTCAATCTATCAGCTGGAAGCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTaga                                                                                                                          <  2:160247/140‑1 (MQ=37)
                          aaTTAACAATATC‑TATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGACGATGATATAATAGTAGAGTCGC‑CTAtc                                                                                                               >  1:135504/1‑140 (MQ=255)
                                        tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCaa                                                                                                  >  1:63304/1‑140 (MQ=255)
                                                   aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCag                                                                                       <  1:287278/140‑1 (MQ=255)
                                                      ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATCATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAAACGCacaga                                                                                    >  2:71650/1‑140 (MQ=255)
                                                        gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgagg                                                                                  >  1:319219/1‑140 (MQ=255)
                                                               tAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACTCAGAGAGGAGGTGta                                                                           <  1:238047/140‑1 (MQ=255)
                                                                      aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCGAGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTCTCTCTCAGGCGTCAATTTCGACGCAGAGAGGAGGGGTA‑TAA‑GG‑Tg                                                                 >  2:152419/1‑140 (MQ=255)
                                                                                   tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCg                                                    <  1:51659/140‑1 (MQ=255)
                                                                                     tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGtt                                                  >  1:178377/1‑140 (MQ=255)
                                                                                             aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCAcatcat                                          <  2:297655/140‑1 (MQ=255)
                                                                                                     aaCTTATCAGAAGCAATTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTAGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAg                                  <  1:273982/140‑1 (MQ=255)
                                                                                                            cAGAAGCAAGTGACTGTGCTTGTATATTTACCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGAACCAGTCATCAg                           >  2:86032/1‑140 (MQ=255)
                                                                                                                 ggaagtagatatgCTTGTATGTTCACAGCTGGTATAATAATAGAATCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCtg                     >  2:175165/12‑140 (MQ=255)
                                                                                                                 ggaagtagatatgCTTGTATGTTCACAGATGGTATAATAATAGATTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTGTCGTTCACAGCATAGCACCAGTCATCAGTGGcgg                     >  2:175161/12‑138 (MQ=255)
                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATCGAGTCGC‑ATATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCGGTGGCtg                     >  1:84332/12‑140 (MQ=255)
                                                                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                     <  2:251844/129‑1 (MQ=255)
                                                                                                                 ggaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                     <  2:122064/128‑1 (MQ=255)
                                                                                                                 gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTACGTTCACATCATAGCACCAGTCATCAGTGGCt                      >  2:178940/1‑140 (MQ=255)
                                                                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                    >  2:109822/11‑140 (MQ=255)
                                                                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                    <  1:250714/130‑1 (MQ=255)
                                                                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAAGAATAGAGGCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                    <  2:44671/130‑1 (MQ=255)
                                                                                                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                   <  2:221723/131‑1 (MQ=255)
                                                                                                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTCCGTTCACATCATACCACCAGTCATCAGTGGCcgtg                   >  1:28533/10‑140 (MQ=255)
                                                                                                                   aagtagatatgCTTGTATATTCACAGATGGTATAAAAATAGAG‑CGCTGAAGCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                   <  1:24688/131‑1 (MQ=255)
                                                                                                                   aagtagatatgCTTGTATATTCACAGACGGTATAATAATAGAGGCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                   <  1:295832/131‑1 (MQ=255)
                                                                                                                     gtaggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                 >  2:194116/8‑139 (MQ=255)
                                                                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                 <  1:40409/133‑2 (MQ=255)
                                                                                                                      tatctatgCTTGTATATTCACACATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                <  2:289206/134‑3 (MQ=255)
                                                                                                                      tatatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                <  2:57970/134‑3 (MQ=255)
                                                                                                                      tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGGTTATTTTCGTTCACATCATAGCACCAGTCATCAATGGCTGTGcc                 >  2:51782/1‑140 (MQ=255)
                                                                                                                         atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg             >  2:299884/4‑135 (MQ=255)
                                                                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc            >  2:297183/3‑134 (MQ=255)
                                                                                                                          aatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCAACAGTCATCAGTGGCTGTGctgtttc            >  2:133674/3‑134 (MQ=255)
                                                                                                                             gtttctataCTCGCAGATGCTAAAATCGTAGTGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑AG‑TGTTTCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt          <  1:139485/135‑1 (MQ=255)
                                                                                                                              cTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGGGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCCGTCATCATTGGCTGTGCCATTGCGtt         >  2:106957/1‑140 (MQ=255)
                                                                                                                              cTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGAGGATGTGCCATTGCGtt         >  2:264549/1‑140 (MQ=255)
                                                                                                                               ttGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt       <  2:164329/140‑4 (MQ=255)
                                                                                                                                 gTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt      >  1:251582/1‑140 (MQ=255)
                                                                                                                                     atTTGCCGATGATATAATAGTATAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑CGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt  <  1:81347/140‑1 (MQ=255)
                                                                                                                                     atTTGCCGAGGATACAATAGTAGAGTCGC‑CTATTTCACAGGCGTCAACTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑GGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt  <  1:81341/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TCATAATAAAATCGCCTACTTTCTTAAATTAACAATATC‑TATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCT  >  CP000730/1929661‑1929932

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