Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,510,547 G→A 13.8% intergenic (+669/‑404) USA300HOU_2370 → / → USA300HOU_2372 hypothetical protein/MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,510,5470GA13.8% 78.4 / 4.8 29intergenic (+669/‑404)USA300HOU_2370/USA300HOU_2372hypothetical protein/MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base A (2/2);  total (15/14)
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

CTTGGAGATCGCCAAAAACCACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTGGACCGTGGGTAACA  >  CP000730/2510412‑2510684
                                                                                                                                       |                                                                                                                                         
cTTGGAGATCGCCAAAAACCACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTa                                                                                                                                      >  1:591/1‑141 (MQ=255)
             aaaaaCCACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATgc                                                                                                                         <  2:155982/141‑1 (MQ=255)
             aaaaaCCACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATgc                                                                                                                         <  2:100877/141‑1 (MQ=255)
                  ccACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTAt                                                                                                                    <  2:61494/141‑1 (MQ=255)
                    aCTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCa                                                                                                                  >  2:79488/1‑141 (MQ=255)
                     cTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCaa                                                                                                                 >  2:145640/1‑141 (MQ=255)
                               taaataaaAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATattgatt                                                                                                       >  2:27603/1‑141 (MQ=255)
                                                 gtgATGGTAGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTttga                                                                                     <  1:12786/141‑1 (MQ=255)
                                                  tgATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATAAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTttgat                                                                                    >  2:71257/1‑141 (MQ=255)
                                                  tgATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATAAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTttgat                                                                                    >  2:226371/1‑141 (MQ=255)
                                                            aTACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGatt                                                                          >  2:249318/1‑141 (MQ=255)
                                                                             gATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGc                                                         >  1:93308/1‑141 (MQ=255)
                                                                              atttatttAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCt                                                        <  1:96186/141‑1 (MQ=255)
                                                                                      aGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGACGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTa                                                <  1:193836/141‑1 (MQ=255)
                                                                                      aGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTa                                                <  1:258520/141‑1 (MQ=255)
                                                                                        tGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATg                                              >  2:206368/1‑141 (MQ=255)
                                                                                        tGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATg                                              >  2:263380/1‑141 (MQ=255)
                                                                                         gCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGt                                             <  2:102171/141‑1 (MQ=255)
                                                                                         gCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGt                                             <  2:34718/141‑1 (MQ=255)
                                                                                         gCCGAGTACAAGTACAATTAATTATATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGt                                             <  1:120381/141‑1 (MQ=255)
                                                                                                  aagtacaaTTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCAtt                                    >  2:195303/1‑141 (MQ=255)
                                                                                                       caaTTAATTCTATGAAAAGCTTAATTGCACATAAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATg                               <  1:71257/141‑1 (MQ=255)
                                                                                                       caaTTAATTCTATGAAAAGCTTAATTGCACATAAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATg                               <  1:226371/141‑1 (MQ=255)
                                                                                                        aattaattCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATgg                              <  2:591/141‑1 (MQ=255)
                                                                                                                      aaaGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTg                >  1:218894/1‑141 (MQ=255)
                                                                                                                          cTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATCTCCTGACTGGGATGGCGGAAGTCTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAGAGTGAATAGTGGAcc            >  2:24036/1‑141 (MQ=255)
                                                                                                                             aaTTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTGGACCGTg         <  2:221684/141‑1 (MQ=255)
                                                                                                                                    caTGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTGGACCGTGGGTAACa  >  1:86591/1‑141 (MQ=255)
                                                                                                                                    caTGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTGGACCGTGGGTAACa  >  1:31096/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CTTGGAGATCGCCAAAAACCACTTATGTATCTAAATAAAAATCGCGGTGGTGATGGTGGCATACCAGGTTATATTATGATTTATTTAGTGCCGAGTACAAGTACAATTAATTCTATGAAAAGCTTAATTGCACATGAAGTAAATCATAATATGCGCTATCAATATATTGATTGGGATGGCGGAAGTTTGATTGAATTGATTATTGCAGAAGGATTAGCTGAAAACTATGTAGAATCATTGTATGGTAAAGCGCATATTGGACCGTGGGTAACA  >  CP000730/2510412‑2510684

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: