Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,804 A→G 50.0% intergenic (‑887/‑281) 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,8040AG50.0% ‑1.4 / 50.3 41intergenic (‑887/‑281)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/8);  new base G (9/11);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  CP000730/1929668‑1929941
                                                                                                                                        |                                                                                                                                             
aaaaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGatat                                                                                                                                            <  1:129154/140‑1 (MQ=255)
         tACTTGTTTAAATTAACAAGATCTATTCTCATAGAATTTGTCCATTTAAGTGTAGACGATTCAATCTATCAGCTGGAAGCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTaga                                                                                                                                   <  2:160247/140‑1 (MQ=37)
                   aaTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGACGATGATATAATAGTAGAGTCGC‑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)
                                                                                                                               tcacagatggtATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt        >  1:128815/6‑140 (MQ=255)
                                                                                                                                 acagatggtATAATAATAGAGACGC‑CTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTCTGCTGTTGCGTATTCTAATTATTg      >  2:266247/4‑140 (MQ=255)
                                                                                                                                    gatggtATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtttactaagccaataag   <  2:196909/140‑15 (MQ=255)
                                                                                                                                    gatggtATAATAATAGAGTCGC‑CTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt   <  1:170002/140‑1 (MQ=255)
                                                                                                                                    gatgatATAATAGTAGAGGCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc    <  1:38697/140‑1 (MQ=255)
                                                                                                                                     atggtATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATCCTTAATGGCtt  >  2:11187/1‑139 (MQ=255)
                                                                                                                                      tgatATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa  >  1:292239/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                             
AAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  CP000730/1929668‑1929941

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: ‑