Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,559,840 T→G 10.5% intergenic (+34/+89) bcr → / ← USA300HOU_2421 MFS family major facilitator transporter, bicyclomycin:cation symporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,559,8400TG10.5% 95.5 / 5.0 38intergenic (+34/+89)bcr/USA300HOU_2421MFS family major facilitator transporter, bicyclomycin:cation symporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/19);  new base G (2/2);  total (17/21)
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.94e-01

TTAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTGTAGT  >  CP000730/2559703‑2559979
                                                                                                                                          |                                                                                                                                           
ttAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGAtatcta                                                                                                                                          >  2:265746/1‑141 (MQ=255)
ttAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGAtatcta                                                                                                                                          >  2:127013/1‑141 (MQ=255)
            cATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAAGAAGAACTCTGACAAAGATAGATATCTTTgctagg                                                                                                                               <  2:384784/141‑6 (MQ=37)
            cATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                                                              <  1:87745/141‑1 (MQ=255)
               atatGATTATTATCTTCATTACAGCCATTCTATTAGTCATTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttc                                                                                                                           <  1:452545/141‑1 (MQ=255)
               atatGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttc                                                                                                                           <  1:203153/141‑1 (MQ=38)
                      tattatCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                                                                    <  1:176700/141‑1 (MQ=255)
                                  tACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                                                        <  2:286180/141‑1 (MQ=255)
                                  tACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                                                        <  1:236682/141‑1 (MQ=255)
                                     aGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                                                                     >  1:13893/1‑141 (MQ=255)
                                     aGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                                                                     >  1:168462/1‑141 (MQ=255)
                                               ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                                                                           <  1:223662/141‑1 (MQ=255)
                                               ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                                                                           <  1:35523/141‑1 (MQ=255)
                                                       tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                                                                                   >  1:92644/1‑141 (MQ=255)
                                                          aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                                                                >  1:160151/1‑141 (MQ=255)
                                                                 attttaaagtaaaTGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATa                                                                         >  2:318486/11‑141 (MQ=255)
                                                                  tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                                                        >  2:410138/1‑141 (MQ=255)
                                                                    tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                                                      >  1:399500/1‑141 (MQ=255)
                                                                       tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                                                                   <  2:70014/141‑1 (MQ=255)
                                                                       tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                                                                   <  1:305002/141‑1 (MQ=255)
                                                                         tAAAAAGATTAAAAACCAACATGTCACATAACACTACAACAT‑TATTAGACCCCTAGCAAAGATATCTATCTTTGTCAGGTTTCTTCTTTATGAATTATAAGATCGAATCTTCAACTAAAATTACGCCTTCCTAGCAAAGaa                                                                 <  2:271034/141‑2 (MQ=255)
                                                                                aTTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                                                          >  1:273222/1‑141 (MQ=255)
                                                                                aTTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                                                          >  1:300156/1‑141 (MQ=255)
                                                                                 ttAAAAATCAACATTTCACATAACACTTCAACAT‑AATTAGTACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATTGAAGCTGCAACTAATTTTACGCCTTCATAGCAAGGACATTTCTAt                                                         <  2:339710/141‑1 (MQ=255)
                                                                                         cAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAc                                                 >  2:23195/1‑141 (MQ=255)
                                                                                          aaCATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc                                                >  1:83329/1‑141 (MQ=255)
                                                                                                 cGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc                                         <  2:154374/141‑1 (MQ=255)
                                                                                                     tAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGc                                     <  2:312362/141‑1 (MQ=255)
                                                                                                      aaCACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCa                                    >  1:410522/1‑141 (MQ=255)
                                                                                                       acacTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCAt                                   <  1:143637/141‑1 (MQ=255)
                                                                                                          cTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCa                                >  2:117044/1‑141 (MQ=255)
                                                                                                                     aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTa                      >  1:263315/1‑141 (MQ=255)
                                                                                                                      ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTaa                     <  2:227239/141‑1 (MQ=255)
                                                                                                                      ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTaa                     <  2:227231/141‑1 (MQ=255)
                                                                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa                                          <  2:147307/119‑1 (MQ=255)
                                                                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa                                          >  1:147307/1‑119 (MQ=255)
                                                                                                                             cctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATAttt              >  2:103547/6‑141 (MQ=255)
                                                                                                                              ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATAtttt             >  2:418122/5‑141 (MQ=255)
                                                                                                                                     aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTg      <  1:301212/141‑1 (MQ=255)
                                                                                                                                     aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGAGCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTg      <  1:86030/141‑1 (MQ=255)
                                                                                                                                         atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTgtagt  <  1:67787/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TTAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACAT‑AATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTGTAGT  >  CP000730/2559703‑2559979

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