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

AAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  CP000730/1929668‑1929941
                                                                                                                                        |                                                                                                                                             
aaaaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGatat                                                                                                                                            <  1:129154/140‑1 (MQ=255)
         tACTTGTTTAAATTAACAAGATCTATTCTCATAGAATTTGTCCATTTAAGTGTAGACGATTCAATCTATCAGCTGGAAGCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTaga                                                                                                                                   <  2:160247/140‑1 (MQ=37)
                   aaTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGACGATGATATAATAGTAGAGTCGCCTAtc                                                                                                                        >  1:135504/1‑140 (MQ=255)
                                tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCaa                                                                                                           >  1:63304/1‑140 (MQ=255)
                                           aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                <  1:287278/140‑1 (MQ=255)
                                              ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATCATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAAACGCacaga                                                                                             >  2:71650/1‑140 (MQ=255)
                                                gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagg                                                                                           >  1:319219/1‑140 (MQ=255)
                                                       tAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACTCAGAGAGGAGGTGta                                                                                    <  1:238047/140‑1 (MQ=255)
                                                              aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCGAGTATATTTGCCGATGATATAATAGTAGAGTCGCCTCTCTCTCAGGCGTCAATTTCGACGCAGAGAGGAGGGGTATAAGGTg                                                                          >  2:152419/1‑140 (MQ=255)
                                                                           tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCg                                                             <  1:51659/140‑1 (MQ=255)
                                                                             tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGtt                                                           >  1:178377/1‑140 (MQ=255)
                                                                                     aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCAcatcat                                                   <  2:297655/140‑1 (MQ=255)
                                                                                             aaCTTATCAGAAGCAATTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTAGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAg                                           <  1:273982/140‑1 (MQ=255)
                                                                                                    cAGAAGCAAGTGACTGTGCTTGTATATTTACCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGAACCAGTCATCAg                                    >  2:86032/1‑140 (MQ=255)
                                                                                                         ggaagtagatatgCTTGTATGTTCACAGCTGGTATAATAATAGAATCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCtg                              >  2:175165/12‑140 (MQ=255)
                                                                                                         ggaagtagatatgCTTGTATGTTCACAGATGGTATAATAATAGATTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTGTCGTTCACAGCATAGCACCAGTCATCAGTGGcgg                              >  2:175161/12‑138 (MQ=255)
                                                                                                         ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATCGAGTCGCATATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCGGTGGCtg                              >  1:84332/12‑140 (MQ=255)
                                                                                                         ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                              <  2:251844/129‑1 (MQ=255)
                                                                                                         ggaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                              <  2:122064/128‑1 (MQ=255)
                                                                                                         gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTACGTTCACATCATAGCACCAGTCATCAGTGGCt                               >  2:178940/1‑140 (MQ=255)
                                                                                                          gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                             >  2:109822/11‑140 (MQ=255)
                                                                                                          gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                             <  1:250714/130‑1 (MQ=255)
                                                                                                          gaagtagatatgCTTGTATATTCACAGATGGTATAAGAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                             <  2:44671/130‑1 (MQ=255)
                                                                                                           aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                            <  2:221723/131‑1 (MQ=255)
                                                                                                           aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTCCGTTCACATCATACCACCAGTCATCAGTGGCcgtg                            >  1:28533/10‑140 (MQ=255)
                                                                                                           aagtagatatgCTTGTATATTCACAGATGGTATAAAAATAGAGCGCTGAAGCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                            <  1:24688/131‑1 (MQ=255)
                                                                                                           aagtagatatgCTTGTATATTCACAGACGGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                            <  1:295832/131‑1 (MQ=255)
                                                                                                             gtaggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                          >  2:194116/8‑139 (MQ=255)
                                                                                                             gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                          <  1:40409/133‑2 (MQ=255)
                                                                                                              tatctatgCTTGTATATTCACACATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                         <  2:289206/134‑3 (MQ=255)
                                                                                                              tatatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                         <  2:57970/134‑3 (MQ=255)
                                                                                                              tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGGTTATTTTCGTTCACATCATAGCACCAGTCATCAATGGCTGTGcc                          >  2:51782/1‑140 (MQ=255)
                                                                                                                 atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                      >  2:299884/4‑135 (MQ=255)
                                                                                                                  tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                     >  2:297183/3‑134 (MQ=255)
                                                                                                                  aatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCAACAGTCATCAGTGGCTGTGctgtttc                     >  2:133674/3‑134 (MQ=255)
                                                                                                                     gtttctataCTCGCAGATGCTAAAATCGTAGTGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAAGTGTTTCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt                   <  1:139485/135‑1 (MQ=255)
                                                                                                                      cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGGGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCCGTCATCATTGGCTGTGCCATTGCGtt                  >  2:106957/1‑140 (MQ=255)
                                                                                                                      cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGAGGATGTGCCATTGCGtt                  >  2:264549/1‑140 (MQ=255)
                                                                                                                       ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                <  2:164329/140‑4 (MQ=255)
                                                                                                                         gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt               >  1:251582/1‑140 (MQ=255)
                                                                                                                             atTTGCCGATGATATAATAGTATAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGCGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt           <  1:81347/140‑1 (MQ=255)
                                                                                                                             atTTGCCGAGGATACAATAGTAGAGTCGCCTATTTCACAGGCGTCAACTTAGACGCAGAGAGGAGGTGTATAAGGGGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt           <  1:81341/140‑1 (MQ=255)
                                                                                                                               tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt        >  1:128815/6‑140 (MQ=255)
                                                                                                                                 acagatggtATAATAATAGAGACGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTCTGCTGTTGCGTATTCTAATTATTg      >  2:266247/4‑140 (MQ=255)
                                                                                                                                    gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtttactaagccaataag   <  2:196909/140‑15 (MQ=255)
                                                                                                                                    gatggtATAATAATAGAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt   <  1:170002/140‑1 (MQ=255)
                                                                                                                                    gatgatATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc    <  1:38697/140‑1 (MQ=255)
                                                                                                                                     atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATCCTTAATGGCtt  >  2:11187/1‑139 (MQ=255)
                                                                                                                                      tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa  >  1:292239/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                             
AAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTA  >  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: