Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,002,941 | G→A | intergenic (‑140/+45) | aspC ← / ← ompF | aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,002,941 | 0 | G | A | 96.7% | 225.1 / ‑2.3 | 60 | intergenic (‑140/+45) | aspC/ompF | aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F) |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (27/31); total (28/32) | |||||||||||
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.99e-01 |
AACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGC > REL606/1002870‑1003012 | aaCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGAcc > 1:951507/1‑76 (MQ=255) aaCTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGAcc < 1:26924/76‑1 (MQ=255) cTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCaa > 2:1365370/1‑76 (MQ=255) tATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCaaa > 2:279262/1‑76 (MQ=255) aTTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAg > 1:868288/1‑76 (MQ=255) aTTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAg > 1:707534/1‑76 (MQ=255) tgatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTg > 1:642584/1‑76 (MQ=255) gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGt > 1:912392/1‑76 (MQ=255) gatgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGt > 1:864295/1‑76 (MQ=255) tgaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGttt < 2:1207632/76‑1 (MQ=255) gaAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGtttt > 1:106979/1‑76 (MQ=255) aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc > 2:35988/1‑76 (MQ=255) aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc < 2:130552/76‑1 (MQ=255) aaCTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTc > 2:235114/1‑76 (MQ=255) cTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCgg < 2:657366/76‑1 (MQ=255) cTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCgg < 1:973904/76‑1 (MQ=255) cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa > 1:1494005/1‑76 (MQ=255) cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa > 1:869695/1‑76 (MQ=255) cAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCa > 1:966544/1‑76 (MQ=255) ggCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCAtt < 2:1420854/76‑1 (MQ=255) gCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCAttt < 2:154280/76‑1 (MQ=255) cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTa > 2:429362/1‑76 (MQ=255) gTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAAc > 2:1047959/1‑76 (MQ=255) tttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACa < 2:1498814/76‑1 (MQ=255) gttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACa < 2:51420/75‑1 (MQ=255) ttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAag < 2:1249318/76‑1 (MQ=255) tGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAaga > 2:1276888/1‑76 (MQ=255) cAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAgg > 1:1251591/1‑76 (MQ=255) aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt > 1:949834/1‑76 (MQ=255) aaGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt > 1:1127117/1‑76 (MQ=255) aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg < 2:288141/76‑1 (MQ=255) aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg < 2:1040779/76‑1 (MQ=255) gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGATgtgt < 2:1431193/76‑1 (MQ=255) gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt < 2:1227455/76‑1 (MQ=255) gACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtgt < 2:705597/76‑1 (MQ=255) aCGTGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGtgtg < 1:327866/76‑1 (MQ=255) tGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTa > 2:920836/1‑76 (MQ=255) tGAGATTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTa < 1:982541/76‑1 (MQ=255) agaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAtt < 1:429362/76‑1 (MQ=255) aTTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg > 1:596439/1‑76 (MQ=255) aTTGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg < 1:1467789/76‑1 (MQ=255) tGCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa < 2:1251138/76‑1 (MQ=255) gCTCTGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc > 1:481939/1‑76 (MQ=255) ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt < 2:1068025/76‑1 (MQ=255) ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt < 2:1065100/76‑1 (MQ=255) ctctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt < 2:840681/76‑1 (MQ=255) ctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg > 1:109379/1‑76 (MQ=255) ctGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg > 2:760853/1‑76 (MQ=255) tGGAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGt < 1:514262/76‑1 (MQ=255) ggAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTa < 2:375838/76‑1 (MQ=255) gAAGGTATAAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa < 2:823920/76‑1 (MQ=255) tataAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGAt > 1:310972/1‑76 (MQ=255) tataAAAAAAACAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGAt > 1:1207612/1‑76 (MQ=255) aaaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAcc < 1:643014/76‑1 (MQ=255) aaaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAcc < 1:9324/76‑1 (MQ=255) aaaaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAccc > 1:975822/1‑76 (MQ=255) aaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcac < 2:671941/76‑1 (MQ=255) aaaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcac < 2:842212/76‑1 (MQ=255) aaaaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca < 1:1140833/76‑1 (MQ=255) aaCAAGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACAATACCCACAGc < 2:794365/76‑1 (MQ=255) | AACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGC > REL606/1002870‑1003012 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |