Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,730,462 | 0 | C | G | 35.7% | 29.1 / 10.5 | 28 | A608A (GCG→GCC) | hypF | carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases |
Reads supporting (aligned to +/- strand): ref base C (13/5); new base G (10/0); total (23/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACTCCGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGC > REL606/2730390‑2730537 | cACTCCGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGca > 1:569540/1‑76 (MQ=255) tgcgtggcgtggGGCTGCGAGCCCCTCGAAAACCAAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAg < 2:4620721/74‑1 (MQ=255) tCCGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGcgcac > 2:2612881/1‑75 (MQ=255) cccGGGGCATGGGGCTGCGAGCGCCTCCAGAGCAAAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAg < 1:2922708/75‑1 (MQ=255) ccGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGcacaac > 2:4950863/1‑74 (MQ=255) ccGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACCGc > 2:2501659/1‑76 (MQ=255) ggggCATGAGTGGGCGAGCCCCTCAAGAGAAAAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGccc < 2:167066/74‑1 (MQ=255) gggCAGGGGGCTGGGAGCCCCTCCGGAGCAAAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCa < 1:2267895/75‑1 (MQ=255) gAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCCCTgccc > 1:4708163/1‑75 (MQ=255) ggCTGCGCGCCCCCCAAGCGCAAAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTgccgcc < 1:1952990/76‑1 (MQ=255) gCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCAGTgccgcccccc > 1:2134005/1‑75 (MQ=255) cGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCGCAGCCCCGCgccgccgcccc > 1:3879992/1‑74 (MQ=255) cGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCCGTGCCCCCGCCAc > 1:847247/1‑76 (MQ=255) gAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCAcc > 1:2335528/1‑75 (MQ=255) gAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAACCCACTGCCGCCCCCCCt > 1:2515786/1‑76 (MQ=255) aGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACACACCCGGGCCGCCGCCACTg > 2:3296597/1‑76 (MQ=255) gcCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCAGATCCGCCGCCCCTGCAt > 2:4347917/1‑76 (MQ=255) gcCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGGCCCCGCCCCTGCAt > 1:4771510/1‑76 (MQ=255) cTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCCGTGCCGCCGCCACCGCATCg > 1:2093364/1‑76 (MQ=255) cTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCCCCGCATCg > 1:2832880/1‑76 (MQ=255) tCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGGACAACCCCGGGGCACCCCCACCCCATCGa > 1:2842865/1‑76 (MQ=255) tCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACACACCACTGCCGCCGCCCCTGTCTCGa > 1:4082901/1‑76 (MQ=255) ccAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCCGCGCCGCCGCCCCCGCCTCgga > 2:3570246/1‑74 (MQ=255) gCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAAc > 2:2344851/1‑76 (MQ=255) aaaaaGCCGCTTCACTTCA‑ATAACTAACCTGGGCGGGCCCAAACCCCAGGGCCCCCGCCACTGCATCGAAAAAACg < 2:1068521/73‑1 (MQ=255) aaaaaGCCGCTCCCCCTTC‑ATACCTTAACGTGGCTGGCGCACGGCCCGGTGCCGCCGCCACTGCATCGAAAAAACg < 1:3853339/73‑1 (MQ=255) caAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAACCCAGTGCCGCCGCCACCGCATCGAAAAAACGcc > 2:1092184/1‑76 (MQ=255) aacgCCGCTTCACCTAC‑ATAACTCAACGTGGCGGGCGCCCAGCCCCGTCCCCCCGCCACTGCATCGAAAAAACGcc < 2:5118751/73‑1 (MQ=255) aaGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCGCAGCCCAGAGCCGCCGCCACCGACTCGAAAAACAGccc > 2:4304735/1‑76 (MQ=255) gCCGCTTCCCGTTC‑ATAACTTAAGGGGGCTGCCGAACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCcac < 1:5063184/76‑1 (MQ=255) gCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAAAGCCcac > 1:5107656/1‑76 (MQ=255) ccccTTCCCCTTC‑ATACCTTACCGTGGCGGGCCCACACCCCAGTCCCGCCGCCACTGCATCGAAAAAACGCCcaca < 1:582056/73‑1 (MQ=255) ccGCTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACACCCCAGGGCCGCCCCCCCCGCATCCAAAAAAAGCCccca > 1:1687031/1‑76 (MQ=255) cGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGGCCCGCGCCGCCGCCCCTGCATCGGAAAAATGCCCacac > 2:815545/1‑76 (MQ=255) cGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCcaccc > 1:1456634/1‑74 (MQ=255) cTTCACCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGa > 2:2683403/1‑76 (MQ=255) ttCACCTTC‑TCAACTAAAGGAGGCAGGCGACCACCCCAGAGCCGCCGCCACTGCATCGAAAAAACGCCCACACGAt < 1:3773157/76‑1 (MQ=255) ttCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACCGCATCGAAAAAACGCCCACACCAt > 1:207441/1‑76 (MQ=255) ttCACCTTC‑ATAACTTAACGTGGCTGGCGCAAAGCCCAGTGCCGCCGCCACCGCATCCAAAAAACGCCCACACGAt > 1:701130/1‑76 (MQ=255) tCACCTTC‑ATAACTTAACGTGGCTGGGGCCCAGCCCAGTGCCGCCCCCCCCGCATCGAAAAAACGCCCACACGATg > 2:268540/1‑76 (MQ=255) tCACCTTC‑ATAACTTAACGTGGCTGGCGCCCAGCCCAGTGCCGCCGCCCCCCCATCGAAAAAACGCCCACACGATg > 2:3354218/1‑76 (MQ=255) ccccTTCAATACCTTAAGGTGGCTGCCCCCGCGCCCGGTCCCGCCGCCACTGCATCGAAAAAACGCCCACACGATg < 2:4846809/74‑1 (MQ=255) ccccTTC‑ATAACTTAAGGGGGCTGGCGCACAGCCCGGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGa < 1:2513853/74‑1 (MQ=255) cccTAC‑ATACCTAAAGGGGGCTGGCCCACGCCCCAGTCCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGAc < 1:1916631/75‑1 (MQ=255) aCCTTC‑ATAACTTAACGTGGCTGGGGCACAGCCCAGTGCCGCCGCCACCGCATCGAAAAAACGCCCCCACGATGAc > 2:2746981/1‑76 (MQ=255) ccTTC‑ATAACTTAGCGTGGCGGGCCCCCACCCCAGCCCCCCCGCCACTGCATCGAAAAAACGCCCACACGATGACg < 2:4136633/76‑1 (MQ=255) ccTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGGCCCTGCATAGAAAAAACGCCCACACGATGACg > 1:1006776/1‑76 (MQ=255) cTTC‑ATAACTTAACGTGGCTGGGGCCCAGCCCAGTGCCCCCGCCACCGCATCGAAAAAACGCCCACACCATGaacc > 2:3752813/1‑73 (MQ=255) ttC‑AAAACTAAACGGGGCTGGCGCACAGCCCAGTCCCCCCGCCACTGCATCGAAAAAACGCCCACACGATGACGcc < 1:4568634/76‑1 (MQ=255) tac‑ATATCTAAGCGTGGCTGCCGCCGACCCCAGCCCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGcc < 2:469503/74‑1 (MQ=255) tC‑ATAACTTAACGTGGCTGGGGCCCAGCCCAGGGCCGCCGCCACTGCATCGAAAAAACGCCCCCACGATGACGCCa > 2:1975785/1‑76 (MQ=255) tC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCCCCACTGCATCGAAAAAAAGCCCACACGATGACGCCa > 1:3803182/1‑76 (MQ=255) c‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCCCCCGATGACGCCAg > 1:189594/1‑76 (MQ=255) c‑ATAACTTAACGTGGCGGGCGCACACCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAg < 1:2501659/76‑1 (MQ=255) aaCTTAACGTGGCTGGGGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACCATGACGCCCgcgg > 2:3940120/1‑76 (MQ=255) aaCTTAACGTGGCTGGCGCACAGCCCCGTGCCGCCGCCACCGCATCGAAAAAAACCCCCCACGATGAAGCCAccgg > 2:71923/1‑76 (MQ=255) cTTAACGTGGCTGGGGCACAGCCCAGTGCCGCCGCCCCTGCATCGAAAAAACGCCCACACCAAGACGCCAccgggg > 2:835747/1‑74 (MQ=255) ttAACGTGGCTGGGGCACAGCCCAGTGCCGCCCCCCCTGCATCGGAAAAAAGCCCACACGATGACGACAAGGgcgc > 2:2351304/1‑76 (MQ=255) ttAACGTGGCTGGCGCACAGCCCCGTGCCGCCGCCACTGCATCGAAAAAAAGCCCACACGATGACCCCAGCGgcgc > 2:961235/1‑76 (MQ=255) ttAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCCCTGCATCGAAAAAACGCCCACACGATGACGCCAGCGgcgc > 2:995990/1‑76 (MQ=255) tAACGTGGCTGGGGCCCAGCCCAGTGCCGCCGCCACCGCCTCGAAAAAACGCCCACACGATGACGCCACGGgcgcc > 2:2059470/1‑75 (MQ=255) aaCGTGGCTGGGGCACAACCCCGGGCCGCCGCCACTGCATCGAAAAAACGCCCACACCAAGACGCCAgcggcgggg > 2:982061/1‑73 (MQ=255) aaCGTGGCTGGGGCAAAGCCCCCTGCCGCCGCCCCGGCACCGGAAAAATGCCCACACGATGACCCCCGCGGGGCGt > 2:5185229/1‑76 (MQ=255) aaCGTGGCTGGCGCACCACCCAGTGCCCCCGCCACCGCATCGAAAAAAAGCCCACACGAAGACCCCAGCcgcgcgg > 2:555527/1‑75 (MQ=255) aaCGTGGCTGGCGCACAGCCCCGTGCCGCCGCCACCTCATCGAAAAAAAGCCCACACGATGACGCCAGCGGCGGGt > 2:4989200/1‑76 (MQ=255) gggggcgggcgCCCCGCCCGCTCCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTa < 1:4192260/73‑1 (MQ=255) ggTGGCGGCCCCACGGCCCGGCGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTa < 2:3505405/75‑1 (MQ=255) ggctggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAAtt > 2:3203493/1‑76 (MQ=255) ggctggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAAtt > 2:2859841/1‑76 (MQ=255) ggctggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAAtt > 2:2702642/1‑76 (MQ=255) ggtggcCCCCACCCCGGTCCCCCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTc < 1:5044112/74‑1 (MQ=255) tggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCa > 1:66342/1‑76 (MQ=255) tggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCa > 1:3469200/1‑76 (MQ=255) tggcGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCa > 1:2103136/1‑76 (MQ=255) ggcGCACCGCCCGGTCCCCCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCAc < 2:2134794/76‑1 (MQ=255) gcgcACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACg > 1:1080418/1‑76 (MQ=255) gcgcACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACg > 2:3454566/1‑76 (MQ=255) cgcACAGCCCAGTGCCGCCGCCACTGCATCGAAAATCAGCCCACACGATGACGCCAGCGGCGCATTAATTCCACGc > 2:1076403/1‑76 (MQ=255) | CACTCCGTGGCATGAGGCTGCGAGCGCCTCCAGAGCACAAGCCGCTTCACCTTC‑ATAACTTAACGTGGCTGGCGCACAGCCCAGTGCCGCCGCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGC > REL606/2730390‑2730537 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |