| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 566,205 | 0 | T | C | 38.1% | 7.4 / 16.1 | 21 | intergenic (‑227/‑479) | intD/renD | putative integrase/protein RenD |
| Reads supporting (aligned to +/- strand): ref base T (7/6); new base C (5/3); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCTCTGTTATGCCATCACCCGCGCTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGG > NC_000913/565986‑566398 | 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) 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) gTCAANACTCGNCNNGTAGTTTTCATTCCGCTCANTCACGACATAATGCANGCCTNCACGCntantn > 2:97153/1‑63 (MQ=255) | CCTCTGTTATGCCATCACCCGCGCTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTT‑CGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGG > NC_000913/565986‑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: ‑ |