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% 72.5 / NA 25intergenic (+2563/‑1229)ftsQ/USA300HOU_1125cell division protein FtsQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/12);  total (13/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAA  >  CP000730/1196375‑1196642
                                                                                                                                     |                                                                                                                                      
ttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCCGAACAAGCGTAAATACTTATAGCAAt                                                                                                                                 <  2:247398/141‑1 (MQ=255)
                           cacgATGAGATTGTTGTAACAGTTATGGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACt                                                                                                      <  1:67676/138‑1 (MQ=255)
                              gttGAGATTGTTGGAACAGTTATTTCAACTTGTTTTGATGACAAACCAACAACACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCa                                                                                                   <  2:91583/139‑1 (MQ=255)
                                            aaCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCAAATGCTTGATGGTGCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGc                                                                                     <  1:289945/141‑1 (MQ=255)
                                               aGTTATTGTAGTTTGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTCATAGCTATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGcaca                                                                                  <  2:304129/141‑1 (MQ=255)
                                               aGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATAATTCCCTTCAAATTCATCAAATGcaca                                                                                  >  2:97584/1‑141 (MQ=255)
                                               aGTTATGGCAACTGGTTTGGATGACAAGCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGcaca                                                                                  <  2:274616/141‑1 (MQ=255)
                                                    ttGCAACTGGTTTTGATAACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCaa                                                                             <  1:324973/141‑1 (MQ=255)
                                                     tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGAAACTTCTAAAGATGAATCATTCAATTCAAACTCATCAAATGCAAAAGCAAc                                                                            >  2:259897/1‑141 (MQ=255)
                                                     tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGAGTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                            >  2:219557/1‑141 (MQ=255)
                                                                           cgaaCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCAAAAGCAACTGATAGTGTAAGTGAAAGAaca                                                      >  2:201683/3‑141 (MQ=255)
                                                                              aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACAt                                                   >  2:162778/1‑141 (MQ=255)
                                                                               cATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATa                                                  <  2:148080/141‑1 (MQ=255)
                                                                                 tCACATGGTCGTAAAACTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAACAACTGGTAGTGTAAGTAAAAGAACACATaaa                                                >  2:107233/1‑139 (MQ=255)
                                                                                     aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGCCTCATTCTCGTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTa                                            >  1:290239/1‑141 (MQ=255)
                                                                                     aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTACAGCCATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATGGTGTAAGTGAAAGAACGCATACagcca                                            >  1:173683/1‑137 (MQ=255)
                                                                                      tGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAANTGATAGTGTAAGTGAAAGNACACATACAACTaa                                           <  2:201333/141‑1 (MQ=255)
                                                                                                  tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGatat                               <  1:162778/141‑1 (MQ=255)
                                                                                                         aCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCCCAAGCAACTGATAGTGTAAg                                                                <  1:242450/101‑1 (MQ=255)
                                                                                                         aCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAg                                                                >  2:242450/1‑101 (MQ=255)
                                                                                                             gATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGCTAGTGTAAGTGAAAGAACACATACAACTAAAGAAGAGGATATTCCTAGCTTCa                    >  1:351024/1‑141 (MQ=255)
                                                                                                             gATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACCCATACAACTAAAGAAGATGATATTCCTAGCTTCa                    >  1:171194/1‑141 (MQ=255)
                                                                                                                          cGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAgaga       >  2:154282/1‑141 (MQ=255)
                                                                                                                          cGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAgaga       >  1:306228/1‑141 (MQ=255)
                                                                                                                               aTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaa  <  1:338818/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAA  >  CP000730/1196375‑1196642

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