Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,747,903 | 0 | G | T | 21.7% | 188.5 / 33.4 | 69 | A8S (GCA→TCA) | ydiL | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (28/26); new base T (15/0); total (43/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.95e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG > REL606/1747829‑1747976 | agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc > 2:484820/1‑76 (MQ=255) agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc > 2:1257852/1‑76 (MQ=255) ggCTGATAATGTTATTTATATTTTATTNTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa > 1:1523776/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:1086334/1‑76 (MQ=255) gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:980955/1‑76 (MQ=255) gCTGATAATGTTAGTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc > 1:110161/1‑76 (MQ=255) cTGATAATGTTATTTATGTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCTAGGCACt > 1:83922/1‑76 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:582152/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc > 1:944467/1‑75 (MQ=255) cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt > 1:575343/1‑76 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgg > 1:981315/1‑74 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt > 1:965065/1‑75 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt > 1:935996/1‑75 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt > 1:1086415/1‑75 (MQ=255) tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTc > 1:218527/1‑76 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt > 1:749422/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccag > 1:475212/1‑73 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc > 1:1220324/1‑76 (MQ=255) gATAATGTTATTTATATTTTATTGTTTATGTCATGCTGCCATTACACGATGAACGCTNATGAACTCCAGGCACTcc < 1:1465814/76‑1 (MQ=255) tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc > 2:624646/1‑76 (MQ=255) aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc > 1:280068/1‑76 (MQ=255) aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc > 2:867178/1‑76 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCNCTCCGCCat < 1:39095/76‑1 (MQ=255) tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat < 1:1209193/76‑1 (MQ=255) ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat > 2:959818/1‑76 (MQ=255) atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt < 2:366464/76‑1 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:164638/76‑1 (MQ=255) tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt < 2:536421/76‑1 (MQ=255) ttatATTTTATTGNTTATTTCATGCAGCNATNACACGATGAACGCTTAAGAACTCCAGGCACTCCGCCATAttttt < 2:1381228/76‑1 (MQ=255) tataCTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg > 1:1240786/1‑76 (MQ=255) atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa < 2:939926/76‑1 (MQ=255) ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt > 2:269954/1‑76 (MQ=255) tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg < 2:1421847/76‑1 (MQ=255) ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa < 2:285521/76‑1 (MQ=255) ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa > 1:476432/1‑76 (MQ=255) ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa < 1:999935/76‑1 (MQ=255) tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc > 2:387990/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGACCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt < 1:1221295/76‑1 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCNATGACTAt > 1:1339096/1‑76 (MQ=255) tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt > 1:702122/1‑76 (MQ=255) gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt > 1:400968/1‑76 (MQ=255) gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt < 2:837446/76‑1 (MQ=255) tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg > 1:1495374/1‑76 (MQ=255) tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg < 1:969059/76‑1 (MQ=255) ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga > 2:515786/1‑76 (MQ=255) tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat > 1:366479/1‑76 (MQ=255) tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat > 1:1172606/1‑76 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga < 2:463413/74‑1 (MQ=255) aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga > 1:463413/1‑74 (MQ=255) tttCATGCAGCCATTACACGATGACCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga < 1:652535/76‑1 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:1275830/1‑76 (MQ=255) ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa > 1:676782/1‑76 (MQ=255) tCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt < 2:757791/76‑1 (MQ=255) cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg < 1:1324440/76‑1 (MQ=255) cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg < 1:588323/76‑1 (MQ=255) ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGANTGCGCGAcc < 1:37863/76‑1 (MQ=255) ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc > 1:1028928/1‑76 (MQ=255) aTTACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg < 2:824930/76‑1 (MQ=255) aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg < 1:384050/76‑1 (MQ=255) aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg < 2:141189/76‑1 (MQ=255) cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt < 2:70685/76‑1 (MQ=255) cacGATGAACGCGTATGAACTCNANGCACACCGCCNTATTTTTGCCATGANGATAGATGAATGCGCGTCCGGGaat < 2:27698/76‑3 (MQ=255) aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc > 1:1258361/1‑76 (MQ=255) tGACCGTTTATAACCCCCGGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa < 1:343714/76‑1 (MQ=255) agcGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa < 2:719727/74‑1 (MQ=255) aaCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa > 2:686731/1‑76 (MQ=255) cGCTTATGACCCCCGGGCGCTGCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc < 1:1453701/76‑1 (MQ=255) gCTTATGAACTCCAGGCACTCCGCCACATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg > 1:918001/1‑76 (MQ=255) ttATGACCCCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt < 1:796339/76‑1 (MQ=255) tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg > 1:446277/1‑76 (MQ=255) gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATGGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa < 1:379950/76‑1 (MQ=255) tCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa > 1:550439/1‑76 (MQ=255) cccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa < 1:844379/75‑1 (MQ=255) cAGNCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGNGACCTGGATTGCCCAGACCGGTGACAGTGaaa < 2:1434239/76‑1 (MQ=255) cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa > 2:343909/1‑76 (MQ=255) aGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg < 2:867354/76‑1 (MQ=255) | AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG > REL606/1747829‑1747976 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |