| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,171,384 | 0 | A | C | 52.4% | 14.8 / 24.8 | 21 | K285Q (AAG→CAG) | insF5 | IS3 element protein InsF |
| Reads supporting (aligned to +/- strand): ref base A (7/3); new base C (4/7); total (11/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGTTCAGCAGATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGTTCCCTGTTC > NC_000913/2171131‑2171682 | cTGTTCAGCAGATCATCAAGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAAGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATATATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTCGAAAACCAGAACCTCGCTTCGGCCTGTGTCCCTATTACGTGGGTAGGATCAga > 2:781049/1‑299 (MQ=11) tGTTCAGCAGATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGTAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCNACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgag < 2:598024/300‑3 (MQ=11) gTTCAGCAGATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagt < 1:262835/300‑4 (MQ=11) gTTCAGCAGATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAAAGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAAACTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagt > 2:437456/1‑297 (MQ=2) gcagATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTAGCTTAGGCCTGTGTCCATATTAAGTGGGTAGGATCAgagtggtgt > 2:770809/1‑292 (MQ=2) gcagATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGATGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgt < 1:607157/300‑9 (MQ=2) agATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGttcggtaaaa > 2:691166/1‑291 (MQ=11) ttATCATGCGCAACTGAATCGGCATAATCTGCTTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTACGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCTGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatg < 2:33988/300‑14 (MQ=11) ttATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGTTTGCTGCTACGATAATGCCGGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCAGGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCAAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatg < 1:56498/300‑14 (MQ=11) ttATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatg < 1:4499/300‑14 (MQ=2) aaCTGAAGCTGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatgcaggccagtgg < 1:428394/300‑25 (MQ=11) ggCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTAGTGGATCATTTTGACTGCCGCCt > 1:136463/1‑300 (MQ=255) ggCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCAGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGATTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCt > 2:110115/1‑300 (MQ=255) gcgcAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGgcgcgc > 1:193957/1‑300 (MQ=255) tcttTCAGTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGGTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtct < 2:136463/300‑1 (MQ=255) acTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCNTCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGc < 2:591316/300‑1 (MQ=255) gtgtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCNCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTc > 2:592963/1‑299 (MQ=255) tAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGTCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGaaa > 1:227252/1‑300 (MQ=255) agcGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCg < 1:321775/298‑1 (MQ=255) gTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATAATCCCGCCAGATCCTGaa > 2:168879/1‑300 (MQ=255) gTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTACGAGGCAGCAAATACTCCCGCCAGATCATGaa > 2:169238/1‑300 (MQ=255) aaGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGTTCCCTGttc < 1:451712/299‑1 (MQ=255) | CTGTTCAGCAGATTATCAGGCGCAACTGAAGCGGCATAATCTGCGTGGAAGTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGTTCCCTGTTC > NC_000913/2171131‑2171682 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |