Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 50.0% | 35.7 / 21.1 | 40 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (18/2); new base C (0/20); total (18/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.35e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 | gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca > 2:943531/1‑76 (MQ=255) tCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccac > 2:534379/1‑76 (MQ=255) cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc > 1:1330106/1‑76 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 1:798247/76‑1 (MQ=255) tGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccacta > 1:1373364/1‑76 (MQ=255) tGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccacta > 2:851875/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 1:235850/1‑76 (MQ=255) aTTTCTTTTGCCACGACTAGCTTGACCCCGTTTTTCCCCTTCATTTGCAGCGCTTAGACTCCCCCAGccactacca < 2:247977/76‑1 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATATGCAGCGCATCGATTCCACCAGccactaccact > 2:85904/1‑76 (MQ=255) tCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccacta > 1:1558334/1‑76 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 1:697901/1‑76 (MQ=255) tGCCCCTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGt < 2:533716/76‑1 (MQ=255) tGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGt < 1:801274/76‑1 (MQ=255) tGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGt > 2:1089272/1‑76 (MQ=255) cccACTACTAGCGTCACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGTTTCCCCCAGCCACTACCACTAATAGTg < 1:1220864/75‑1 (MQ=255) cccctactaGCGTGACACCGCTTTTACCCTTAATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgc < 1:304682/73‑1 (MQ=255) ccACTACTAGCGTGACCCCGCTTTTTCCNTTAATTTGCCGCTCATAGCTTCCCCCAGCCACTACCACTAATAGTgc < 1:140185/76‑1 (MQ=255) cACTACTGGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATNCCCCCAGCCCCTACCACTAATAGTgcg < 1:1614127/76‑1 (MQ=255) cctactaGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTgcgc < 1:558983/75‑1 (MQ=255) tactaGCGTGCCCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc < 1:1246877/76‑1 (MQ=255) taccaGCGTGCCACCCCTTTTACCCTCCATTTGCAGCCCATAAATTCCCCCAGCCACTACCACTAATAGTGCGccc < 1:1246860/76‑1 (MQ=255) taGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGt < 1:322007/76‑1 (MQ=255) taGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGt < 1:708210/76‑1 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:737740/1‑76 (MQ=255) gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:1277383/76‑1 (MQ=255) gCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:1260802/76‑1 (MQ=255) gcgACACCGCTTTTCCCCTTCATTTGCAGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt < 2:1093645/74‑1 (MQ=255) tGACCCCCCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACGAATAGTGCGCCCAGTAGATa < 1:8516/76‑1 (MQ=255) gACACCGCTTTTNCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 1:22554/1‑76 (MQ=255) gCTTTTTCCCTTTTTTCCCAGCGCATAGTCTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg < 1:821514/76‑1 (MQ=255) gCTTTTCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg < 1:872455/76‑1 (MQ=255) gCTTTTCCCCTTCATGCGCAGCGCATAGTTTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg < 1:821532/76‑1 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt > 1:519542/1‑76 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt < 2:66986/76‑1 (MQ=255) tccccTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 2:1550788/74‑1 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 2:275334/1‑76 (MQ=255) ccTTCATTTGCAGCACATAGATCCGCCCTGCNACTACCACGAATAGNGCGCCNAGTAGATAAAATGAGAATTTAAt < 2:212267/76‑1 (MQ=255) ccTTCATTGGCAGCGCAAACATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt < 1:577540/76‑1 (MQ=255) ttCATTTGCGCCGCAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 1:759040/76‑1 (MQ=255) cATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt > 1:279355/1‑76 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 2:441196/1‑76 (MQ=255) gCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt < 1:851875/76‑1 (MQ=255) cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa > 1:977253/1‑76 (MQ=255) cAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAACTGTGAATTTAATCAGTCGGATa < 1:8635/76‑1 (MQ=255) gcgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa < 2:607201/76‑1 (MQ=255) gcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc > 1:232221/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:483599/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 1:1550329/1‑76 (MQ=255) tAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAATCAATCAGTCGGATACAGCGAt < 1:4783/76‑1 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAACATGTGAATTTAATCAGTCGGATACAGCGATGAAGt > 1:561062/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 1:541092/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 1:632474/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 < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |