Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 757,280 | 0 | G | T | 27.3% | 31.2 / 12.2 | 22 | A458A (GCG→GCT) | sdhA | succinate dehydrogenase, flavoprotein subunit |
Reads supporting (aligned to +/- strand): ref base G (5/11); new base T (6/0); total (11/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.24e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGCGTCTGAAAAATGCCCGTCTGGATGACACTTC > NC_000913/757141‑757414 | gTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTaaa < 2:20804/149‑1 (MQ=255) caAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATgcagca > 5:168850/1‑149 (MQ=255) agTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTATCCGTAAAACTCAACAAAAAAGAATGCAGCCTAAc > 4:23103/1‑149 (MQ=255) aTCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTATCCGTAAAAATAAAAAAAAATGTATGCAGCATAACTtctc > 6:32880/1‑149 (MQ=255) ccgctcAGGTCGCACTGCGAGATGCCAGCGAGTCTGATGTTGATGCGTCTCGTGATCGCCTGAACCGCTGGAACAATAATCGTAACGGGGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTc < 3:23103/144‑1 (MQ=255) cgcACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAg < 5:32880/149‑1 (MQ=255) gcgATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCg < 5:114930/128‑1 (MQ=255) gcgATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTATCCGTAAAAATAAAAAAAAATGTATGCAGCATAACTTCTCg > 6:114930/1‑128 (MQ=255) aGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTcc < 1:31729/127‑1 (MQ=255) aGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTATCCGTAAAAATAAAAAAAAAAATATGCAGCATAACTTCTCGGTCTTcc > 2:31729/1‑127 (MQ=255) cGAGTCTGTTGTTGATGCGTCTCTGGATCGCCTGAACCGGTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAg < 7:40412/149‑1 (MQ=255) ctGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTg < 3:38020/112‑1 (MQ=255) ctGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTCTCCGTAAAAATAAAAAAAAATGTATGAAGCATAACTTATAGGTCTTCCGTGAAGGTg > 4:38020/1‑112 (MQ=255) gATCGCCTGAACAGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGAt < 6:168850/146‑1 (MQ=255) aTCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGa < 2:150101/130‑1 (MQ=255) aTCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGa > 1:150101/1‑130 (MQ=255) gAACCGCTGGAACAATAATCGTAACGGTGTAGATCACGTTGCTCTCCGTAACAATAAAAAAAAAAGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGt > 6:172540/1‑100 (MQ=255) aaCCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGc > 7:96614/1‑145 (MQ=255) aaCCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGc < 8:96614/145‑1 (MQ=255) taatCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGCGTCTGAAAAATGccc < 4:123481/146‑1 (MQ=255) taatCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGCGTCTGAAAAATGCCCgag > 3:123481/1‑147 (MQ=255) gAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGCGTCTGAAAAATGCCCGTCTGGATGACACTTc > 5:87011/1‑149 (MQ=255) | GTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGAGCAGTTGAAAGTGATCCGCGAGCGTCTGAAAAATGCCCGTCTGGATGACACTTC > NC_000913/757141‑757414 |
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: ‑ |