Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,082,696 | 0 | G | C | 21.6% | 114.8 / 29.4 | 52 | Q22E (CAG→GAG) | fdoH | formate dehydrogenase‑O, Fe‑S subunit |
Reads supporting (aligned to +/- strand): ref base G (21/19); new base C (5/6); total (26/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.43e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCAT > NC_000913/4082553‑4082830 | ggggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcg > 8:22053/1‑149 (MQ=255) gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg < 1:145090/149‑1 (MQ=255) ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg < 7:22053/149‑1 (MQ=255) tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga < 5:139567/149‑2 (MQ=255) cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGCCCCGCGCCTCAGGCGCgggacttc < 5:152664/143‑6 (MQ=255) cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGCCCCGCGCCTCAGGCGCgggacttc > 6:152664/1‑138 (MQ=255) cgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg < 5:1124/112‑1 (MQ=255) cgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg > 6:1124/1‑112 (MQ=255) aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata > 2:9447/1‑149 (MQ=255) cGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtc < 1:9447/149‑1 (MQ=255) gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCCCGGAACGACGAATGATATCTTGAGATATATAACCCAAATgtcg > 6:95716/1‑149 (MQ=38) cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGaa < 3:112864/109‑1 (MQ=255) cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGaa > 4:112864/1‑109 (MQ=255) aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCatat < 3:18787/126‑1 (MQ=255) aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCatat > 4:18787/1‑126 (MQ=255) gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTcc > 8:114624/1‑148 (MQ=255) aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac > 3:129008/1‑142 (MQ=255) aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac < 4:129008/142‑1 (MQ=255) gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt < 6:105466/54‑1 (MQ=255) gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt < 6:103077/54‑1 (MQ=255) gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt > 5:105466/1‑54 (MQ=255) gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt > 5:103077/1‑54 (MQ=255) tttCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCg < 5:95716/149‑1 (MQ=255) gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtct > 3:110748/1‑101 (MQ=255) gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtct < 4:110748/101‑1 (MQ=255) ccACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt < 3:130961/75‑1 (MQ=255) ccACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt > 4:130961/1‑75 (MQ=255) ccACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc > 4:54444/1‑149 (MQ=255) cACTTCTTGCTGTAAGTCCCGCGCGTTTGGCGCGGGTGTCAGACCGTTAGTCGCGGAACGACGAATTATATCTTTAGATTTATAAGcc < 2:139222/88‑1 (MQ=255) cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc > 1:139222/1‑88 (MQ=255) cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata > 5:110420/1‑68 (MQ=255) cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata < 6:110420/68‑1 (MQ=255) cttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt < 3:54444/148‑1 (MQ=255) cTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg > 4:92418/1‑149 (MQ=255) tGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaagga > 6:50186/1‑118 (MQ=255) tGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaagga < 5:50186/118‑1 (MQ=255) ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg > 8:150988/1‑39 (MQ=255) ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg < 7:150988/39‑1 (MQ=255) ctgacccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc > 2:106355/7‑148 (MQ=255) gAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCAt < 7:114624/149‑1 (MQ=255) aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAaggaagga > 3:42801/1‑119 (MQ=255) aaGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAaggaagga < 4:42801/119‑1 (MQ=255) ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt > 6:121032/2‑104 (MQ=255) ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt < 5:121032/103‑1 (MQ=255) cccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt > 1:91736/1‑136 (MQ=255) cccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt < 2:91736/136‑1 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg < 3:30077/57‑1 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg > 4:30077/1‑57 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt < 6:60222/79‑1 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt > 5:60222/1‑79 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaag > 7:113219/1‑112 (MQ=255) ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaag < 8:113219/112‑1 (MQ=255) | GGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCAT > NC_000913/4082553‑4082830 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |