Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,738,310 | C→T | G250D (GGC→GAC) | sufD ← | component of SufBCD complex |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,738,310 | 0 | C | T | 97.6% | 153.2 / ‑3.5 | 41 | G250D (GGC→GAC) | sufD | component of SufBCD complex |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (19/21); total (20/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.88e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01 |
CGCCAGGCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGCCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCAAAGTGG > REL606/1738237‑1738381 | cGCCCGGCTATTGGTCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCa > 2:600633/1‑76 (MQ=255) cAGGCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCAccc > 1:966918/1‑76 (MQ=255) aGGCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCca < 1:1365677/76‑1 (MQ=255) aGGCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCca > 2:1074401/1‑76 (MQ=255) gCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCcagc < 1:1575231/76‑1 (MQ=255) tATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCcagcag > 1:208221/1‑76 (MQ=255) aTTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAgg < 2:1188047/76‑1 (MQ=255) ttGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGa < 2:1064353/76‑1 (MQ=255) tGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGaa < 1:1687000/76‑1 (MQ=255) gATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGaaa > 2:1693394/1‑76 (MQ=255) gATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGaaa > 2:1300765/1‑76 (MQ=255) gATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGaaa < 2:197728/76‑1 (MQ=255) tCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACt < 1:1552815/76‑1 (MQ=255) cGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGCCACCCAGCAGGAAACTATGGCTa > 2:602553/1‑76 (MQ=255) gTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTaa > 2:1385470/1‑76 (MQ=255) gCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAAt < 2:1029815/76‑1 (MQ=255) cTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATg > 1:1560184/1‑76 (MQ=255) ttCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTg < 1:721004/76‑1 (MQ=255) ttCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTg > 2:685299/1‑76 (MQ=255) tCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTgg < 2:1044733/76‑1 (MQ=255) tCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTgg < 2:873342/76‑1 (MQ=255) tCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTgg < 2:43247/76‑1 (MQ=255) cGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGAGGGGGc > 2:441537/1‑76 (MQ=255) cATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATc < 2:927075/76‑1 (MQ=255) aTTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCg < 2:1729592/76‑1 (MQ=255) tGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTc > 1:456117/1‑76 (MQ=255) tGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTc > 1:1656342/1‑76 (MQ=255) gAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTca < 1:358222/76‑1 (MQ=255) aGTTGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTcag < 2:209202/76‑1 (MQ=255) tGCGTACTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTcagcca < 1:348855/76‑1 (MQ=255) cTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGcaaca > 1:1385251/1‑76 (MQ=255) cTGGTGTTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGcaaca > 2:645348/1‑76 (MQ=255) gtgtTGTGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAAt < 2:1277890/76‑1 (MQ=255) tgttgtGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATc > 1:1634446/1‑76 (MQ=255) tgttgtGTCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATc < 2:671135/76‑1 (MQ=255) gtgtCGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTa < 2:1213681/76‑1 (MQ=255) cGTAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATggg > 2:179767/1‑76 (MQ=255) tAACACTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCa > 2:455506/1‑76 (MQ=255) cacTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTCATGGACAAAg > 2:917937/1‑76 (MQ=255) acTGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCAAAgt > 2:569072/1‑76 (MQ=255) tGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCAAAgtgg < 2:412137/76‑1 (MQ=255) tGCGTCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCAAAgtgg < 2:412124/76‑1 (MQ=255) | CGCCAGGCTATTGATCCGCAGCGTGCTGTTTTCGCCATTGAGTTGCGTACTGGTGTTGTGTCGTAACACTGCGCCACCCAGCAGGAAACTATGGCTAAATGCGGTGGCATCGTCAGCCAGCAACAAATCGTTATGGGCAAAGTGG > REL606/1738237‑1738381 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |