Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 47.7% | 44.6 / 33.3 | 44 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (23/0); new base C (0/21); total (23/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.97e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTT > REL606/2296254‑2296402 | gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca < 1:1001907/76‑1 (MQ=255) tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg > 1:1545511/1‑76 (MQ=255) aTTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGc > 1:1709331/1‑76 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc > 1:1666837/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCATGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 1:1319173/1‑76 (MQ=255) tGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccacta > 2:1838253/1‑76 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1006971/1‑76 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 1:437964/1‑76 (MQ=255) ttttGCCGATGCAAGCGTGACCCCGCTTTCTCCCTTAATTGGAAGCGCATAAATTCCACCAGCCACTACCACTAAt < 2:660891/76‑1 (MQ=255) ttttGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAAt < 2:1658721/76‑1 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAAt < 1:1391434/76‑1 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAAt > 1:363058/1‑76 (MQ=255) ttcccACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAg < 2:1454854/73‑1 (MQ=255) tGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGt > 1:736391/1‑76 (MQ=255) cgCCCCTACTGCCGTGACCCCGCTTTTCCCCTTCATTGGCAGCGCATAGTTTCCCCCAGCCACTACCACTAATAGt < 2:907078/75‑1 (MQ=255) gCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTg < 2:79088/76‑1 (MQ=255) gCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTg > 1:1160332/1‑76 (MQ=255) ccACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgc < 1:1737150/76‑1 (MQ=255) ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgc > 2:944989/1‑76 (MQ=255) cACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCTACTTGCCCCTACCCCTAATAGTgcg > 2:995597/1‑76 (MQ=255) cctactgGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgcgc < 2:537895/75‑1 (MQ=255) actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc > 1:1035463/1‑76 (MQ=255) actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc > 1:1035451/1‑76 (MQ=255) tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc > 2:46794/1‑76 (MQ=255) actaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCa < 1:1375925/76‑1 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 2:6204/1‑76 (MQ=255) ctaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAg > 1:723653/1‑76 (MQ=255) ggCGTGACTCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa < 2:1667075/75‑1 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:321328/1‑76 (MQ=255) gCGTGCCCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:896270/76‑1 (MQ=255) gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:592093/76‑1 (MQ=255) gTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt < 1:783074/76‑1 (MQ=255) cacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa > 1:602545/1‑76 (MQ=255) cTTTTCCCCTTCATTGGCGGCGCATAGTTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa < 1:1454174/76‑1 (MQ=255) tttCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt < 1:374578/76‑1 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCAATACCACTAATACTGCGNCCAGTAGATAAAATGTGAAtt > 2:295126/1‑76 (MQ=255) cTTCATTTGCAGCGCATAGATNCCCCCAGCCACTACCAGTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 2:293971/76‑1 (MQ=255) ttCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 2:1545511/76‑1 (MQ=255) cATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt > 2:1392432/1‑76 (MQ=255) tggCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg < 1:1838253/74‑1 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 1:409929/1‑76 (MQ=255) gCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt < 2:1731807/76‑1 (MQ=255) cgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg < 1:1022862/76‑1 (MQ=255) cgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg < 1:839429/76‑1 (MQ=255) gcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc < 1:378326/76‑1 (MQ=255) ccATATATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc < 2:1653154/75‑1 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:602191/1‑76 (MQ=255) aGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg < 1:1762263/76‑1 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:1004217/1‑76 (MQ=255) ccaccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGtt > 1:520476/1‑76 (MQ=255) | GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTT > REL606/2296254‑2296402 |
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: ‑ |