Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 51.9% | 25.1 / 69.5 | 52 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (19/6); new base A (0/27); total (19/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.32e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058870‑4059038 | aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa < 2:309246/90‑1 (MQ=255) aTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa > 2:213183/1‑90 (MQ=255) aTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa < 2:1121750/90‑1 (MQ=255) ttGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg > 1:1387315/1‑90 (MQ=255) gCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGAt > 1:934949/1‑90 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc > 1:1239825/1‑90 (MQ=255) ccAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCa > 1:53493/1‑90 (MQ=255) aaaTTTTTTTCCGCGAGTCGAGAAAAAATCCGGGGAGCATATCGGAAAACGCGGAAAAAAACTTCCCCACCAGAAAGATAATCCCGCCAt < 1:1249902/89‑1 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg > 2:395816/1‑90 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg > 1:1667535/1‑90 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 1:1422960/1‑90 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGACATAGTAAGc < 1:946836/90‑1 (MQ=255) ttcccGCGAGTCGAAAAAAAACCGGGTGAACATATCGGAAAACGCGAAAAAAACCCTCCCCCCAAAAAAAATAACCCCGCCATAGTAAGc < 2:500870/87‑1 (MQ=255) ccGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAgg > 1:1621027/1‑90 (MQ=255) cgggagccgagaAAAAACCGGGTGAGCATATGGGAAACCGGGGAAAAAACCTTCCCCCCCAAAAAAATAATCCCGCCATAGTAAGCAGGc < 1:1017951/87‑1 (MQ=255) cgcgAGTAGAGAAAAAACCCGGTGAGCATTTGGGTAAACGCGGAAAAAAACTCCGCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGc < 1:738142/90‑1 (MQ=255) gggagtcgagaAAAAACCCGGTGAGCATATCGGTAAACGCGGAAAAAAACTCCCCCACCAGAAAGATAACCCCGCCATAGTAAGCAGGCa < 1:1159418/88‑1 (MQ=255) gggagtcgagaAAAAACCCGGTGAGCATATCGGAAAACGCGGAAAAAACCTTCGCCACCAAAAAGAAAACCCCGCCATAGTAAGCAGGCa < 1:31710/88‑1 (MQ=255) gggagtcgagaAAAAACCCGGTGAGCATATCGGAAAACGCGAAAAAAAACTCCCCCACCAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 1:1289471/88‑1 (MQ=255) gcgAGTGGAAAAAAAACCGGGAGAGCATATGGAAAAACGGGGAAAAAACCTTCCCCACAAAAAAGATACCCCCGCCATAGTAAGCAGGCa < 2:1744060/90‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:336331/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:113843/1‑90 (MQ=255) gcgAGTCGAGTAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:1017951/1‑90 (MQ=255) gcgAGTCGAGAAAAAACCGGGGGAGCATACCGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 2:1783786/90‑1 (MQ=255) gcgAGTCGAGAAAAAACCCGGTGAGCATATCGGAAAACGCGGAAAAAACCTTCCCCCCAAAAAAAATAATCCCGCCATAGTAAGCAGGCa < 1:494906/90‑1 (MQ=255) gcgAGTAAAGAAAAAATCCGGGAAGCATACCGGAAAACGCGGAAAAAACCTCCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 2:113843/90‑1 (MQ=255) gcgAGCCGAAAAAAAACCCGGTGAGCATACCGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 2:307093/90‑1 (MQ=255) gcgACTGGAGAAAAAACCCGGGGAGAATATCGAAAACCGCGAAAAAAAACTTCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 1:40792/90‑1 (MQ=255) gagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATc > 2:710736/1‑90 (MQ=255) aaaaCCGGGTGAGCATATGGGAAACCGCGAAAAAAAACTTCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt < 1:252240/90‑1 (MQ=255) aaaCCCGGTGAGCATATCGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:597805/90‑1 (MQ=255) aacccGGTGAGAATATCGGAAAACGCGAAAAAAACCTTCCCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 1:767751/87‑1 (MQ=255) cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat > 1:1621183/1‑90 (MQ=255) gAGCATATCGGTAAACGCGGTAAAAAACTTCGACACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATaaaa > 1:200801/1‑90 (MQ=255) aatatCGGAAACCGCGGAAAAAACCTTCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 1:651697/89‑1 (MQ=255) tCGGAAACCGCGAAAAAAACCTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 1:1270413/90‑1 (MQ=255) tCGGAAAACGCGAAAAAAAACTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 1:1794212/90‑1 (MQ=255) gcGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:1095048/89‑1 (MQ=255) cGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga > 2:1499948/1‑90 (MQ=255) cGGAAACCGGGGAAAAAACCTCCCCCCCAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 1:832569/90‑1 (MQ=255) gaaaaCGGGGAAAAAACCTTCCCCACAAAAAAAAAACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:1401832/88‑1 (MQ=255) gaaaaCGCGGAAAAAAACTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:1265061/88‑1 (MQ=255) gaaaaCGCGAAAAAAAACTCCCCCCCCAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:573877/88‑1 (MQ=255) gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 1:306453/1‑90 (MQ=255) aaCGCGAAAAAAACCTCCGCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 2:870095/90‑1 (MQ=255) aaCGCGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 2:1373610/90‑1 (MQ=255) gaaaaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 1:1053419/87‑1 (MQ=255) aaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 2:1387315/90‑1 (MQ=255) ttCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt > 2:909934/1‑65 (MQ=255) ttCCCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 1:909934/65‑1 (MQ=255) tcccccacaaAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag < 2:336331/86‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 1:1655960/1‑90 (MQ=255) | AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058870‑4059038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |