| 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 | ||
| * | CP000730 | 1,307,499 | 0 | G | T | 16.9% | 37.2 / 4.4 | 24 | K251N (AAG→AAT) | USA300HOU_1213 | hypothetical 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 | |||||||||||
ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA  >  CP000730/1307362‑1307637                                                                                                                                          |                                                                                                                                          aTGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATgaagaa                                                                                                                                          >  2:23506/1‑140 (MQ=255)     gTTACCATTAAACAAAGTCGTTCAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCTAATTGAGAGAAAATGAAAATCAAGAAATACATTACCA‑AGAAATGAAGAAacaga                                                                                                                                     <  2:9500/140‑1 (MQ=255)     gTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACCAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAACAAATACATTACGA‑AGGAATGAAGAAtcaga                                                                                                                                     <  2:176756/140‑1 (MQ=255)         ccATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAATTGAAGAAACAGACAGt                                                                                                                                 <  1:145141/140‑1 (MQ=255)          cATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTc                                                                                                                                >  2:85786/1‑140 (MQ=255)                 caaAGTCGTTGAACATCCAGATTTCTTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTa                                                                                                                         <  2:67182/140‑1 (MQ=255)                                  cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                       >  2:3555/1‑140 (MQ=255)                                  cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg                                                                                                       >  2:236233/1‑140 (MQ=255)                                     atttatttGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGga                                                                                                     >  1:203144/1‑140 (MQ=255)                                                     aaGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAACTCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAATCCAAAGTAAGAATAATTTAGATTTTGGAGATCAGTTtggtggtg                                                                                     >  1:111810/1‑140 (MQ=255)                                                            tGTACAATGGAAGATTTTGAGTATATTATCGAGTTGAGAGAAAATGAAAATCAAGAAATTCATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTGAGATTTTGGAGATCCGTTTGGTGGtggttttgg                                                                              <  1:105824/140‑1 (MQ=255)                                                                             tGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAgtgt                                                             <  2:211742/140‑1 (MQ=255)                                                                                    attatCGAATGCAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAaggagg                                                      <  2:163341/140‑1 (MQ=255)                                                                                      tatCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCt                                                    <  1:34244/140‑1 (MQ=255)                                                                                       atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                                   <  2:284409/140‑1 (MQ=255)                                                                                       atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGGCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa                                                   <  1:23506/140‑1 (MQ=255)                                                                                         cGAATTGAGAGAAAATGAAAATAAAGAAATAAAATACGATATGATTAAAGAAA‑TGACAGAAGGAAATAATTATAACATAAATGTGGGATATCCGTTGGGTGATGGTTTTTGTTTCTTCGAGTAAGTGTAAGGAGGCTATg                                                 <  2:68560/140‑1 (MQ=255)                                                                                                   gaaaatgaaaatCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGAAAGAATAATTTAGATTTTGGAGATCCGGTTGGTGGAGGTGTTGGTTTCTTCGAGCAAGTGTAAGGAGGCAATGAATAGAacac                                       >  2:275712/1‑140 (MQ=255)                                                                                                               cAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTCGGTGGTGGTTTTGGTTTCATCGAGTATGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAGgc                           >  2:158912/1‑140 (MQ=255)                                                                                                                     aTACATTACAA‑AGAAAAGAAGAAACAGACAGTCCAAAGTAAGAATAATGTAGATATTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAGGTGTAAGAAGGCTATGAATTGAAAACGCTCGGTGAAGCGCTaaa                     >  2:61812/1‑140 (MQ=255)                                                                                                                        cagtACGA‑AGAAATGAAGAAACAGACTGTCCAAAGTAATAATGATTTAGATGTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTAGAGTAAGTGTAAGGAGGCTATGTATTGAAAACGGTCGGTGAAGCGCTAAAAgg                  <  2:177005/137‑1 (MQ=255)                                                                                                                         aTTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGGGTAATTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGt                 <  1:184029/140‑1 (MQ=255)                                                                                                                         aTTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa                <  1:3555/140‑1 (MQ=255)                                                                                                                         aTTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa                <  1:236233/140‑1 (MQ=255)                                                                                                                                    aaTGAAGAAACCGCCAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAGCGGTCGGTGAAGCGCTAAAAGGTAGACGTGaaa       <  1:158912/140‑1 (MQ=255)                                                                                                                                         agaaACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTa  >  1:49128/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                          ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA  >  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: ‑ |