Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP006881 | 3,591,188 | A→C | 22.1% | R6R (CGT→CGG) | U712_18465 ← | Molybdenum cofactor biosynthesis protein A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP006881 | 3,591,188 | 0 | A | C | 22.1% | 51.2 / 2.5 | 34 | R6R (CGT→CGG) | U712_18465 | Molybdenum cofactor biosynthesis protein A |
Reads supporting (aligned to +/- strand): ref base A (23/3); new base C (2/6); total (25/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.44e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTCCTCTTTCTTTAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCATTTGAATCGCAAGCTCC > CP006881/3591053‑3591326 | ctcctcTTTCTTTAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCTGTCTGTAACAGAAATGAGAAGGTCTCTCAGCGGTATGTTTCTTTTATCCAACATCTGAcgc > 3:135394/1‑139 (MQ=255) ctcctcTTTCTTTAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGAcgc > 7:60924/1‑139 (MQ=255) ctttctttAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTTCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTAACCAAAATCTTAAGCGActc > 1:286805/1‑139 (MQ=255) tctttAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTg > 1:130965/1‑138 (MQ=255) gATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTa > 5:171901/1‑138 (MQ=255) agccGTGCCTAATAATGGGGCGGCCATACAGTACGTGCATCGAATGTGGCCGGGGTCTTAAACACAAATCCCAAGGTCTCTCCGCGGTCTGTTTCTTTAAACAAAAATCCCCCGCGACTCTGTAATCATAGCTACAGcc < 4:140298/137‑1 (MQ=255) ctaatattGCGGGGGGCATACCGTACTGGCATCGTATGTTGCCGCGGTCTTAAACACAAATGCGACGGTATCTCCGCGGTCTGTTTCTTTTATCCAACATCTCCCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTa < 8:1883/137‑1 (MQ=255) ccggcggGCATACCGAACGTGCAGCGGAAGTGGCCGCGGACTGAAACACAAACGCGACGGTCTCGCCGCGGGCTGTTGCTTTTACCCACCACCTGCCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGa < 7:79995/138‑1 (MQ=255) gggggggAATACAGAACGTGCATCGTATGTGGCTGCGGTTTTTAACATAAATGCGACGGTCTCTCTGCGTTGTGTTTTTTTAAACAAAAATCTCCCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGa < 8:167899/134‑1 (MQ=255) cgggggGCAAACAGTCCGTGCTTCGGACGTTGCTGGGGTCTGAAACATAAAGGCGACGGGCTCGCCGCGGGCTGTTTCTTTTACCCAACATGTCCCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGa < 4:279854/138‑1 (MQ=255) cggGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATtct > 1:223957/1‑139 (MQ=255) attttcGTGCATCGGATTTGGATGCGGTTTTTAACATAATTGCGATGGGCTCTCAGCGGTCTGTTTTTTTTAACAAACATCTCCCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGgt < 8:60924/134‑1 (MQ=255) cGTGCATCGGAAGTTGCAGCGGTCTGTACCAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGta > 4:148346/1‑139 (MQ=255) cGTGCATCGGAAGTTGCAGCGGTCTGTAACAGATATGCGAAGCTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGAAGCGACTCTGTAATCATAGCTACGGCCCCTTCCTACTCCCAGATTCTCTGAGGGTGGGta > 5:85188/1‑139 (MQ=255) cGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTTCAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGta > 3:253328/1‑139 (MQ=255) cGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGta > 2:225599/1‑139 (MQ=255) cGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGta > 6:50305/1‑139 (MQ=255) gTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGta > 3:217829/1‑138 (MQ=255) tGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGtata > 3:125588/1‑139 (MQ=255) cATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGTTGTGtataga > 5:95261/1‑137 (MQ=255) tCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATAcc > 1:88823/1‑138 (MQ=255) ggAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCTACATCTGCCGCGACTCTGTAATCATAGCTACAGCGCCTTCCTACTCCCTGATTCTCTCAGGGTGGGGATATACAt > 3:271344/1‑138 (MQ=255) ggAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGCCGCGACTCTGTAATCATTGCTACAGCCCCTTCCTACTCCCTGAATCTCTCAGGGTGTGTAAATACag > 6:195784/1‑137 (MQ=255) ggAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACAt > 2:4236/1‑138 (MQ=255) ggTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATACGtt > 6:243943/1‑139 (MQ=255) gTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTAATCCCTGATTCTCTCAGGTGTTGTATATAAATTAAATGTTGCGtt > 4:169288/1‑138 (MQ=255) tCTTTAACATAATTGCTATGGTCTCTCCGGGTTATTTTTATTTAAACAAACATCTTACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGttt < 4:135024/138‑1 (MQ=255) ctcAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGa < 7:730/76‑1 (MQ=255) ctcAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGa > 8:730/1‑76 (MQ=255) tttATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCAtt > 5:263622/1‑139 (MQ=255) ttATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCAttt > 3:201096/1‑139 (MQ=255) agacAACCTCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGCTTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCCGCCATTTGa < 4:271344/136‑1 (MQ=255) aTCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCATTTg > 5:131580/1‑138 (MQ=255) cATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCATTTGAATCGc > 7:213773/1‑138 (MQ=255) aCGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCATTTGAATCGCAAGCTcc > 7:274849/1‑139 (MQ=255) | CTCCTCTTTCTTTAAAAAAGGATAATCCGGGCCGAATAATTCGGCGGGCATACAGTACGTGCATCGGAAGTTGCAGCGGTCTGTAACAGAAATGCGAAGGTCTCTCAGCGGTCTGTTTCTTTTATCCAACATCTGACGCGACTCTGTAATCATAGCTACAGCCCCTTCCTACTCCCTGATTCTCTCAGGGTGTGTATATACATTAAATGATCCGTTTCTGACAAATCCAATCGCTGTGATATTCAATTCTTCAGCCATTTGAATCGCAAGCTCC > CP006881/3591053‑3591326 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |