| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 2,604,001 | 0 | T | G | 56.4% | 31.9 / 50.5 | 39 | L179V (TTG→GTG) | hyfD | hydrogenase 4, membrane subunit |
| Reads supporting (aligned to +/- strand): ref base T (6/11); new base G (21/1); total (27/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.12e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGG > NC_007779/2603901‑2604098 | aGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACT‑AAAGCGTTACTTATCACTCATATCGGCTCCGTTGGGGt > 3:1449806/1‑101 (MQ=255) ggcggcTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGtt < 1:712378/101‑1 (MQ=255) gctgctCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTTGGGGTGTAt > 3:1197207/1‑101 (MQ=255) gctCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCtt < 3:417763/101‑1 (MQ=255) tCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGgtgtttgtttc > 3:1268020/1‑94 (MQ=255) ccTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTgcc < 2:138562/101‑1 (MQ=255) ttGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCAc < 3:1027786/101‑1 (MQ=255) ttGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGCTGGGGGTGTATTTTGCCGCCGCCAc > 1:535587/1‑101 (MQ=255) tGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTAGCTGGCCGCCTCCACg < 4:68352/101‑1 (MQ=255) tGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGGGTAGCGGGCCGCCGCCACg < 4:93743/101‑1 (MQ=255) ttattaCCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTTGGGGTGTTTTTTgccgccgcccccccacttt > 3:1548962/1‑91 (MQ=255) taCCCGCGCGATACCGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCa > 4:15207/1‑101 (MQ=255) taCCCGAGCGATATCGCGCCGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCCCGCTGTTTTTGCa > 4:99574/1‑101 (MQ=255) gagCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCg < 3:220418/101‑1 (MQ=255) agCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGGCGGCGCCACGCTGGTTTTGGAGACCgg > 3:1403979/1‑101 (MQ=255) cGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGaa < 3:1326294/101‑1 (MQ=255) aGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGGTGGGGGTGTATCTTTCCGCCGCCACGCTGTTTTTGCAGACCGGAACGCTTg > 1:997484/1‑101 (MQ=255) cgcAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTgcg < 1:149983/101‑1 (MQ=255) cgcAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGGTGGGGGTGGGTTTTGCCGCCGCCCCGCCGTTTTTACCGACCGGGATGTTTgag > 1:630977/1‑99 (MQ=255) gTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAgc > 1:702307/1‑101 (MQ=255) ttCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAgcg < 3:1557708/101‑1 (MQ=255) cAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTTCCGCCGCCCCGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGcgcg > 1:1212114/1‑101 (MQ=255) cAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGGCGCCGCCACGCTGTTTTTTCAGACCGGAACGTTTGCGCCTAGcgcg > 3:1181025/1‑101 (MQ=255) cAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CCTTGGGGGTGGATCTTTCCGCCGCCACGCTGTTTTTTCTGTCCGGCACGTTTTGGCTTTTcgcg > 1:593341/1‑101 (MQ=255) aGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGa > 3:925874/1‑101 (MQ=255) cACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATg < 2:99680/101‑1 (MQ=255) cACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATg < 2:15213/101‑1 (MQ=255) aCTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGa > 1:253124/1‑101 (MQ=255) aCTAAAAGCGATACTTATCACTAATATCGGCT‑CGTTGGGGGTGTATCTTGCCCCCGCCACGCTGTTTTTTCAGACCGGAACGTTTGCGCTTAGCGCGATGa > 1:1589003/1‑101 (MQ=255) cTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGGGCTTTGCGCGGTGAg > 1:1073454/1‑101 (MQ=255) tAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCCCGCTGTTTTTGCCGACCGGAACGTTTGCGCTTAGCGCGATGAGc > 3:1144056/1‑101 (MQ=255) cGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGAGGTTTACGCGATGAGCGAGGTa > 3:654407/1‑101 (MQ=255) ttACTTATCACTCATATCGGCT‑CGTTGGGGGTGTATCTTTCCGCCGCCCCGCCGTTTTTTCAGACCGGAACGCTTGCGCTTTGCGCGATGAGCGAGGTaca > 3:1073944/1‑101 (MQ=255) tACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGGTacac > 1:601543/1‑101 (MQ=255) tACTTATCACTCATATCGGCT‑CCTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGCTacac > 3:110086/1‑101 (MQ=255) tatCGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCCCGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGtt > 3:200684/1‑101 (MQ=255) atCGGCT‑CGTTGGGGGTGTATCTTGTCGCCGCCCCGCCTTTTTTTTACACCCGAAAAATTGGGCTTTGCGCGAATAGGGAGGGTTACGGCGACCCAccgtc > 1:296873/1‑97 (MQ=255) cGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATc > 1:1344054/1‑101 (MQ=255) cGGCT‑CGTTGGGGGTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATc > 1:1113996/1‑101 (MQ=255) gCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTg > 1:1434079/1‑101 (MQ=255) gTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGttt < 1:293398/101‑1 (MQ=255) gggTTGTATCTTGCCGACGCCACGCTGTTTTTGCAGACCGGAACGTTTGCCCTTAGCGCGATGAGCGAGTTACACGGGGCCGCACGTTATCTGGGTTATgg > 1:1244123/1‑101 (MQ=255) | AGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCT‑CGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGG > NC_007779/2603901‑2604098 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |