Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,259,207 | G→A | A33V (GCG→GTG) | glmM ← | phosphoglucosamine mutase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,259,207 | 0 | G | A | 96.3% | 198.0 / ‑2.8 | 54 | A33V (GCG→GTG) | glmM | phosphoglucosamine mutase |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (27/25); total (28/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
ATATAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATC > REL606/3259135‑3259277 | atatAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCcacg > 1:392352/1‑75 (MQ=255) atatAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCcacg > 2:628039/1‑75 (MQ=255) tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg < 2:1128952/75‑1 (MQ=255) tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg > 2:1171209/1‑75 (MQ=255) tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg > 2:1462269/1‑75 (MQ=255) tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg < 2:820490/75‑1 (MQ=255) tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg > 1:742042/1‑75 (MQ=255) aGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGccc > 1:1159877/1‑75 (MQ=255) aGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGccc > 1:439850/1‑75 (MQ=255) gCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCa > 1:1357850/1‑75 (MQ=255) ccAGAAATCCGCGTGTCTTTCCAAATAATGTTCTTCCGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCaa < 1:893568/75‑1 (MQ=255) aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc > 2:14982/1‑75 (MQ=255) aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc > 1:1242394/1‑75 (MQ=255) aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc > 1:1220719/1‑75 (MQ=255) aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc > 1:1164048/1‑75 (MQ=255) aTACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGc > 1:332161/1‑75 (MQ=255) tACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCt > 2:338682/1‑75 (MQ=255) aCGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCCGCtt > 2:1155131/1‑75 (MQ=255) aCGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCACCtt > 2:1507915/1‑75 (MQ=255) gcgTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTaa < 1:881521/75‑1 (MQ=255) cgTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAg < 2:462579/75‑1 (MQ=255) gtgtCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGc < 2:852087/75‑1 (MQ=255) tgtCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGca > 1:587720/1‑75 (MQ=255) tCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGcaca > 1:1131089/1‑75 (MQ=255) tCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGcaca > 2:526773/1‑75 (MQ=255) cTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACaa < 1:1287611/75‑1 (MQ=255) ttACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACaaaa > 1:1345345/1‑75 (MQ=255) cAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAg < 2:1360641/75‑1 (MQ=255) aataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAgg < 1:1462270/75‑1 (MQ=255) aataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAgg > 1:1582133/1‑75 (MQ=255) ataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTAGCCACGGCCCAACCCAGCtt > 2:304768/1‑59 (MQ=255) ataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCtt < 1:304768/59‑1 (MQ=255) aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt < 1:645924/75‑1 (MQ=255) aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt < 1:1333844/75‑1 (MQ=255) aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt > 2:1573082/1‑75 (MQ=255) atGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGtgtg < 1:515217/75‑1 (MQ=255) tGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGa < 1:338682/75‑1 (MQ=255) gATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGAt > 1:908327/1‑75 (MQ=255) tCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGTCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCg < 2:1323282/75‑1 (MQ=255) cTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCgg > 1:1299228/1‑75 (MQ=255) tACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGgcg > 1:1474165/1‑75 (MQ=255) gggAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCAt > 1:487222/1‑75 (MQ=255) ggAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATc > 2:524476/1‑75 (MQ=255) gAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcc < 2:1035980/75‑1 (MQ=255) gCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcccc < 1:1398066/75‑1 (MQ=255) gCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcccc > 2:1241956/1‑75 (MQ=255) gTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTac < 2:112418/75‑1 (MQ=255) gTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTac > 1:316547/1‑75 (MQ=255) ggCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTacac < 1:1495635/75‑1 (MQ=255) cgcgcCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGa < 1:188162/75‑1 (MQ=255) ccAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCAc < 2:1535295/75‑1 (MQ=255) ccAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCAc < 2:1158987/75‑1 (MQ=255) aCTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATc < 1:1608809/75‑1 (MQ=255) cTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCCTCCCCTACACGACCACGAATcc < 2:1315879/75‑1 (MQ=255) cTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATcc < 1:100556/75‑1 (MQ=255) ttACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCa < 2:12435/75‑1 (MQ=255) tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt < 1:429282/75‑1 (MQ=255) tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt < 1:1405231/75‑1 (MQ=255) tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt > 1:167262/1‑75 (MQ=255) cccccCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATc < 1:698254/74‑1 (MQ=255) | ATATAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATC > REL606/3259135‑3259277 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |