| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 566,245 | 0 | G | A | 44.4% | 6.1 / 13.6 | 18 | intergenic (‑267/‑439) | intD/renD | putative integrase/protein RenD |
| Reads supporting (aligned to +/- strand): ref base G (5/5); new base A (5/3); total (10/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.64e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGC > NC_000913/566011‑566494 | cACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTg < 1:78110/250‑1 (MQ=255) cTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGACGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCACTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCa > 1:125248/1‑250 (MQ=255) aCAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGAGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGt > 1:173539/1‑250 (MQ=255) gaaagtgattgcgcctACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGGCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑cgc < 1:123009/213‑1 (MQ=255) gCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCAc < 1:108622/250‑1 (MQ=255) gcgcctACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAACACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATTACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATg > 1:227856/6‑251 (MQ=255) gTGAGGATGCGTCATCGCCCTTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGAACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCCTGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTTCTCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt > 1:156541/1‑251 (MQ=255) aTCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACAGAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGTGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc < 1:160536/251‑1 (MQ=255) ccAAATACAAAACCAATTTCAGCCAGTGCCTCGTACATTTTTTCGATGAACTCCGGCACGATCTCGTAAAAACTCGCCATGTACTTTTCATCCCGCTCACTCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑CGCGTCa > 1:201635/1‑175 (MQ=255) nnnnnncNCGNNCATTNTTTCGATGAACTCCGGCACCATCTNGTCAAAACTCGCCATGTACTTTTCATTCCGCNCNATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGNCAAAGTACCANGCANNTTTT‑NGNGNCACCCACATGCTGTACTGCACCTGGGCCANGTAAGCCGATTTTATTGCCTCTANNCCNNNGANNCGGaa < 2:173539/213‑1 (MQ=255) ccTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAAGGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTGCCAGGCATGTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATTTAGGCCGTCTTTATGGCCTCGAAACCACCGAGCCGGaactt < 1:209918/217‑4 (MQ=255) ccATTTTTTCGATGAACTACGGCACCCTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATCATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGCGTATTGCCGGAAGCGTgg > 1:80063/1‑251 (MQ=255) ngATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑CGCGTCAc > 1:141364/2‑133 (MQ=255) cGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑CGCGTCAc < 2:141364/133‑1 (MQ=255) aaCTCCGGAACGATCTCGTCAACACTCGCCATGTACTTTTCATCCCGCTCAATCCCGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAg > 1:196999/1‑111 (MQ=255) gCACGATCTCGTCAACACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAg > 1:178655/1‑104 (MQ=255) cATCTCGTCAAAACTCGCCNTNTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCNTGCGCGGGTCNTAGNNGGCANANTNCCAGGCATCTTTT‑CGCGTCACCCACATGCNGTACTGCACCTGGGCCATGTAAGNNGANNNTANNGCCTc < 2:173648/165‑1 (MQ=255) cATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGAGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTc > 1:173648/1‑165 (MQ=255) gcgcGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATACTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGGGGAAATCGAACGCGCATTACTGGTTGGTTCCCAACTTGTCCCAGAACATGCGGCCCAATGCgctggc < 1:183670/250‑1 (MQ=255) | CACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGC > NC_000913/566011‑566494 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |