| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,114,469 | C→A | A2S (GCT→TCT) | fpr ← | ferredoxin‑NADP reductase; flavodoxin reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,114,469 | 0 | C | A | 97.5% | 108.4 / ‑5.9 | 40 | A2S (GCT→TCT) | fpr | ferredoxin‑NADP reductase; flavodoxin reductase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (29/10); total (30/10) | |||||||||||
| 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 | |||||||||||
ACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCC > NC_000913/4114402‑4114530 | aCGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGc > 1:1447103/1‑68 (MQ=255) gagaCTAAACAGGGCGGCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGt < 2:1038785/68‑1 (MQ=255) agaCTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGtt < 1:1296675/68‑1 (MQ=255) tAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTtc > 2:242785/1‑68 (MQ=255) aCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTcc > 2:95771/1‑68 (MQ=255) aCAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTcc > 2:318729/1‑68 (MQ=255) aCAGGGCGTCGGTAAAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGATATTCTcc > 1:1026194/1‑68 (MQ=255) cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCTATCAGACATGTTTTTCTCCt > 1:284054/1‑68 (MQ=255) cAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCt > 2:1271631/1‑68 (MQ=255) aGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTg > 1:810808/1‑68 (MQ=255) gggCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGt > 2:189279/1‑68 (MQ=255) gCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTTCCCAATCTGACAAGTTTTTCTCCTGttt > 1:999720/1‑68 (MQ=255) gTCGGTCCAGTTCTGCGCTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTttg > 2:475367/1‑68 (MQ=255) gTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTttg > 1:469744/1‑68 (MQ=255) tCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTttgattga < 1:999982/68‑1 (MQ=255) gTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTg < 1:1064378/68‑1 (MQ=255) tCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGcc > 1:562539/1‑68 (MQ=255) tCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGcc > 2:1204764/1‑68 (MQ=255) tCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGcc > 1:209542/1‑68 (MQ=255) aCTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGGCTTGCCTTATc > 2:932144/1‑68 (MQ=255) aCTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATc > 2:933892/1‑68 (MQ=255) aCTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATc > 2:71274/1‑68 (MQ=255) aCTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATc > 2:1265028/1‑68 (MQ=255) aCTTTAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATc < 2:179402/68‑1 (MQ=255) ttAGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTc > 1:280045/1‑68 (MQ=255) atAGTGACTTTGCATGTTACCAAATCAGACATGTTTATCTCATGTTTTGATAGACATGCCTTATCTTc > 1:1427204/2‑68 (MQ=255) aGTGACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGt > 2:216909/1‑68 (MQ=255) gACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTg < 1:1113320/53‑1 (MQ=255) gACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTg > 2:1113316/1‑53 (MQ=255) gACTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTtct > 2:758075/1‑68 (MQ=255) cTTTGCCTGTTGCCCGATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTcc > 2:153006/1‑68 (MQ=255) cTTTGCCTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTtctca > 2:1443446/1‑67 (MQ=255) ccTGTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACa > 1:655819/1‑68 (MQ=255) cagTTAACCAATCAGACATGTATTTCTCAAGTTTTGATTGACTAGACTTAACATCGTTCTCCGAACaa > 2:202736/3‑68 (MQ=255) tGTTACCCAATCAGACATGTTTTTCTCCTGCTTTGATTGACTTGCCTTATCTTCGTCCTCCGAACaaa < 2:7916/68‑1 (MQ=255) gTTACCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACaaaa < 2:1082385/68‑1 (MQ=255) aCCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACaaaaaaa > 2:133451/1‑68 (MQ=255) aCCCAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACaaaaaaa > 1:769319/1‑68 (MQ=255) ccAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAAtt < 1:1224034/68‑1 (MQ=255) ccAATCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAAtt < 1:1204768/68‑1 (MQ=255) aaTCAGACATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTcc > 1:263901/1‑68 (MQ=255) | ACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCTGTTTTGATTGACTTGCCTTATCTTCGTTCTCCGAACAAAAAAATTCC > NC_000913/4114402‑4114530 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |