Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,804 A→G 34.6% 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,8040AG34.6% 21.0 / 18.4 23intergenic (‑887/‑281)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/4);  new base G (3/5);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAG  >  CP000730/1929656‑1929946
                                                                                                                                                    |                                                                                                                                                 
tctTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatga                                                                                                                                                   <  1:167696/149‑1 (MQ=255)
                         ttCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                                          >  2:180660/1‑149 (MQ=255)
                                    aCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAt                                                                                                               <  1:51315/149‑1 (MQ=255)
                                     cAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAtt                                                                                                              <  2:175705/149‑1 (MQ=255)
                                            tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGc                                                                                                       >  1:127436/1‑149 (MQ=255)
                                            tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGc                                                                                                       >  1:299464/1‑149 (MQ=255)
                                             aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                      <  2:41006/149‑1 (MQ=255)
                                                 tcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                  >  2:144614/1‑149 (MQ=255)
                                                 tcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                  >  1:34623/1‑149 (MQ=255)
                                                          ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGta                                                                                        <  1:6580/149‑1 (MQ=255)
                                                                                              tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCAcc                                                  >  1:101770/1‑149 (MQ=255)
                                                                                                  aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGtc                                              >  2:36890/1‑149 (MQ=255)
                                                                                                                       aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                        <  2:256425/140‑1 (MQ=255)
                                                                                                                          tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                     <  1:16143/143‑1 (MQ=255)
                                                                                                                            gatgtgCTTGGAGATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                   <  1:61186/147‑6 (MQ=255)
                                                                                                                            gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                   <  1:47198/145‑6 (MQ=255)
                                                                                                                              tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGACTTCCTTTGCTAATCATTTTCGATCACATGATAGCCACAGTAATCAGAGGCTGTGATGTTGCGtattatac                 >  1:212493/3‑145 (MQ=255)
                                                                                                                              tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                 <  2:106189/147‑8 (MQ=255)
                                                                                                                              tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGACTGTGCTGTTGCGtattatac                 >  1:293503/3‑145 (MQ=255)
                                                                                                                                tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt               >  2:113632/1‑140 (MQ=255)
                                                                                                                                      tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc         <  2:64192/139‑1 (MQ=255)
                                                                                                                                      tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg          >  1:75537/1‑149 (MQ=255)
                                                                                                                                           tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt    >  2:244270/6‑149 (MQ=255)
                                                                                                                                            cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAAAT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa   >  2:134451/5‑149 (MQ=255)
                                                                                                                                             gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa   >  2:111276/1‑149 (MQ=255)
                                                                                                                                             gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa   >  2:173132/1‑149 (MQ=255)
                                                                                                                                             gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa   >  2:111281/1‑149 (MQ=255)
                                                                                                                                             acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGTGCATTACTTTGCTAGTCATTTTGGATCGCAGCATAGCACCATGCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa  >  1:255763/4‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAG  >  CP000730/1929656‑1929946

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