Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,918 A→G 26.1% intergenic (‑1001/‑167) 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,9180AG26.1% 30.6 / 12.8 23intergenic (‑1001/‑167)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/9);  new base G (2/4);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTC  >  CP000730/1929787‑1930044
                                                                                                                                     |                                                                                                                              
ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                                                       <  1:99412/141‑1 (MQ=255)
ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                                                                       <  1:86022/141‑1 (MQ=255)
ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                      <  1:45677/141‑5 (MQ=255)
  gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                    >  1:79753/12‑138 (MQ=255)
    atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc                                                                                                                   >  1:15642/1‑141 (MQ=255)
       tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTa                                                                                                                <  1:17412/141‑1 (MQ=255)
                      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                <  1:11683/137‑1 (MQ=255)
                      ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                 <  1:53373/141‑1 (MQ=255)
                         gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                                                                                              >  2:48661/1‑141 (MQ=255)
                          tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                             >  2:104924/1‑141 (MQ=255)
                            gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                           >  1:79094/1‑141 (MQ=255)
                            gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                           <  1:104924/141‑1 (MQ=255)
                            gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                           >  1:93301/1‑141 (MQ=255)
                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                  >  2:64517/1‑129 (MQ=255)
                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                 <  2:18056/141‑14 (MQ=255)
                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                 <  2:109048/141‑14 (MQ=255)
                                           cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                            <  1:4359/141‑1 (MQ=255)
                                                        aGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                               <  2:70661/141‑1 (MQ=255)
                                                                                                            aGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaaa             >  2:16685/1‑141 (MQ=255)
                                                                                                            aGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaaa             >  2:126219/1‑141 (MQ=255)
                                                                                                                 cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt        <  1:25685/141‑1 (MQ=255)
                                                                                                                 cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt        <  1:25661/141‑1 (MQ=255)
                                                                                                                       tcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTc  >  1:29710/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                              
TTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTC  >  CP000730/1929787‑1930044

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