Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 3,737,443 | 0 | T | C | 36.1% | 46.1 / 17.0 | 38 | I229I (ATT→ATC) | paoB | Aldehyde oxidoreductase FAD‑binding subunitPaoB |
Reads supporting (aligned to +/- strand): ref base T (19/4); new base C (0/13); total (19/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTT > ACB143/3737316‑3737572 | ttGCGGCATACCATCGCCCCCCCCAGGACACCCCGGAGCTACAAAATCAGCTGGAAATCTATCAAGTGGTGGGCGCAATCGAGCTGCCGCCCCCGGTGTTCCAAAAACATTACGTTTATATTAAAATTCGCGAcc < 4:235607/135‑1 (MQ=255) gCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACcgcg > 3:357047/1‑136 (MQ=255) ccATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTAc > 8:328721/1‑136 (MQ=255) cGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTg > 1:72337/1‑136 (MQ=255) ccGGACCGATACATTCCGCTGCCATTCGATTAACTGGTGGGCGCAATGGAGGGGCCGCCCCCGTGGTCCAAAAAACAATACCTTAACCTTAAAACTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTcg < 8:1312/136‑1 (MQ=255) ggAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTcgcc > 3:29969/1‑136 (MQ=255) gCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTc > 6:309486/1‑136 (MQ=255) ccAATCGATGAAGTGGTGGGCGCCATCGAGCTGCCGCCCCGGGTGTTCTAAAAAAATTACCTTTATCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATg < 5:99860/135‑1 (MQ=255) cAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTACTACGCGTTTGCGCTGTTATCTGTCGCCGCAGCCTTGAGCTTCGATgg > 2:97010/1‑136 (MQ=255) gATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAgg > 7:15017/1‑136 (MQ=255) tGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGAGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTAGACCTTCGATGGCCAGGt > 8:151237/1‑135 (MQ=255) aaCTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGAt > 1:240486/1‑136 (MQ=255) gCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCACGCTc > 2:267353/1‑136 (MQ=255) gCCAGCGTGCTGCCCCCCCCGGTTTAATAAAAAAATAACCTTTAACTAAAATTACGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTc < 7:279741/136‑1 (MQ=255) cATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGc > 6:346778/1‑136 (MQ=255) cATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGc > 4:356166/1‑136 (MQ=255) cgcgcggcccccCCCGGCGTACCAAAAACAAACCCGTAACCAAAAAACCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTg < 3:270239/130‑1 (MQ=255) cGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTg > 5:297112/1‑136 (MQ=255) ccccgGGTTTACAAAAAACATAACCTATATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcg < 6:21612/136‑1 (MQ=255) cccgGGTGTCCTAAAAACACACCCGTTCTCTTACATCCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcgg < 3:221152/136‑1 (MQ=255) ccccGGTGTTCCAAAAACATAACCGTTATCTCAAAAACCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGggcgg < 6:126169/136‑1 (MQ=255) ttCCAAAAAAATACCCTTAATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCt < 4:357047/135‑1 (MQ=255) aaaaaCATACCCTTAATCTAAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACaa < 3:66686/136‑1 (MQ=255) aaaaCATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAg > 4:318072/1‑136 (MQ=255) ataaCCTTTATCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTg < 2:19818/133‑1 (MQ=255) aTTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTg > 1:88209/1‑136 (MQ=255) aCCTTTATATCAAAAAACGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg < 3:132636/136‑1 (MQ=255) aCCTTTAGCTGAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg < 5:122028/136‑1 (MQ=255) aCCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGc > 8:341605/1‑135 (MQ=255) ccGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCg < 8:222040/135‑1 (MQ=255) gtttttctcAAATTCCGCGACCGCGCCTCCTCCGCGTTTGCGCTGGTATCGGGCGCCGCAGCCTTGAGCTTCGATGGCCATGGGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag < 5:60312/131‑1 (MQ=255) cttaaTCTCAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag < 4:318410/132‑1 (MQ=255) cGTTTTCTTAAAATCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag < 3:97638/136‑1 (MQ=255) cGTTATCTCAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCgag < 7:328721/136‑1 (MQ=255) tCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCg > 1:161943/1‑136 (MQ=255) tGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGt > 3:357709/1‑136 (MQ=255) tGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGt > 2:346369/1‑136 (MQ=255) aaaaTCCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGtt < 4:203686/135‑1 (MQ=255) | TTGCGGCATTCCATCGCCTCCCCGAGGACACCCCGGAGCGAGACAATCAGCTGGCAATCGATGAACTGGTGGTCGCCATCGAGCTGCCGCCCCCGGTGTTCGAAAAACATTACCGTTATCTGAAAATTCGCGACCGCGCCTCCTACGCGTTTGCGCTGGTATCGGTCGCCGCAGCCTTGAGCTTCGATGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTT > ACB143/3737316‑3737572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |