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

TTAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTGTAGT  >  CP000730/2559703‑2559979
                                                                                                                                          |                                                                                                                                           
ttAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                                                                          >  2:265746/1‑141 (MQ=255)
ttAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                                                                          >  2:127013/1‑141 (MQ=255)
            cATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAAGAAGAACTCTGACAAAGATAGATATCTTTgctagg                                                                                                                               <  2:384784/141‑6 (MQ=37)
            cATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                                                              <  1:87745/141‑1 (MQ=255)
               atatGATTATTATCTTCATTACAGCCATTCTATTAGTCATTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttc                                                                                                                           <  1:452545/141‑1 (MQ=255)
               atatGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttc                                                                                                                           <  1:203153/141‑1 (MQ=38)
                      tattatCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                                                                    <  1:176700/141‑1 (MQ=255)
                                  tACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                                                        <  2:286180/141‑1 (MQ=255)
                                  tACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                                                        <  1:236682/141‑1 (MQ=255)
                                     aGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                                                                     >  1:13893/1‑141 (MQ=255)
                                     aGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                                                                     >  1:168462/1‑141 (MQ=255)
                                               ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                                                                           <  1:223662/141‑1 (MQ=255)
                                               ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                                                                           <  1:35523/141‑1 (MQ=255)
                                                       tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                                                                                   >  1:92644/1‑141 (MQ=255)
                                                          aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                                                                >  1:160151/1‑141 (MQ=255)
                                                                 attttaaagtaaaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATa                                                                         >  2:318486/11‑141 (MQ=255)
                                                                  tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                                                        >  2:410138/1‑141 (MQ=255)
                                                                    tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                                                      >  1:399500/1‑141 (MQ=255)
                                                                       tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                                                                   <  2:70014/141‑1 (MQ=255)
                                                                       tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                                                                   <  1:305002/141‑1 (MQ=255)
                                                                         tAAAAAGATTAAAAACCAACATGTCACATAACACTACAACATTATTAGACCCCTAGCAAAGATATCTATCTTTGTCAGGTTTCTTCTTTATGAATTATAAGATCGAATCTTCAACTAAAATTACGCCTTCCTAGCAAAGaa                                                                 <  2:271034/141‑2 (MQ=255)
                                                                                aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                                                          >  1:273222/1‑141 (MQ=255)
                                                                                aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                                                          >  1:300156/1‑141 (MQ=255)
                                                                                 ttAAAAATCAACATTTCACATAACACTTCAACATAATTAGTACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATTGAAGCTGCAACTAATTTTACGCCTTCATAGCAAGGACATTTCTAt                                                         <  2:339710/141‑1 (MQ=255)
                                                                                         cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAc                                                 >  2:23195/1‑141 (MQ=255)
                                                                                          aaCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc                                                >  1:83329/1‑141 (MQ=255)
                                                                                                 cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc                                         <  2:154374/141‑1 (MQ=255)
                                                                                                     tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGc                                     <  2:312362/141‑1 (MQ=255)
                                                                                                      aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCa                                    >  1:410522/1‑141 (MQ=255)
                                                                                                       acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCAt                                   <  1:143637/141‑1 (MQ=255)
                                                                                                          cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCa                                >  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)
                                                                                                                                          |                                                                                                                                           
TTAATACATCACCATATATGATTATTATCTTCATTACAGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACAGGCATCCAAATTTCTGTAATATATTTTTCACTTGTAGT  >  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: