| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,493,045 | 0 | A | G | 27.3% | 28.2 / 11.5 | 22 | I70I (ATT→ATC) | frc | formyl‑CoA transferase, NAD(P)‑binding |
| Reads supporting (aligned to +/- strand): ref base A (8/8); new base G (3/3); total (11/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGG > NC_000913/2492904‑2493139 | tGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAAc < 3:55415/149‑1 (MQ=255) ccAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGttt > 1:139709/1‑129 (MQ=255) ccAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGttt < 2:139709/129‑1 (MQ=255) gTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAAt > 4:50323/1‑116 (MQ=255) gTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGt > 5:98005/1‑148 (MQ=255) gCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTaa > 2:110446/1‑131 (MQ=255) gCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTaa < 1:110446/131‑1 (MQ=255) aTGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATc < 6:98005/149‑1 (MQ=255) tCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCTTTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACt < 1:60283/116‑1 (MQ=255) tCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACt > 2:60283/1‑116 (MQ=255) gAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATgg > 3:138682/1‑111 (MQ=255) gAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATgg > 3:140914/1‑111 (MQ=255) gAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATgg < 4:138682/111‑1 (MQ=255) gAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATgg < 4:140914/111‑1 (MQ=255) aGATATCAGCTTCGCTTATCTGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTTGTATTTAACTAAATAGAACTTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATAtct < 7:66174/128‑1 (MQ=255) aGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCTATATAAGGAAAAtct > 8:66174/1‑128 (MQ=255) ttCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGatat > 4:124686/1‑74 (MQ=255) ttCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGatat < 3:124686/74‑1 (MQ=255) gTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATc > 4:46975/1‑62 (MQ=255) gTATTTAACTCGATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATc < 3:46975/62‑1 (MQ=255) gTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCgg < 4:37117/106‑1 (MQ=255) gTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCgg > 3:37117/1‑106 (MQ=255) | TGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGG > NC_000913/2492904‑2493139 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |