Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 56.7% | 7.9 / 33.7 | 30 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (7/6); new base A (0/17); total (7/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.43e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA > NZ_CP009273/4058869‑4059044 | aaaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg < 2:669852/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa < 2:209784/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa < 1:828953/90‑1 (MQ=255) aTTAGCATTTAGGCCCAATCTTTTTCCGCGAGTCGAGTAAAAATCCGGCGAGCATCTCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa < 1:509675/90‑1 (MQ=255) tGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGa > 1:307329/1‑90 (MQ=255) cAATTTTTTTCCGCGAGTCGAGTAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt < 2:321820/90‑1 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 1:599099/1‑90 (MQ=255) tttCCGCGAGTCGAGAAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc < 2:100391/90‑1 (MQ=255) ttcccGGGAGTAGAAAAAAAATCGGGTGAGAATATGGGAAACCGCGGAAAAAACCTCCCCCCCAAAAAAGATAACCCCGCCATAGTAAGc < 1:755149/87‑1 (MQ=255) cgggagtcgagaAAAAACCCGGTGAGCATATGGGTAACCGCGGAAAAAAACTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGc < 2:782658/87‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:54531/1‑90 (MQ=255) agtcgaaaAAAAACCGGGTGAGAATATGGGAAAACGGGAAAAAAACCTCCGCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATcc < 1:215731/90‑1 (MQ=255) aaccccGGGGAGAATATCGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:226565/86‑1 (MQ=255) aaaCCCGGTGAGCATATGGGTAACCGCGGAAAAAACCTTCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 2:843361/90‑1 (MQ=255) aacccGGTGAGCATAGCGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 1:311589/87‑1 (MQ=255) aaTCCGGTAAACATATTGATAAACGCGGAAAAAACCTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 2:168986/90‑1 (MQ=255) taccGGAAACCGGGGAAAAAACCTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:319605/87‑1 (MQ=255) tCGGAAAACGGGGAAAAAACCTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:722964/90‑1 (MQ=255) tCGGAAAACGCGATAAAAAACTTCGCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:599099/90‑1 (MQ=255) gaaaaCGGGGAAAAAACCTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:753362/88‑1 (MQ=255) aaaaaCGCGAAAAAAAACTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:681550/88‑1 (MQ=255) accgcgGAAAAAACCTCCCCCCCCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 2:173819/88‑1 (MQ=255) ggTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt > 2:669158/1‑90 (MQ=255) gaaaaaaaaCTTCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 2:202354/87‑1 (MQ=255) aaaaacCTCCGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 1:628366/90‑1 (MQ=255) aaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 1:23125/90‑1 (MQ=255) cTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa > 1:627252/1‑90 (MQ=255) cTTCGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa < 1:336737/90‑1 (MQ=255) tcccccacaaAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag < 1:355921/86‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 2:817225/1‑90 (MQ=255) ccaGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTa > 2:368956/1‑90 (MQ=255) | AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA > NZ_CP009273/4058869‑4059044 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |