Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 24.2% | 161.4 / 12.4 | 62 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (45/2); new base C (0/15); total (45/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.46e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 | gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca > 2:818767/1‑76 (MQ=255) gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca < 1:1474078/76‑1 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 1:374197/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 1:374191/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 1:1365327/1‑76 (MQ=255) aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact > 2:557285/1‑76 (MQ=255) aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacc > 2:1187390/1‑76 (MQ=255) tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac > 1:660572/1‑76 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 2:1248193/1‑76 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAAt > 1:1402565/1‑76 (MQ=255) ttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAg > 2:253519/1‑76 (MQ=255) tGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGt < 2:187849/76‑1 (MQ=255) ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACCACCACTAATAGTgc > 1:1150805/1‑76 (MQ=255) cACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTgcg < 2:1167830/76‑1 (MQ=255) cACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcg > 1:1058757/1‑76 (MQ=255) ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 2:312435/1‑76 (MQ=255) tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc > 2:1224994/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:848920/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:1078048/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:25214/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:1158952/1‑76 (MQ=255) aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa < 1:141235/76‑1 (MQ=255) aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTa < 2:268767/76‑1 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:65713/1‑76 (MQ=255) gCGTGCCCCCGCTTTTACCCTTCATTTGCACCCAATGGATCCCACCAGCCACTGCCACTAATAGTGCGCCCAGTAg < 2:1312752/76‑1 (MQ=255) gCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:945075/76‑1 (MQ=255) cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa > 1:1169248/1‑76 (MQ=255) gTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt < 2:1127149/76‑1 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCCGCGCATAGATTCCACCAGCCCCTACCACTAATAGTGCGCCCAGtatata > 2:647961/1‑76 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa > 2:677885/1‑76 (MQ=255) gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 1:292141/1‑76 (MQ=255) acacCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa < 1:1117650/76‑1 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTCCCACTAATAGTGCGCCCAGTAGATaaa > 1:799065/1‑76 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa > 1:799049/1‑76 (MQ=255) ccccGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa < 1:1245720/74‑1 (MQ=255) ccGTTTTTCCCCTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg < 1:955610/76‑1 (MQ=255) ccGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg < 2:754493/76‑1 (MQ=255) ccGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg > 2:632883/1‑76 (MQ=255) ccGCTTTTACCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg < 1:350153/76‑1 (MQ=255) cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt > 1:790/1‑76 (MQ=255) gCTTTTCCCCTTCATTTGCAGCGAATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg < 2:1319489/76‑1 (MQ=255) ttttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa > 1:769725/1‑76 (MQ=255) ttttACCCCTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa > 1:1181439/1‑76 (MQ=255) tttCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt < 2:223630/76‑1 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt > 1:1279346/1‑76 (MQ=255) tccccTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 1:1052812/74‑1 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 2:807251/1‑76 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa > 1:781710/1‑76 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGACCAGTAGATAAAATGTGAATTTa > 1:169070/1‑76 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 1:1271044/1‑76 (MQ=255) cTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 1:327857/76‑1 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 1:1503447/1‑76 (MQ=255) tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg > 1:555059/1‑76 (MQ=255) tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg > 2:939219/1‑76 (MQ=255) tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg > 2:939203/1‑76 (MQ=255) cgcATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg < 1:138162/76‑1 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:65791/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:65802/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:1025314/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 1:1244307/1‑76 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:1344988/1‑76 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:1344993/1‑76 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:1521724/1‑76 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt > 1:751903/1‑76 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt > 2:973108/1‑76 (MQ=255) caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGttt > 1:454288/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 2:1493749/1‑76 (MQ=255) | GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |