Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,247,298 | A→G | D238D (GAT→GAC) | truB ← | tRNA pseudouridine synthase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,247,298 | 0 | A | G | 91.4% | 204.4 / 8.3 | 58 | D238D (GAT→GAC) | truB | tRNA pseudouridine synthase B |
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (31/22); total (34/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01 |
CAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGT > REL606/3247228‑3247369 | cAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCaa < 2:254554/75‑1 (MQ=255) gaagtCGTTAACGGAAGATTCACCCCCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGt < 1:1190531/75‑1 (MQ=255) aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTa < 2:695308/75‑1 (MQ=255) aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTa < 2:631535/75‑1 (MQ=255) aagaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTa < 2:230969/75‑1 (MQ=255) gaCGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAAc < 2:582591/75‑1 (MQ=255) cGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTc < 1:634568/75‑1 (MQ=255) tAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCcgc > 1:1079553/1‑75 (MQ=255) tAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCcgc < 2:590451/75‑1 (MQ=255) gATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGaa < 2:1066653/75‑1 (MQ=255) gATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGaa > 1:388237/1‑75 (MQ=255) aTTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAat > 2:923088/1‑75 (MQ=255) aTTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAat < 2:545315/75‑1 (MQ=255) tCACCGCCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAatat < 2:115187/75‑1 (MQ=255) caccacCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATc > 2:478479/1‑75 (MQ=255) accaccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATcc > 1:1350044/1‑75 (MQ=255) ccaccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCct < 2:516673/75‑1 (MQ=255) caccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCctg > 2:1285882/1‑75 (MQ=255) caccGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAACATcccg > 2:396746/1‑73 (MQ=255) accGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCctgc < 1:3261/75‑1 (MQ=255) ggATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATTGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGt > 1:679031/1‑75 (MQ=255) ggATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGt > 1:1469817/1‑75 (MQ=255) gATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGtt < 1:1418211/75‑1 (MQ=255) tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa < 2:1378971/75‑1 (MQ=255) tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa < 2:373530/75‑1 (MQ=255) tAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCa > 1:739969/1‑75 (MQ=255) aGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAg > 2:221993/1‑75 (MQ=255) tCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCt > 2:1466258/1‑75 (MQ=255) aGCTGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTc > 2:1244223/1‑75 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac > 1:772107/1‑75 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac > 1:1112390/1‑75 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac > 2:1578446/1‑75 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac < 1:408652/75‑1 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac < 1:408637/75‑1 (MQ=255) tGGACTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTcaac > 1:206029/1‑75 (MQ=255) cTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGt < 2:1585565/75‑1 (MQ=255) cTGTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGt > 1:766074/1‑75 (MQ=255) gTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTc < 2:1586580/75‑1 (MQ=255) gTCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTc < 1:32538/75‑1 (MQ=255) tCCATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCa > 2:1391892/1‑75 (MQ=255) ccATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCAc < 1:1511923/75‑1 (MQ=255) ccATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCAc > 2:416591/1‑75 (MQ=255) cATTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACg > 1:1590772/1‑75 (MQ=255) aTTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGc < 2:109749/75‑1 (MQ=255) aTTGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGc > 1:926592/1‑75 (MQ=255) ttGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCa > 2:665289/1‑75 (MQ=255) ttGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCa > 2:434581/1‑75 (MQ=255) tGGCATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAg > 2:1123195/1‑75 (MQ=255) cATCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTg > 1:1052921/1‑75 (MQ=255) tCAGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCt > 2:1060669/1‑75 (MQ=255) aGTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTcc > 1:1133810/1‑75 (MQ=255) gTAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCa > 1:1523320/1‑75 (MQ=255) tAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAg > 1:878604/1‑75 (MQ=255) tAATGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAg > 1:1198574/1‑75 (MQ=255) aaTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAgg > 1:775640/1‑75 (MQ=255) aTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAggg > 1:1246130/1‑75 (MQ=255) aTGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAggg < 1:1027142/75‑1 (MQ=255) tGGGTCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGt > 1:1474935/1‑75 (MQ=255) | CAGAAGACGTTAACGGAAGATTCACCACCGGATAGTCCGAAGCTGGACTGTCCATTGGCATCAGTAATGGATCAAGTAACTCCGCGGCTGGAATATCCTGCTGTTCAGCTTGCTCAACAAGTTCACGCAGGTGCTCCAGGGT > REL606/3247228‑3247369 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |