Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 133,296 | C→G | 14.3% | P25R (CCG→CGG) | PP_0128 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 133,296 | 0 | C | G | 14.3% | 70.3 / 5.8 | 28 | P25R (CCG→CGG) | PP_0128 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (17/7); new base G (2/2); total (19/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 |
GTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCA > NC_002947/133160‑133445 | gTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCCTTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAg > 1:829297/1‑150 (MQ=255) gACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcag > 1:229614/1‑130 (MQ=255) gACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcag < 2:229614/132‑3 (MQ=255) cAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCAt < 2:357932/150‑1 (MQ=255) ccGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcag < 1:6722/79‑3 (MQ=255) ccGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCGGGGCCTGcag > 2:6722/1‑77 (MQ=255) ccGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa < 1:380543/104‑1 (MQ=255) ccGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa > 2:380543/1‑104 (MQ=255) cAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCATGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGAc < 2:313435/150‑1 (MQ=255) agagCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATg > 2:548205/1‑79 (MQ=255) agagCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATg < 1:548205/79‑1 (MQ=255) aCGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGt < 2:327093/147‑1 (MQ=255) aCGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGt > 1:327093/1‑147 (MQ=255) tGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCATCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGgca < 1:746794/150‑1 (MQ=255) tGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACct > 2:14778/1‑150 (MQ=255) ccTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgct > 1:760261/1‑150 (MQ=255) cTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgctg > 1:139128/1‑150 (MQ=255) cTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgctg > 2:62493/1‑150 (MQ=255) cGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACg > 2:514460/1‑150 (MQ=255) cGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACg > 1:473632/1‑150 (MQ=255) gtaGGTCAACGAGCATCACCTGCAGGCACCGGGCCTGCCCGAGGATGGACGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTAACTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg > 2:509921/3‑150 (MQ=255) aaCGAGCATCACCTGAAGGGCGCGGGCCTGCCCGAGGATGGCGGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGcc > 2:847076/1‑150 (MQ=255) cGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCgg > 1:671575/1‑150 (MQ=255) tCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCt > 1:631108/1‑150 (MQ=255) gCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCagaa < 1:62493/150‑1 (MQ=255) gCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCagaa > 1:14647/1‑150 (MQ=255) aGGCCCCGGGCCTGCCCGAGGATGGACCTACTCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGccggtct > 1:450767/1‑144 (MQ=255) gCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATc > 2:57346/1‑150 (MQ=255) cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAgc > 1:57913/1‑71 (MQ=255) cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAgc < 2:57913/71‑1 (MQ=255) cccGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCgg > 1:462264/1‑150 (MQ=255) cGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCa < 1:14778/150‑1 (MQ=255) | GTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCA > NC_002947/133160‑133445 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |