Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,647,572 | Δ1 bp | coding (1150/1335 nt) | oprB‑II ← | carbohydrate‑selective porin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,647,572 | 0 | C | . | 97.2% | 152.1 / ‑2.0 | 36 | coding (1150/1335 nt) | oprB‑II | carbohydrate‑selective porin |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base . (16/19); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.72e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01 |
TCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCC > NC_002947/1647430‑1647705 | tCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATagag > 3:109618/1‑148 (MQ=255) tGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCcggc < 1:108195/148‑1 (MQ=255) acacCCCGCCTGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGcc > 4:37380/1‑148 (MQ=255) ccgcccgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGt > 4:58844/1‑147 (MQ=255) cgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCa < 1:20377/148‑1 (MQ=255) gTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCt < 7:61286/148‑1 (MQ=255) cTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTg > 2:39417/1‑147 (MQ=255) tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGc < 2:31424/147‑1 (MQ=255) tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCg > 3:131723/1‑148 (MQ=255) tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCg > 5:35404/1‑148 (MQ=255) aGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATa < 1:84712/148‑1 (MQ=255) ggTTGGGGCGTACCGTGAGCCAATCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAg < 3:133646/148‑1 (MQ=255) gCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGAcc > 4:77130/1‑120 (MQ=255) gCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGt > 5:47318/1‑148 (MQ=255) cGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCa < 3:77130/120‑1 (MQ=255) cGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGt < 8:97897/148‑1 (MQ=255) tGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGGTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGc < 5:127790/148‑1 (MQ=255) ccAAGGGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGCCCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAgg > 7:138202/1‑148 (MQ=255) tGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCg < 6:47318/148‑1 (MQ=255) gCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa > 8:120650/1‑148 (MQ=255) gCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa < 7:133770/148‑1 (MQ=255) cGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGTGTTGACGTGCACGCTGGCCAGGGCGAAACCGAt < 8:101722/147‑1 (MQ=255) cGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATg < 6:35404/148‑1 (MQ=255) gCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTg < 8:2144/148‑1 (MQ=255) tCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGCGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGgc > 4:21378/1‑148 (MQ=255) cTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCgg > 4:36096/1‑147 (MQ=255) aCTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATc < 8:103835/147‑1 (MQ=255) tGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCttt < 1:31280/148‑1 (MQ=255) tGCACTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCtt > 3:151088/1‑147 (MQ=255) aCTGGCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGaa < 5:117697/148‑1 (MQ=255) gCAGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATa > 1:141493/1‑146 (MQ=255) aGGAAGCCCGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGACACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAg > 6:116862/1‑148 (MQ=255) cccGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgc < 4:27760/148‑1 (MQ=255) ccGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc > 8:104155/1‑148 (MQ=255) ccGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc < 4:151088/148‑1 (MQ=255) ccGGGTTGTCGTAGT‑ATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc > 1:136521/1‑148 (MQ=255) | TCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCC > NC_002947/1647430‑1647705 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |