Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 48.6% | 29.3 / 19.9 | 35 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (14/4); new base A (0/17); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058869‑4059038 | aaaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAg < 1:788497/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGa < 1:831476/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGa < 1:337161/90‑1 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGaa > 2:257472/1‑82 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGaa < 1:257472/82‑1 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATc > 1:1087120/1‑90 (MQ=255) aGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATccc > 2:763449/1‑90 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 1:718687/1‑90 (MQ=255) ttcccGCGAGTCGAGAAAAAATCCGGTGAGAATATCGGAAAACGCGGAAAAAAA‑CTTCCCCACAAAAAAGATAACCCCGCCATAGTAAGc < 1:641110/87‑1 (MQ=255) cgcgAGTCGAGAAAAAACCCGGTGAGAATATGGGAAACCGTGGAAAAAAA‑CTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGc < 1:1078809/90‑1 (MQ=255) gcgAGTCGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAAA‑CTCCGCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCa < 1:763449/90‑1 (MQ=255) gcgAGTCGAGAAAAAACCGGGTGAGAATATCGGAAAACGCGAAAAAAAC‑CTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa < 2:186949/90‑1 (MQ=255) gagagtcgagaAAAAACCGGGTGAGCATACCGGAAAACGCGAAAAAAAC‑CTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa < 2:921814/88‑1 (MQ=255) aaaTCCGGTGAGCATATGGGTAAACGGGGTAAAAAA‑CTTCCCCCCCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:173896/90‑1 (MQ=255) aaaTCCGGTGAGCATATCGGAAAACGCGGAAAAAAA‑CTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:1149494/90‑1 (MQ=255) ggTGAGCATATGGGAAAACGCGGAAAAAAA‑CTCCCCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 1:820592/90‑1 (MQ=255) ggTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata > 1:67885/1‑90 (MQ=255) gCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAtt > 2:23847/1‑90 (MQ=255) cATATCGAAAAACGCGAAAAAAAC‑CTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 2:718687/90‑1 (MQ=255) aatatCGGAAAAGGGGGAAAAAAC‑CTCCCCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 2:1058533/89‑1 (MQ=255) aatatCGGAAAACGGGAAAAAAAA‑CTCCGCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 1:764676/89‑1 (MQ=255) atatCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg > 2:691192/1‑90 (MQ=255) tCGGAAACCGCGAAAAAAAA‑CTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 1:1183216/90‑1 (MQ=255) cGGAAACCGCGAAAAAAAC‑CTCCGCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 1:791014/90‑1 (MQ=255) gTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 2:576333/1‑90 (MQ=255) gTAAACGCGGTAAAAAA‑CTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:23847/90‑1 (MQ=255) aaCGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc > 1:377816/1‑90 (MQ=255) accgggTAAAAAATCTCCCCCACGAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 1:1086587/86‑1 (MQ=255) gcgGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga > 2:195457/1‑90 (MQ=255) gaaaaaaac‑CTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 2:67885/87‑1 (MQ=255) aaaaaaa‑CTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc < 2:1138510/89‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 2:397464/1‑90 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 2:612738/1‑90 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 1:1174374/1‑90 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 1:731075/1‑90 (MQ=255) cccccacaaAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga < 1:1003850/87‑1 (MQ=255) | AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAA‑CTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058869‑4059038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |