Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,797 G→A 48.6% intergenic (‑880/‑288) 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,7970GA48.6% 10.5 / 41.4 43intergenic (‑880/‑288)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base A (10/11);  total (21/22)
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.57e-01

TCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCT  >  CP000730/1929661‑1929932
                                                                                                                                         |                                                                                                                                           
tCATAATAAAATCGCCTACTTACTCAAACTAACAATATCTTCTTATCATGAATTTTGCCAATTAAGTGTCGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCg                                                                                                                                          <  1:231249/140‑1 (MQ=255)
       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)
                                                                                                                                         |                                                                                                                                           
TCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCT  >  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: