Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,747,903 | 0 | G | T | 43.6% | 73.0 / 90.0 | 78 | A8S (GCA→TCA) | ydiL | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (20/24); new base T (34/0); total (54/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.90e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA > REL606/1747828‑1747974 | cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg > 2:1507333/1‑76 (MQ=255) ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa > 1:458616/1‑76 (MQ=255) ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa < 2:1588810/76‑1 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:62875/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:480250/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:424280/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGNTTATTTCACGCCGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:1496384/1‑76 (MQ=255) cTGATAATGTTCTTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACACTTATGAACTCCAGTCAcc > 1:46720/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg > 1:277713/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:197272/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca > 1:362361/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:1034855/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:840101/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:643963/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTCATGAACTCCAGTCACt > 1:40857/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTCTATGTCNTGCAACCNCTNCCNGATGAACGCTTATGAACTCCAGTCACt > 1:1488412/1‑76 (MQ=255) cTGATAATGTTAGTTCTATTTTATTGTTTATTTGATGCAGCCATTACACGATGAACGCTTANGAACTCCAGTCACt > 1:1469244/1‑76 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcga > 1:66542/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct > 1:1550499/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac > 1:196265/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt > 1:965936/1‑75 (MQ=255) tGATAACGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt > 1:311553/1‑75 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc > 1:70451/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc > 1:1347645/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta > 1:746356/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc > 1:721573/1‑76 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc < 2:964922/76‑1 (MQ=255) gATAATGTTATGTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt > 1:43582/1‑73 (MQ=255) gATAATGTTATCTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcact > 1:50202/1‑73 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag > 1:1125907/1‑73 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag > 1:549218/1‑73 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac > 1:1633566/1‑73 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac > 1:506619/1‑73 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcca > 1:1339632/1‑75 (MQ=255) aTAATGTTAGTTACATTTTCTTGTTTATTTCNTGCAGCCGTCACACGATGAACGCTTATGAACTNCAGTCACttag > 1:202060/1‑73 (MQ=255) aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCAGTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc < 2:1088455/76‑1 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:1353657/76‑1 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:223443/76‑1 (MQ=255) tttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata < 1:409950/75‑1 (MQ=255) gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACATGATGAACGCTTATGAACTCCAGGCACTCCGCCata > 1:850550/1‑76 (MQ=255) ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat > 2:1028120/1‑76 (MQ=255) tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatatc > 1:37564/1‑75 (MQ=255) atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt > 1:644313/1‑76 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt > 1:1078323/1‑76 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:557748/76‑1 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:661576/76‑1 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:814768/76‑1 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt > 2:710251/1‑76 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt > 1:1365461/1‑76 (MQ=255) tataTTTTATTGTTTATTTCATGCAGCCATTCCANGATGGACNCTNATNAACTACAGTCACTCCGCCATATTTTTg > 1:1490763/1‑76 (MQ=255) tataTTTTATTATTTATTTCATGCAGCCATTACACGATGAAGGCTTATGAACTCCAGGCGATCCGCCATATTTTTg > 1:57998/1‑76 (MQ=255) ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt > 2:1426231/1‑76 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg < 2:722896/76‑1 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGAGGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAtn < 1:672936/76‑2 (MQ=255) aTTGTTTATTTCATGCAGCCATTACGCGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt > 1:118912/1‑76 (MQ=255) aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt > 1:102588/1‑76 (MQ=255) ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCAGTCAGCCATATTTTTGCCATGACTATtga > 1:267363/1‑76 (MQ=255) tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat < 1:998888/76‑1 (MQ=255) tATTTCATGCAGACATTACCCGATGAACGATTAGGAACACCAGACACCCAGCCATATTTTTGCCATTGCTATtgat > 1:1659244/1‑76 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg < 2:1659242/76‑1 (MQ=255) tttCATGCAGCCATTACACGATGAACGCTTTTGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga < 1:366222/76‑1 (MQ=255) ttnCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACNATtganga < 2:1397212/76‑1 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:844881/1‑76 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:941078/1‑76 (MQ=255) ttCATGCAGCCACTACACGCTGAACGCTTATGACCTCCTGACAGACCCGCATATTTTTGCCATGACTGTTGATGaa > 1:1659569/1‑76 (MQ=255) cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg < 1:411501/76‑1 (MQ=255) cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt > 1:1312591/1‑76 (MQ=255) acacGATGAACGCTTATAAACTCCAGTCACTCCGCCATATTTTTGCCATGACTATCGATGAATGCGCGACCTGGAt > 1:1709131/1‑76 (MQ=255) cacGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt < 1:531281/76‑1 (MQ=255) acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg > 1:1071485/1‑76 (MQ=255) aTGACCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATACGCGACCTGGATTGccc < 2:580784/76‑1 (MQ=255) aTGAACGCTTATGTACTCCANGCACTCCGCCATATTNTTGCCATGACTATGGGTGAATGCGCGACCTGGATTGccc < 2:1693445/76‑1 (MQ=255) aTGAACGCTTATGAACTCCAGGTACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc < 2:354107/76‑1 (MQ=255) aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc < 2:1527633/76‑1 (MQ=255) tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa < 2:1418123/76‑1 (MQ=255) gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg < 1:1368326/76‑1 (MQ=255) aCGCTTATGACNTCCAGGCACTCCGCCATATTTTTGTCACGACTNTTGATGAATGCGCGACCTGGACTGCCCAGAc < 1:207086/76‑1 (MQ=255) tATGAACTCCAGGCACTGCGCCATAATTTTGCCATGACTATTGCTGAATGCGCGACCTGGCTTGTCCAGACCGGTg > 1:21501/1‑76 (MQ=255) tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg > 1:839921/1‑76 (MQ=255) tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg > 1:451111/1‑76 (MQ=255) tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg < 2:1354552/76‑1 (MQ=255) aTGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGa < 2:1078321/76‑1 (MQ=255) taaCTCCAGGCACTCCGCCATATTTTTGGCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa < 1:257388/75‑1 (MQ=255) aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt > 1:278646/1‑76 (MQ=255) tCCAGGCACACCGCCATGGTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCANACNGGTGACAGTGa > 2:268213/1‑76 (MQ=255) ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa > 2:961403/1‑76 (MQ=255) | CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA > REL606/1747828‑1747974 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |