Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,824,380 | 2 bp→AA | 100% | pseudogene (513‑514/687 nt) | rph → | ECK3633:JW3618:b3643; ribonuclease PH |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,824,380 | 0 | G | A | 98.5% | 191.9 / ‑4.2 | 67 | pseudogene (513/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (30/36); total (31/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
* | W3110S.gb | 3,824,381 | 0 | G | A | 98.5% | 190.7 / ‑4.3 | 67 | pseudogene (514/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (30/36); total (31/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG > W3110S.gb/3824318‑3824444 || aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTGTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa > 1:3548455/1‑71 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg < 1:571632/70‑1 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg < 1:3990401/70‑1 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa > 1:321600/1‑71 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGGTATTGTAAACGGTGAAGCAATTTGCg < 1:3861988/70‑1 (MQ=255) aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt > 1:3728049/1‑71 (MQ=255) aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc > 1:4313045/1‑71 (MQ=255) aCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCt > 1:1343977/1‑71 (MQ=255) ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg > 1:4496380/1‑71 (MQ=255) aacccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg > 1:2132856/4‑70 (MQ=255) acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa > 1:3589874/3‑71 (MQ=255) cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc > 1:1737716/1‑71 (MQ=255) gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc < 1:1601840/58‑1 (MQ=255) gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt > 1:1134665/1‑61 (MQ=255) gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg > 1:1423959/1‑71 (MQ=255) aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATAc > 1:2140326/1‑69 (MQ=255) aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt > 1:3969904/1‑71 (MQ=255) aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt < 1:2433068/71‑1 (MQ=255) tGAAAGGGATGGTAGCCGCGGTCTCTGTCGGGATTGTGAACGGTGAAGCGGTTTGCGATCTGGaa > 1:925511/1‑65 (MQ=255) tGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt < 1:4399799/71‑1 (MQ=255) aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg < 1:989373/70‑1 (MQ=255) aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa < 1:1993110/71‑1 (MQ=255) aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa < 1:4453496/71‑1 (MQ=255) aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg < 1:3655848/71‑1 (MQ=255) aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAg < 1:2673858/71‑1 (MQ=255) aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct > 1:201971/1‑70 (MQ=255) tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg < 1:1255633/63‑1 (MQ=255) gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa < 1:4397379/63‑1 (MQ=255) tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc < 1:2969009/70‑1 (MQ=255) aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg < 1:3273058/44‑1 (MQ=21) ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg < 1:3709661/54‑1 (MQ=255) ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa < 1:3312934/71‑1 (MQ=255) ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa < 1:2591528/71‑1 (MQ=255) ccGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa < 1:4256024/71‑1 (MQ=255) gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc < 1:3982888/43‑1 (MQ=19) gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg > 1:3146170/1‑45 (MQ=25) gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg < 1:2066202/45‑1 (MQ=25) gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa > 1:795706/1‑69 (MQ=255) cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa > 1:4011926/1‑68 (MQ=255) aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg < 1:4030100/38‑1 (MQ=16) aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa > 1:3578420/1‑46 (MQ=25) gTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa < 1:1526810/66‑1 (MQ=255) tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa < 1:4288044/53‑1 (MQ=255) tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc < 1:1725134/56‑1 (MQ=255) ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa > 1:3666359/1‑43 (MQ=19) tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa > 1:1048061/1‑44 (MQ=21) tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga > 1:1660460/1‑65 (MQ=255) tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAccta > 1:3942372/1‑69 (MQ=255) tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa > 1:3277817/1‑71 (MQ=255) cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa < 1:4522083/70‑1 (MQ=255) tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAc < 1:3965933/50‑1 (MQ=255) tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg < 1:3782366/71‑1 (MQ=255) tCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg < 1:773061/71‑1 (MQ=255) gtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc < 1:538908/59‑1 (MQ=255) taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga > 1:1362981/2‑57 (MQ=255) taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg < 1:829932/70‑1 (MQ=255) aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc < 1:3675257/47‑1 (MQ=37) aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt > 1:2758869/1‑65 (MQ=255) aTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg > 1:2579504/1‑70 (MQ=255) ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATg < 1:3587321/65‑1 (MQ=255) ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg > 1:1788059/1‑69 (MQ=255) aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatg < 1:1483511/71‑1 (MQ=255) aCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat > 1:1049470/1‑70 (MQ=255) gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa < 1:1668691/42‑1 (MQ=25) gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg < 1:3152273/58‑1 (MQ=255) tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc > 1:1161518/2‑70 (MQ=255) gAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCg > 1:4245443/1‑70 (MQ=255) || AAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCG > W3110S.gb/3824318‑3824444 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |