| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AE015451_phaM1 | 5,223,013 | 0 | C | A | 42.9% | ‑6.4 / 24.9 | 35 | L428F (TTG→TTT) | PP_4603 CDS (transposase) | |
| Reads supporting (aligned to +/- strand): ref base C (1/19); new base A (0/15); total (1/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG > AE015451_phaM1/5222912‑5223140 | gggagtgCTGGGTGATTTCTTATCAAACAAGCTGCGACCCA‑TTGTTCTCTAAGCACTC‑A‑CGCGGCCACCAGTT‑GATACCCCAGTTTCCTGGCTTTGCGAAACAAAGCCCGGACCTGCCGGTCTTTGCTTCTGGTTTc < 2:2891880/131‑1 (MQ=1) ggttGGGTGATCTTTCATCAAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCTGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAAct < 1:1702106/134‑1 (MQ=0) cgggttGATCTCTCATCAAACAAAATGAAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAActcc < 1:11859020/132‑1 (MQ=0) tctccccaaaaaaaATGCACCCCA‑TTTTTTTCGAAGCACTA‑A‑AGTGGCCACCATTG‑GAGACCCCATTTTCCGGGTTTTGCGAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 1:1549291/131‑1 (MQ=1) tcgcatcaAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GAGACCCCATTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 2:13976992/134‑1 (MQ=0) tcgcaccaaacaaaATGCAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCCGGGCTTTGCGAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 2:5811174/131‑1 (MQ=0) tcccatcaAACAAAAGACAACCCA‑TTTTTTCCTAAACACTA‑A‑AGCGGCCACAGGTT‑GGCACCCCATTTTCTGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 1:9140338/134‑1 (MQ=1) ccccaccaaacaaaAGACAACCCA‑TGGTTTTCTAAGCACCC‑A‑AGTGGCCACCAGTT‑GATACCCCAGTTTCCGGGTTTGGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 1:6039184/131‑1 (MQ=1) cccaacaaaaaaaAAGAAAACCAGTTTTTTTCGAACCACTA‑A‑AGAGGCCACAGGTG‑GATACCCCGTTTTCTGGGTTTGGAAAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt < 1:689980/132‑1 (MQ=0) tcctcaCACAAAATACAACCCA‑TTGTTTTCTCAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGt > 2:9513128/4‑137 (MQ=0) aagtGCAACCCG‑TTGTTTTCGAAGCACTCTAGAGCGGCCACCGTTT‑GATACCCCGTTTTCCGGGCTTTGAGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATa < 1:6577975/134‑1 (MQ=1) tgaaaCCCA‑TTTTTTTCTAAGCACCC‑A‑AGCGGCCACCAGTG‑AATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCt < 1:7945393/134‑1 (MQ=1) tttttttCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCATTTTCCGGGTTTGGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg < 1:15193842/134‑1 (MQ=0) tttttttCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCAATTTTCTGGGCTTTGCAAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg < 1:4614453/134‑1 (MQ=1) tttttttCTAAGCACCC‑A‑AGCGGCCACCATTG‑GATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAATCCTGCTTCTTgg < 1:9323057/134‑1 (MQ=1) tttttttCGAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCATTTTCTGGGTTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg < 2:14169810/134‑1 (MQ=1) tttttttCGAAGCACCA‑A‑AGCGGCCACCAGTG‑GATACCCCATTTTCCGGGCTTGGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg < 1:4614499/134‑1 (MQ=1) tttttctCTAGGCACTC‑A‑AGCGGCCACCAGTT‑GAGACCCCATTTTCCGGGCTTTGAGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg < 2:5890914/131‑1 (MQ=1) ttttAAGCACTC‑A‑AGCGGCCACCAGTT‑GATCCCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAAc < 2:1668616/134‑1 (MQ=1) tcgaaGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCGTTTTCTGGGTTTGGAGAAAAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 1:4772237/134‑1 (MQ=1) tcgaaGCACCC‑A‑AGCGACCACAGGTT‑GACCCCCCGTTTTCTGGGCTTTGAGAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 1:6457846/134‑1 (MQ=1) caaaaccACTC‑A‑AGCGGCCA‑CAGTTGGATACCCCATTTTCCGGGCTGTGCGAGAAAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTGCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 1:7331908/131‑1 (MQ=1) gaaacACTC‑A‑ATCGGCCACCAGTG‑GATACCCCAGTTTCCGGGCTTGGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg < 1:5172927/133‑1 (MQ=0) gaaacACCC‑A‑AGCGGCCACCGGTG‑GAGCCCCCGTTTTCTGGTTTTGGAAAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg < 1:6096966/133‑1 (MQ=1) aaccGGCCACCAGTT‑GATACCCCATTTTCTTGGTTTTGCGAAGCAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAATCCTGCTTCTTGGTTAACAGGTTGTACAcc < 1:1060420/134‑1 (MQ=1) tggacaCCCCGGTTTCGGGGCTTGAAAAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACt < 1:12765822/132‑1 (MQ=1) atcccccATTTTCTGGGTTTTGAGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGAt < 1:4556220/134‑1 (MQ=1) actccAGTCTCCGGGCTTTGCGAAGCAAAGCCCGGATCTACCGTTTTTGGCTTCTGGTTGCGAACTCCTCAGGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgc < 2:6902743/134‑1 (MQ=1) ggttttGCGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGa < 1:5281629/134‑1 (MQ=0) tttggCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACCCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg < 1:14640908/133‑1 (MQ=0) ttggCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 2:12743762/134‑1 (MQ=1) ttggCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 1:13867702/134‑1 (MQ=1) ttggCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 2:7633015/134‑1 (MQ=1) ttggCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 1:5637951/134‑1 (MQ=1) ttggCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 1:341086/134‑1 (MQ=1) | GGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTC‑A‑AGCGGCCACCAGTT‑GATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG > AE015451_phaM1/5222912‑5223140 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |