Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,738,402 | 0 | A | C | 52.9% | 21.3 / 32.2 | 35 | D80A (GAC→GCC) | yiaN | 2,3‑diketo‑L‑gulonate transporter large permease YiaN |
Reads supporting (aligned to +/- strand): ref base A (9/7); new base C (0/18); total (10/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.18e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTC > NZ_CP009273/3738314‑3738486 | gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAAAGCGGGGGGGCTGTTCAAAAGCAATGTTGAc > 1:747652/1‑90 (MQ=255) gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTg < 2:746594/90‑1 (MQ=255) gATAGCTTCCCCCCGCGGGCTTTCCCTTTTTTTGTGCGGGCGGGGAAAATAATGAATGCGGGGGGCTTCTAAAAACGCATTGTTCCCCTg < 1:749341/90‑1 (MQ=255) cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAAAGCATTGTTGACCTGCCGAt > 1:633492/1‑90 (MQ=255) ttCCCCCTGCTGGCGATTCCGTTTTTTGTGCGGGGGGGTAAAATCATGAATGGGGGGGGCTTGTAAAAACGCATTGTTGACCTGCCGATg < 1:1354222/90‑1 (MQ=255) tctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGCGGGGGGGGTTTCCAAAAGCATTGTTGGCCTGCCGATGa > 1:1343704/1‑90 (MQ=255) cccTGCTGGCGTTTCCGTTCTTTGTGGGGGGGGGTGAAATCATGAGTGCGGGGGGTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAAc < 1:1339849/90‑1 (MQ=255) cccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAAAGCATTGTTGACCTGCCGATGAAAc > 2:817659/1‑90 (MQ=255) ccTGCTGGCGTTTCCGTTTTTTGTGTGGGGGGGTAAAATCATAAAGGGGGGGGGTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAACt < 1:1162626/90‑1 (MQ=255) ctgGCGATTCCTTTTTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggtg < 2:946163/90‑1 (MQ=255) cgggCGATTCCGTTTTTTGTGCGGGGGGGTAAAATCATGAAGGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACtggtg < 2:1453481/88‑1 (MQ=255) tttCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCAt < 2:913099/89‑1 (MQ=255) aTTCGGTTCTTTGTGCGGGGGGGTAAAATCATGAAGGGGGGCGGCTTGTCAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCAt < 1:1572979/90‑1 (MQ=255) aTTCCTTTCTTTGTGCTGGGGGGTGAAATCATGAATGGGGGCGGTTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCAt < 2:1218109/90‑1 (MQ=255) ttttGTGCTGGCGGGTAAAATCATGAAGGGGGGCGGTTTGCAAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAgg < 1:1570430/89‑1 (MQ=255) cTTTGTGCTGGGGGGGGAAGTTATGAAGGCGGGCGGCTTGCAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAgg < 2:1581767/90‑1 (MQ=255) tttGTGCTGGCGGGTAAAATCATGAAGGGGGGGGGCTTGTCAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGc < 1:1018703/90‑1 (MQ=255) tgtgggCGGGTAAAATCATAAATGCGGGGGGCTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc < 1:1225186/86‑1 (MQ=255) tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc < 2:140947/90‑1 (MQ=255) gCTGGGGGGTGAAATCATGAAGGGGGGGGGTTTGTAAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggct < 1:1114497/90‑1 (MQ=255) gCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggct < 2:256058/90‑1 (MQ=255) aTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGc > 2:1557935/1‑90 (MQ=255) tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGc > 2:855713/1‑90 (MQ=255) ggcgggcggCTTGAAAAACCGATTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGc < 2:1461241/89‑1 (MQ=255) ggggggggCTTGCCAAACCGTTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCa < 1:1163475/88‑1 (MQ=255) gcgggcggCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCa < 2:997321/90‑1 (MQ=255) gggcggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc > 2:1565850/1‑90 (MQ=255) cggcTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTAt < 1:855713/90‑1 (MQ=255) accGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGa < 2:717086/88‑1 (MQ=255) accGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGa < 2:1172726/88‑1 (MQ=255) gCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGcc < 2:936635/90‑1 (MQ=255) cATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCt > 1:411021/1‑90 (MQ=255) aTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCtt > 2:217218/1‑90 (MQ=255) ttgttgACCTGCCGATGAAACTGGTGGGGCATAAACCA‑GCGGGCTGGGCTACCTCGGCGTGCTGGGAGCTATGATTATGGGGAGCCcttt > 2:571759/1‑89 (MQ=255) tgttgCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTc < 1:1136294/90‑1 (MQ=255) | GCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTC > NZ_CP009273/3738314‑3738486 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |