Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 23.8% | 216.6 / 22.0 | 84 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (59/5); new base C (0/20); total (59/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.95e-15 | |||||||||||
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 > 1:875459/1‑76 (MQ=255) tCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccac < 2:2076947/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 1:1627103/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 1:1865586/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTAACCTTCATTTGCAGCGCATAGATTccacca > 2:376002/1‑76 (MQ=255) tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg > 1:183018/1‑76 (MQ=255) tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg > 1:1677815/1‑76 (MQ=255) aTTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGc < 1:1630983/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc < 1:2216201/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc < 1:1274955/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc > 1:1131480/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 1:837777/1‑76 (MQ=255) aaagAATTTCTTTTGCCACTACTAACGTGACACCTCTTTTTCCCTTCATTTGCAGCGCATAGATTCCACCAGccac < 2:195103/73‑1 (MQ=255) aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccac > 1:115732/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 1:1906107/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 2:605113/1‑76 (MQ=255) aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGANTAGATTCCACCAGccactacc > 2:345088/1‑76 (MQ=255) aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACACTTCATTTGCAGCGCATAGATTCCACCAGccactacc > 1:803794/1‑76 (MQ=255) aTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacca > 2:1669299/1‑76 (MQ=255) tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac > 2:1757515/1‑76 (MQ=255) tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac > 2:43009/1‑76 (MQ=255) tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac > 1:662219/1‑76 (MQ=255) tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac > 2:119433/1‑76 (MQ=255) ttCTTTTGCCACTACTAGCGTNACCCAGCTTTTACCCTACATTTGCAGCGCATAGATTCCCCCAGccactaccact < 2:349266/76‑1 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 1:1520854/1‑76 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAAt > 2:1832990/1‑76 (MQ=255) ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAAt > 2:745763/1‑76 (MQ=255) tttGCCACTACTAGCGTGACAACGCTTTTACCTTTCATTTGTAGAGAATAGAATACAAAACCCACTAACACTAATa > 2:1215808/1‑76 (MQ=255) gCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTg < 1:1284156/76‑1 (MQ=255) ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgc > 2:933574/1‑76 (MQ=255) cACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTgcg < 2:1728062/76‑1 (MQ=255) cACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcg > 1:2214803/1‑76 (MQ=255) actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc > 2:751549/1‑76 (MQ=255) actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc > 1:2111778/1‑76 (MQ=255) tactaGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACCACCACTAATAGTGCGccc < 1:330486/76‑1 (MQ=255) actaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCa < 2:1220896/76‑1 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:104137/1‑76 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:1984023/1‑76 (MQ=255) ctaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAg > 1:766771/1‑76 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 2:240747/1‑76 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:2110701/1‑76 (MQ=255) gCGTGCCCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:756887/76‑1 (MQ=255) gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:341193/76‑1 (MQ=255) gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 2:775781/76‑1 (MQ=255) cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa > 1:108090/1‑76 (MQ=255) gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa > 1:308867/1‑76 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa > 2:1885337/1‑76 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa > 2:1593860/1‑76 (MQ=255) ccccGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa < 1:1160112/74‑1 (MQ=255) cacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa > 2:1116208/1‑76 (MQ=255) cTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa > 2:12327/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:453846/1‑76 (MQ=255) ttccccTTCATTTGCAGCGAATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt < 1:661815/73‑1 (MQ=255) ttccccTTCATTTGCAGCCCAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt < 2:1161621/73‑1 (MQ=255) tccccTTCATGTGCAGCGCAGAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 2:2196629/74‑1 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa > 1:177942/1‑76 (MQ=255) cccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa < 2:1584662/76‑1 (MQ=255) ccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt < 1:369765/76‑1 (MQ=255) cTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 2:2085682/76‑1 (MQ=255) cTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc > 1:1599336/1‑76 (MQ=255) cTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 1:710349/76‑1 (MQ=255) ttCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTCATAGTGCGCCCAGTAGATAAACTGTGAATTTAATCa > 1:946650/1‑76 (MQ=255) ttCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa > 1:946632/1‑76 (MQ=255) ttCATTTGCAGCGCAAAAATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 2:1437215/76‑1 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 1:1587971/1‑76 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 1:2027727/1‑76 (MQ=255) tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg > 2:369393/1‑76 (MQ=255) tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg > 1:281156/1‑76 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 2:1315417/1‑76 (MQ=255) cAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa < 1:2216072/76‑1 (MQ=255) aGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc > 2:1714445/1‑76 (MQ=255) aGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc < 2:1697214/76‑1 (MQ=255) aGCGCAAAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc < 1:2166346/76‑1 (MQ=255) cgcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg > 2:233009/1‑76 (MQ=255) gcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc < 2:200186/76‑1 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 2:208158/1‑76 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 1:1172677/1‑76 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 1:346213/1‑76 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 1:81311/1‑76 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 1:1319579/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:135130/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:1373334/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:1628551/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:1830344/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 2:1607273/1‑76 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:1607923/1‑76 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:423019/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTGATCAGTCGGATACAGCGATg > 1:208061/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:1061699/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:213260/1‑76 (MQ=255) gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa < 2:2096553/76‑1 (MQ=255) caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAGAATGTGAATTTCATCCNTCGGATACAGCGATGCAGttt > 2:342939/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 2:2041152/1‑76 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa > 1:792721/1‑76 (MQ=255) | GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA > REL606/2296254‑2296404 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |