| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 2,096,684 | A→G | 22.1% | S246P (TCT→CCT) | USA300HOU_1958 ← | hypothetical bacteriophage protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,096,684 | 0 | A | G | 22.1% | 61.2 / 16.7 | 36 | S246P (TCT→CCT) | USA300HOU_1958 | hypothetical bacteriophage protein |
| Reads supporting (aligned to +/- strand): ref base A (16/12); new base G (4/4); total (20/16) | |||||||||||
| 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.38e-01 | |||||||||||
ATCATAAATATGTTGTGCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGAAATCACGAAAC > CP000730/2096547‑2096811 | aTCATAAATATGTTGTGCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGaaaa < 1:188387/141‑1 (MQ=255) aaaTATGTTGTGCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAg < 1:163440/141‑1 (MQ=255) aaTATGTTGTGCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGt < 1:97775/141‑1 (MQ=255) tgtgCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTta > 1:79356/1‑141 (MQ=255) tGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTACTTGCTTTATATcc > 2:81724/1‑141 (MQ=255) tcCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGcc > 2:193859/1‑141 (MQ=255) tcCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGcc > 2:93890/1‑141 (MQ=255) tcCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGcc > 2:99938/1‑141 (MQ=255) tcCGGCACCTTTACCTTTGTGTTCAACAATACATTGATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTTGTTGCTTTATATCCTTCGcc > 2:38644/1‑141 (MQ=255) acctttacctttTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCCTCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGAtt < 1:206762/141‑1 (MQ=255) tacctttTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGata < 2:69647/141‑1 (MQ=255) tacctttTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCAAATGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGata < 2:93464/141‑1 (MQ=255) acctttTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGatat > 2:46648/1‑141 (MQ=255) tttGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACaa > 2:24338/1‑141 (MQ=255) aTTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACtt > 2:152172/1‑141 (MQ=255) gTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCt > 1:200636/1‑141 (MQ=255) tGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCa > 2:238512/1‑141 (MQ=255) tGAAAATCTGTCATCGATTTAGGGAGCCCTCGTTTCGTGCCCGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGCTTGGATCTCACAAAGTCACCGCCGACTCTCCCCCCTAACTCa > 2:51915/1‑141 (MQ=255) gAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCAt < 2:112667/141‑1 (MQ=255) gAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCAt < 2:208136/141‑1 (MQ=255) tCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGaa > 1:165920/1‑141 (MQ=255) tCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATAGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGaa > 2:222431/1‑141 (MQ=255) cGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATc < 2:6448/141‑1 (MQ=255) ggAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTaa > 1:130013/1‑141 (MQ=255) tCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCtt > 1:223155/1‑141 (MQ=255) tCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCtt > 1:160261/1‑141 (MQ=255) tCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCtt > 1:160270/1‑141 (MQ=255) tCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCtt > 1:211985/1‑141 (MQ=255) cGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTa > 2:2379/1‑141 (MQ=255) tcctGCAAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCATCCCCAAGCATGGCAATCACAAAGTAACCGCCGCCCTTACCAACTCAGTCATTACTTGGAATATCTTCAGTAATCATCTTAGTcc < 1:38644/138‑1 (MQ=255) tCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACCTTACCACCTAAGTCATTACTTGGAATATCTTCAGTAATCATCTTAGTcc < 1:99938/141‑1 (MQ=255) tCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACCTTACCACCTAAGTCATTACTTGGAATATCTTCAGTAATCATCTTAGTcc < 1:93890/141‑1 (MQ=255) tCGTGCCAGCACCAACCCAACCTTTTGCATCAGGAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACCTTACCACCCAATTCATTACTTGGAATATCTTCAGTAATCATCTTAGTcc < 1:193859/141‑1 (MQ=255) aCCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGaa < 1:43961/141‑1 (MQ=255) acccaaccTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGAAATCa < 1:28636/141‑1 (MQ=255) accTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGAAATCACGaaa < 1:152172/141‑1 (MQ=255) ccTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGAAATCACGAAAc < 2:222729/141‑1 (MQ=255) | ATCATAAATATGTTGTGCTGTTCTTCCGGCACCTTTACCTTTTTGTTCAACAATACATTTATAGGTAATTTGAAAATCTGTCATCGCTTTAGGGAGCCCTCGTTTCGTGCCAGCACCAACCCAACCTTTTGCATCAGAAAAATTAGTTGCTTTATATCCTTCGCCAAGATTGGATATCACAAAGTCACCGCCGACTTTACCACCTAAATCATTACTTGGAATATCTTCAGTAATCATCTTAGTCCAACCTTTGAAATCACGAAAC > CP000730/2096547‑2096811 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |