Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,196,508 C→A 100% intergenic (+2563/‑1229) ftsQ → / → USA300HOU_1125 cell division protein FtsQ/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,196,5080CA100.0% 58.7 / NA 21intergenic (+2563/‑1229)ftsQ/USA300HOU_1125cell division protein FtsQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (12/9);  total (12/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAG  >  CP000730/1196380‑1196639
                                                                                                                                |                                                                                                                                   
tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTACATACTTATAGCAATGCacc                                                                                                                         >  2:127194/1‑139 (MQ=255)
tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                         >  1:2905/1‑141 (MQ=255)
tACAGTTATTAATCCTGAATCACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTCTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACTAGCGTAAATACTTATAGCAATGCAAc                                                                                                                         >  1:163715/1‑141 (MQ=255)
tACAGTTATTAATCATGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGGAAATACTTATAGCAATGCAAc                                                                                                                         >  1:147918/1‑141 (MQ=255)
        ttaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAg                                                                                                                 >  1:162552/1‑141 (MQ=255)
           aTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATg                                                                                                              <  1:137481/141‑1 (MQ=255)
                     aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                    >  1:24253/1‑141 (MQ=255)
                     aCAAGATGAGAGTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGCAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                    <  2:112061/141‑1 (MQ=255)
                      cAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACt                                                                                                   <  2:146952/141‑1 (MQ=255)
                                    tgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTAACTTCAAATTCATCaaa                                                                                     >  2:100077/1‑141 (MQ=255)
                                        aCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAACTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATCca                                                                                 >  1:172944/1‑141 (MQ=255)
                                                     cTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATa                                                                    <  1:245411/141‑1 (MQ=255)
                                                     cTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAACCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATa                                                                    <  1:245427/141‑1 (MQ=255)
                                                                  cAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAg                                                       >  2:126498/1‑141 (MQ=255)
                                                                    aaCCAACATTACATGGTCGTATGTCTTTAAGCACTGGATTCGCAACAATCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGaa                                                     <  2:3470/141‑1 (MQ=255)
                                                                       cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAacac                                                  >  2:65033/1‑141 (MQ=255)
                                                                                                       cgagggagAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTc                  <  1:126702/134‑1 (MQ=255)
                                                                                                          ttCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCAtt               <  1:126498/141‑1 (MQ=255)
                                                                                                            cGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTAATTAg             >  1:137451/1‑141 (MQ=255)
                                                                                                             ggAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGa            <  1:47216/141‑1 (MQ=255)
                                                                                                                       tAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGagaag  >  2:228389/1‑141 (MQ=255)
                                                                                                                                |                                                                                                                                   
TACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAG  >  CP000730/1196380‑1196639

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