Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,747,903 | 0 | G | T | 28.4% | 223.0 / 68.8 | 88 | A8S (GCA→TCA) | ydiL | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (38/25); new base T (25/0); total (63/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.93e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA > REL606/1747831‑1747974 | gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:545459/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:2070394/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:2018408/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:1981300/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg > 1:361648/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg > 1:321676/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:434624/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:2016684/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:1607781/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca > 1:769827/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca > 1:1875386/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:935365/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:747244/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:331642/1‑76 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgc > 1:1053030/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca > 1:1398085/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac > 1:387615/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc > 1:1807960/1‑76 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATAAACGCTTATGAACTCCAGTCAcca > 1:1219763/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgat > 1:1848667/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcact > 1:967114/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta > 1:1203207/1‑74 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc > 1:1377929/1‑76 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc > 1:898201/1‑76 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCCTTACACGATGAACGCTTATGAACTACAGTCACTCCg > 2:1502513/1‑76 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcca > 1:283772/1‑75 (MQ=255) aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg > 2:106753/1‑76 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaa > 1:280031/1‑74 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc > 2:1177603/1‑76 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc > 2:917942/1‑76 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc < 2:1816351/76‑1 (MQ=255) aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat > 1:1832016/1‑73 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:141614/76‑1 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:200236/76‑1 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 2:561321/76‑1 (MQ=255) ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat > 2:2068639/1‑76 (MQ=255) ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat > 1:75877/1‑76 (MQ=255) ttgttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt < 1:1171435/73‑1 (MQ=255) atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt > 1:90527/1‑76 (MQ=255) atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt > 1:1779510/1‑76 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:215645/76‑1 (MQ=255) ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt > 1:1324056/1‑76 (MQ=255) ntatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc > 1:5928/2‑76 (MQ=255) atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc > 2:1793708/1‑76 (MQ=255) atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc > 2:1984676/1‑76 (MQ=255) atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc > 2:1557016/1‑76 (MQ=255) tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc < 2:941622/76‑1 (MQ=255) tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc < 2:1214695/76‑1 (MQ=255) atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa > 1:1909456/1‑76 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg < 1:1363716/76‑1 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg > 1:816639/1‑76 (MQ=255) tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc > 2:475545/1‑76 (MQ=255) tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc > 2:1188009/1‑76 (MQ=255) tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc > 2:1837692/1‑76 (MQ=255) aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt > 2:1180272/1‑76 (MQ=255) ttttttATTTCTTGCAGCCATTACACGATGAACGCTTATGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTa < 1:1710552/73‑1 (MQ=255) ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa > 1:771236/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt > 1:1166133/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt < 2:1017059/76‑1 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt > 2:767881/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt > 1:341572/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt > 2:1858327/1‑76 (MQ=255) gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt > 1:329987/1‑76 (MQ=255) tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg > 2:1102157/1‑76 (MQ=255) ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga > 2:1748238/1‑76 (MQ=255) ttATTTCATGCAGCCATCACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga > 1:161621/1‑76 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg > 2:9742/1‑76 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg > 2:830459/1‑76 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg > 1:1395729/1‑76 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa < 2:966245/76‑1 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:1536792/1‑76 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa < 1:1807546/76‑1 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:1967652/1‑76 (MQ=255) atccAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc < 1:1205854/73‑1 (MQ=255) tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg > 1:1670330/1‑76 (MQ=255) cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt > 1:95989/1‑76 (MQ=255) cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt < 1:367104/76‑1 (MQ=255) cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt < 1:1953294/76‑1 (MQ=255) ttACACGATGAACGTTTATAAACTCCGGGCTCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg < 1:467559/76‑1 (MQ=255) ttACACGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg < 1:1137650/76‑1 (MQ=255) ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg < 1:1405177/76‑1 (MQ=255) ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg < 1:1432971/76‑1 (MQ=255) tACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa < 1:1633082/76‑1 (MQ=255) tACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa < 1:1591915/76‑1 (MQ=255) cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt < 2:1395729/76‑1 (MQ=255) cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt < 2:496559/76‑1 (MQ=255) acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg > 1:1936678/1‑76 (MQ=255) cGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGc < 2:1685980/76‑1 (MQ=255) tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa < 2:1214407/76‑1 (MQ=255) gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg < 1:49823/76‑1 (MQ=255) aCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAc < 2:1795817/76‑1 (MQ=255) cGCTTATGATCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc < 1:1355147/76‑1 (MQ=255) tGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc < 2:1346517/76‑1 (MQ=255) tGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc < 1:2007984/76‑1 (MQ=255) gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa < 2:1575060/76‑1 (MQ=255) ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa < 2:1943253/76‑1 (MQ=255) ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa > 1:1323074/1‑76 (MQ=255) | GCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA > REL606/1747831‑1747974 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |