Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 950,280 | T→G | intergenic (+133/‑2) | lrp → / → ftsK | DNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 950,280 | 0 | T | G | 96.8% | 232.5 / ‑1.2 | 63 | intergenic (+133/‑2) | lrp/ftsK | DNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning |
Reads supporting (aligned to +/- strand): ref base T (2/0); new base G (29/32); total (31/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01 |
TCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTA > REL606/950234‑950329 | tCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCcttgcttg < 1:1532421/51‑1 (MQ=255) ccGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGa > 1:7223454/1‑51 (MQ=255) cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGa < 1:7016021/50‑1 (MQ=255) cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAg < 1:4311424/51‑1 (MQ=255) gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGc > 1:6345755/1‑51 (MQ=255) ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGcc > 1:5147063/1‑51 (MQ=255) tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa > 1:6233808/1‑51 (MQ=255) tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa > 1:5693527/1‑51 (MQ=255) tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCa > 1:2261468/1‑51 (MQ=255) ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAg > 1:1362611/1‑51 (MQ=255) ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAg < 1:2430001/51‑1 (MQ=255) tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAgg > 1:2900950/1‑51 (MQ=255) aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGa > 1:331241/1‑51 (MQ=255) gCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGaa > 1:4886809/1‑51 (MQ=255) cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt < 1:2563388/51‑1 (MQ=255) cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt > 1:6407312/1‑51 (MQ=255) cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt > 1:7117742/1‑51 (MQ=255) cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt < 1:2746836/51‑1 (MQ=255) aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATa > 1:969801/1‑51 (MQ=255) aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATa > 1:5708003/1‑51 (MQ=255) tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATac < 1:6938219/51‑1 (MQ=255) tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATac < 1:4331214/51‑1 (MQ=255) cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAAt < 1:2743553/48‑1 (MQ=255) gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATacac > 1:6984412/1‑51 (MQ=255) gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac < 1:2530290/51‑1 (MQ=255) gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac > 1:5091955/1‑51 (MQ=255) gggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac < 1:785812/51‑1 (MQ=255) ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACt < 1:3557081/51‑1 (MQ=255) ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACt > 1:3230470/1‑51 (MQ=255) gCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTg > 1:1270137/1‑51 (MQ=255) cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa > 1:2669661/1‑51 (MQ=255) cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa < 1:894068/51‑1 (MQ=255) cAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGa < 1:2917695/51‑1 (MQ=255) aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa > 1:6105426/1‑51 (MQ=255) aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa > 1:1768534/1‑51 (MQ=255) aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa < 1:5074066/51‑1 (MQ=255) aGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGaa < 1:6621145/51‑1 (MQ=255) ggAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAg < 1:566521/51‑1 (MQ=255) ggAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAg < 1:66377/51‑1 (MQ=255) aaaGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACa < 1:2590768/51‑1 (MQ=255) aaGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACaa < 1:5675704/51‑1 (MQ=255) cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa < 1:1073189/51‑1 (MQ=255) cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa < 1:3546612/51‑1 (MQ=255) cTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaa > 1:3788703/1‑51 (MQ=255) tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATTCACTGAAGACAaagaag < 1:6000948/51‑1 (MQ=39) tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaag < 1:6263659/51‑1 (MQ=255) tGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAaagaag > 1:2137648/1‑51 (MQ=255) gTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGt > 1:6981551/1‑51 (MQ=255) aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac > 1:2438379/1‑51 (MQ=255) aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac < 1:3702726/51‑1 (MQ=255) aCCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcac < 1:3682939/51‑1 (MQ=255) ccTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcaca < 1:2886511/51‑1 (MQ=255) ccTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTcaca > 1:1473235/1‑51 (MQ=255) cTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGCAGTCACAt > 1:4601571/1‑51 (MQ=255) cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt > 1:4436778/1‑51 (MQ=255) cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt > 1:5231927/1‑51 (MQ=255) cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt > 1:4825188/1‑51 (MQ=255) cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt < 1:7238482/51‑1 (MQ=255) cTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAt > 1:1575407/1‑51 (MQ=255) tGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACAtt < 1:4830650/51‑1 (MQ=255) ggAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTg < 1:1510737/51‑1 (MQ=255) gagaGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGa < 1:6051914/51‑1 (MQ=255) gagCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACa < 1:6208542/51‑1 (MQ=255) ccTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAg < 1:4061096/51‑1 (MQ=255) cttgcttgAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGt > 1:2772943/1‑51 (MQ=255) ttCTTGAGCCAGGAATACACTGAAGACAAAGCAGTCACATTGACAAAGTTa > 1:3004462/1‑51 (MQ=255) | TCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTA > REL606/950234‑950329 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |