Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 281,050 | (G)5→6 | coding (934/1155 nt) | yagA ← | CP4‑6 prophage; integrase core domain‑containing protein YagA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 281,045 | 1 | . | G | 92.9% | 42.7 / ‑2.9 | 14 | coding (939/1155 nt) | yagA | CP4‑6 prophage; integrase core domain‑containing protein YagA |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (11/2); total (12/2) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTG > NC_000913/280906‑281173 | cGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTc‑c > 2:94603/1‑140 (MQ=255) tGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑‑g‑ggggg > 1:396501/1‑139 (MQ=255) tctcCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCt > 1:694382/1‑140 (MQ=255) tcCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGt > 1:183028/1‑140 (MQ=255) aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTg < 1:94603/140‑1 (MQ=255) aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTg > 2:622049/1‑140 (MQ=255) ccTTGCCTGCGCTCAGACTTCCCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCCCCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGcc < 2:87232/140‑1 (MQ=255) cgcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACGGCCGCGCGGACAGCGGATACCGCGAGCCCGGTACCGc > 1:193413/1‑140 (MQ=255) ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑‑GAGGGGCGTTGTGTTGCCGCTGTACTGCCGCCCTGACGGCCGATCCCGCGAGCCCCGCACCGCCATATCCAGCGCCTCGt > 2:291504/1‑140 (MQ=255) tttcccGCTGATATCCTCTTCCCTGCCCCCCCCCCCTTCATCGTATCCGG‑‑G‑GGGGCTTTGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGGGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa < 2:183028/136‑1 (MQ=255) tttcccGCTGATACCCACTTTCCTCCCCATCCCCCCTTCATCGTATCCGG‑‑G‑GGGGCGTTGTGTTGCCGCTGTACTCCCGCGCGGACGGTTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa < 2:396501/136‑1 (MQ=255) cTTTCCGCTGATATCCTCTTCCCTGCCCACCCCCCCTTCATCGTATTCGG‑‑G‑GGGGCGTTGTGTGGCCGCTGTCCTGCCGCGCTGACGGTTGATCCCGGGGGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa < 2:82652/140‑1 (MQ=255) atCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGAcc > 2:428430/1‑140 (MQ=255) cACTTTCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTc > 1:411258/1‑140 (MQ=255) tCCTGACCATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGACATATCCAGCGCATCGAGCGGGCGTTCAAGGTTATAGACCGTCCGCCa > 2:270067/1‑140 (MQ=255) aCCATCACCCCTTCATCGTATTC‑CGGG‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTgg > 1:158177/1‑140 (MQ=255) ccATCACCCCTTCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTc > 2:498775/1‑140 (MQ=255) tCATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCt > 1:385423/1‑140 (MQ=255) cATCGTATTC‑CG‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTg > 2:478749/1‑140 (MQ=255) | CGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑‑G‑GGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTG > NC_000913/280906‑281173 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |