Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,307,499 G→T 16.9% K251N (AAG→AAT USA300HOU_1213 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,307,4990GT16.9% 37.2 / 4.4 24K251N (AAG→AATUSA300HOU_1213hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (6/14);  new base T (2/2);  total (8/16)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-02

ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA  >  CP000730/1307362‑1307637
                                                                                                                                          |                                                                                                                                          
aTGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATgaagaa                                                                                                                                          >  2:23506/1‑140 (MQ=255)
     gTTACCATTAAACAAAGTCGTTCAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCTAATTGAGAGAAAATGAAAATCAAGAAATACATTACCAAGAAATGAAGAAacaga                                                                                                                                     <  2:9500/140‑1 (MQ=255)
     gTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACCAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAACAAATACATTACGAAGGAATGAAGAAtcaga                                                                                                                                     <  2:176756/140‑1 (MQ=255)
         ccATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAATTGAAGAAACAGACAGt                                                                                                                                 <  1:145141/140‑1 (MQ=255)
          cATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTc                                                                                                                                >  2:85786/1‑140 (MQ=255)
                 caaAGTCGTTGAACATCCAGATTTCTTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTa                                                                                                                         <  2:67182/140‑1 (MQ=255)
                                  cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                       >  2:3555/1‑140 (MQ=255)
                                  cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                       >  2:236233/1‑140 (MQ=255)
                                     atttatttGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGga                                                                                                     >  1:203144/1‑140 (MQ=255)
                                                     aaGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAACTCAAGAAATACATTACGAAGAAATGAAGAAACAGACAATCCAAAGTAAGAATAATTTAGATTTTGGAGATCAGTTtggtggtg                                                                                     >  1:111810/1‑140 (MQ=255)
                                                            tGTACAATGGAAGATTTTGAGTATATTATCGAGTTGAGAGAAAATGAAAATCAAGAAATTCATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTGAGATTTTGGAGATCCGTTTGGTGGtggttttgg                                                                              <  1:105824/140‑1 (MQ=255)
                                                                             tGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAgtgt                                                             <  2:211742/140‑1 (MQ=255)
                                                                                    attatCGAATGCAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAaggagg                                                      <  2:163341/140‑1 (MQ=255)
                                                                                      tatCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCt                                                    <  1:34244/140‑1 (MQ=255)
                                                                                       atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                                   <  2:284409/140‑1 (MQ=255)
                                                                                       atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGGCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                                   <  1:23506/140‑1 (MQ=255)
                                                                                         cGAATTGAGAGAAAATGAAAATAAAGAAATAAAATACGATATGATTAAAGAAATGACAGAAGGAAATAATTATAACATAAATGTGGGATATCCGTTGGGTGATGGTTTTTGTTTCTTCGAGTAAGTGTAAGGAGGCTATg                                                 <  2:68560/140‑1 (MQ=255)
                                                                                                   gaaaatgaaaatCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGAAAGAATAATTTAGATTTTGGAGATCCGGTTGGTGGAGGTGTTGGTTTCTTCGAGCAAGTGTAAGGAGGCAATGAATAGAacac                                       >  2:275712/1‑140 (MQ=255)
                                                                                                               cAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTCGGTGGTGGTTTTGGTTTCATCGAGTATGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAGgc                           >  2:158912/1‑140 (MQ=255)
                                                                                                                     aTACATTACAAAGAAAAGAAGAAACAGACAGTCCAAAGTAAGAATAATGTAGATATTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAGGTGTAAGAAGGCTATGAATTGAAAACGCTCGGTGAAGCGCTaaa                     >  2:61812/1‑140 (MQ=255)
                                                                                                                        cagtACGAAGAAATGAAGAAACAGACTGTCCAAAGTAATAATGATTTAGATGTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTAGAGTAAGTGTAAGGAGGCTATGTATTGAAAACGGTCGGTGAAGCGCTAAAAgg                  <  2:177005/137‑1 (MQ=255)
                                                                                                                         aTTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGGGTAATTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGt                 <  1:184029/140‑1 (MQ=255)
                                                                                                                         aTTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa                <  1:3555/140‑1 (MQ=255)
                                                                                                                         aTTACGAAGAAATGATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa                <  1:236233/140‑1 (MQ=255)
                                                                                                                                    aaTGAAGAAACCGCCAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAGCGGTCGGTGAAGCGCTAAAAGGTAGACGTGaaa       <  1:158912/140‑1 (MQ=255)
                                                                                                                                         agaaACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTa  >  1:49128/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                          
ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGAAGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA  >  CP000730/1307362‑1307637

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