| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 566,173 | 0 | C | G | 39.1% | 23.6 / 21.0 | 23 | intergenic (‑195/‑511) | intD/renD | putative integrase/protein RenD |
| Reads supporting (aligned to +/- strand): ref base C (7/7); new base G (6/3); total (13/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTCTGCGGAAGAGTGACATTTTCTTTCCTCTGTTATGCCATCACCCGCGCTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGG > NC_000913/565960‑566398 | ttcTGCGGAAGAGTGACATTTTCTTTACTCTGTTATGCTATCCCCCGCGTTCACTTGTTCTGCATGCAACGGCGACTGCAGCGCCGCCATGCTGGCTTGCAGTGTTGTGCGGAACCCGGATATTGTCGTGAGGATGCGTCAACGCCATTGCTCCCCCAATACAACCCCAATTTCAGACAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCCAACCTCGCCATCGACTTTTCAt > 1:115584/1‑245 (MQ=255) ccTCTGTTATGCCATCACCCGCGCTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCTCGCAATCACGACATAATGCAGGCCTTc < 1:69889/250‑1 (MQ=255) cACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTg < 1:78110/250‑1 (MQ=255) cTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGACGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCACTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCa > 1:125248/1‑250 (MQ=255) aCAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGAGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGt > 1:173539/1‑250 (MQ=255) gcgcCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACAAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCtt > 1:171332/1‑202 (MQ=255) cgaaagtgnTNNNGNNTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCNAATACAAAACCAATTNCANCCCGNGCCNCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTAc > 2:161475/6‑142 (MQ=255) gaaagtgattgcgcctACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGGCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTT‑cgc < 1:123009/213‑1 (MQ=255) gCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCAc < 1:108622/250‑1 (MQ=255) cGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTNCCCAAAGACNNAANCAATGNCNNNCNNNGCNTCGTCCATTTTTTCGATGANCTNNGGCACCNTCTCGTCAAAACTCGCCATGTACttt < 2:48108/141‑1 (MQ=255) cGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCCGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACttt > 1:48108/1‑141 (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) ngNGCNTCGTCNATTTNTTCGATGAACTCCGGNACCATCTCGTCAAAACTCGNNANGTACTNTTCAt < 2:89588/66‑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) | TTCTGCGGAAGAGTGACATTTTCTTTCCTCTGTTATGCCATCACCCGCGCTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGG > NC_000913/565960‑566398 |
| 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: ‑ |