Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,553,312 | A→T | 13.6% | L553Q (CTG→CAG) | PP_4041 ← | glycoside hydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,553,312 | 0 | A | T | 13.6% | 114.2 / 10.0 | 44 | L553Q (CTG→CAG) | PP_4041 | glycoside hydrolase |
Reads supporting (aligned to +/- strand): ref base A (14/24); new base T (4/2); total (18/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01 |
ATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCT > NC_002947/4553168‑4553454 | aTTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAgg > 2:722567/1‑150 (MQ=255) ttCGCTACCCACAACTCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGc < 1:489120/150‑1 (MQ=255) gCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGa < 2:596690/150‑1 (MQ=255) cTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAg > 1:45869/1‑150 (MQ=255) tACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaac < 2:105007/150‑3 (MQ=255) gcacagtgcCCGG‑CGCGCGACCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGtt > 1:851134/8‑150 (MQ=255) caACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTg < 1:527190/150‑1 (MQ=255) ggcctcgct‑cCCGGCACTGGTCGCGGTTCTGACGCCTCGCGGTGTTGATGTTGCCTATGTGCCTGAAGGCTTTAGGTAATTTGCCGAGCATGGCCTGGGCGGCCGATGCTTACTGTTCTGCCAGCAGGCCGAGGTCGTTGCACAGGCCgg < 2:333929/145‑1 (MQ=255) g‑cCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGgtgt < 2:45869/150‑1 (MQ=255) cGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTg < 2:693379/150‑1 (MQ=255) cgcgGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGAGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCtt > 1:367792/1‑150 (MQ=255) ggTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGcc > 1:266891/1‑106 (MQ=255) ggTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGcc < 2:266891/106‑1 (MQ=255) gTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGc > 1:294532/1‑150 (MQ=255) cAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCttt < 2:829966/150‑1 (MQ=255) cAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTg < 2:294532/150‑1 (MQ=255) gcgcGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGAACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgc < 1:412493/150‑1 (MQ=255) cgGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgcgca < 2:669301/150‑2 (MQ=255) ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCtt > 1:374731/1‑133 (MQ=255) ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCtt < 2:374731/133‑1 (MQ=255) gatgatGCCTATGTGGCTGAAGGCGTGGGGGAAGTTGACGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGTCAGTAGGCCGAGGTCGTTGCACAGGCAGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTaa < 2:139276/150‑1 (MQ=255) gatgatGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTaa < 2:485695/150‑1 (MQ=255) tgCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACAt > 2:580676/1‑150 (MQ=255) tCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGt > 2:41617/1‑150 (MQ=255) tGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca < 2:204693/150‑1 (MQ=255) cTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagcc > 2:48659/1‑150 (MQ=255) tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagcca < 1:722567/150‑1 (MQ=255) tGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCa < 1:827619/150‑1 (MQ=255) cGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGa > 1:2079/1‑150 (MQ=255) aGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAg < 1:48659/150‑1 (MQ=255) tGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGc < 1:580676/150‑1 (MQ=255) gcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagc > 2:588432/1‑118 (MQ=255) gcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagc < 1:588432/118‑1 (MQ=255) gccggGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGc < 1:250004/150‑1 (MQ=255) cGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc < 2:367792/150‑1 (MQ=255) cGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc < 1:676142/150‑1 (MQ=255) aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca < 1:669306/98‑1 (MQ=255) aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca > 2:669306/1‑98 (MQ=255) aCTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGt > 1:314190/1‑150 (MQ=255) ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca > 2:666487/2‑94 (MQ=255) ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca < 1:666487/93‑1 (MQ=255) ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCt > 2:91522/2‑150 (MQ=255) ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGCTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACCCCAGGACGGTGCCCTCTCCCGGCGTCACGccctccct > 2:344529/2‑150 (MQ=255) ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCACGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCt > 1:509129/2‑150 (MQ=255) | ATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCT > NC_002947/4553168‑4553454 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |