Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,171,107 | 0 | C | A | 51.0% | 24.0 / 68.5 | 49 | E55* (GAG→TAG) | nudF | ADP‑ribose diphosphatase |
Reads supporting (aligned to +/- strand): ref base C (14/10); new base A (0/25); total (14/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.90e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCG > NZ_CP009273/3171021‑3171192 | gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaa > 1:163406/1‑90 (MQ=255) gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaa > 1:1205403/1‑90 (MQ=255) cgacAAACCGAATCTGCTCAATCAGCACACCTTCACCACGCCCTGGGTCAAGGGGTACAAAGACTGCGGCGTGCCCGCGCTCAAAAAttt < 1:919781/87‑1 (MQ=255) cgGCAATCCGAATCTGCTCAATCACCACACCTTCACCACGCACGGGGTCAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAAttt < 1:1541625/90‑1 (MQ=255) gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑cc > 2:87546/1‑90 (MQ=255) cGACTCTTCTCAACCACCACAACTTCACCACCCACGGGGTAAAAGGGTAGCAAGGCTGGGGCGTGACCGCCCTCAAAAATTT‑CCcgccgc < 2:916259/90‑1 (MQ=255) cGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAGGGGTAGAAAGACTGGGGCGTGACCGCGCTAAAAAATTT‑CCcgccgc < 1:1477708/90‑1 (MQ=255) cGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCcgccgc < 2:163406/90‑1 (MQ=255) cGAATCTGCTCAACCACCACAACTCCACCCCCCCGGGGGTCAAGGGGTAGAAAGACGGGGGCGGGCCCGCCCTAAAAAATTT‑CCcgccgc < 1:1311445/90‑1 (MQ=255) cGAATCTGCCCAACCAGCACAGCTTCACCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCcgccgc < 2:354031/90‑1 (MQ=255) aaTCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCAc > 1:1339031/1‑90 (MQ=255) aTCTGCTCAATCAGCACAACTTCATCACGCACGGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTT‑CCCGCCGCAcc < 2:372797/90‑1 (MQ=255) cTCAATCAGCACAACTTCATCCCGCCCTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTT‑CCCGCCGCACCTCATg < 2:498245/90‑1 (MQ=255) cAACCACCACACCTCCACCACGCACTGGGTCAAAGGGTAGCAAGCCGGGGGGGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGAc < 1:1362236/90‑1 (MQ=255) aaaTCAGCACAACTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGAc < 2:952838/89‑1 (MQ=255) aTCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTc > 2:543227/1‑90 (MQ=255) cacaACTTCACCCCGCACTGGGTAAAAGGGTAGAAAGACGGCGGCGTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTg < 2:1339031/90‑1 (MQ=255) aaCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGGGTGACGGTGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGccc < 2:578457/90‑1 (MQ=255) aCTTCACCCCCCACGGGGTCAAGGGGTACCAAGCCGGCGGCGTGCCCGCGCAAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCg < 1:589610/90‑1 (MQ=255) ctccatccCGCACTGGGTCAAGGGGTAGCAAGCCGGCGGCGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGt < 1:1186335/87‑1 (MQ=255) cTTCATCCCGCACTGGGTAAAGGGGAAGCAAGACGGCGGCGTGCCCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGt < 1:1517501/90‑1 (MQ=255) cTTCACCCCCCACGGGGTAAAAGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGt < 2:350662/90‑1 (MQ=255) cTTCACCACCCCCGGGGCAAAAGGGAACCAAGACTGCGGCGTGCCCGCGCAAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGt < 1:1217153/90‑1 (MQ=255) tcatcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTg > 1:1516278/1‑90 (MQ=255) atcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGaa > 2:1103277/1‑90 (MQ=255) atcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGaa > 1:1040083/1‑90 (MQ=255) gCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc > 1:1326607/1‑90 (MQ=255) ccTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:156782/89‑1 (MQ=255) acggggTAAAAGGGTAGCAAGACTGGGGCGTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:249513/87‑1 (MQ=255) ggggTCAAAGGCTACCAAGCTTGCGGCGTCCCCGCGCAAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg < 1:1119180/89‑1 (MQ=255) ccAAAGGGAACAAAGCCTGCGGCGTGCCCGCCCAAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa < 1:460140/89‑1 (MQ=255) aaaaGGGTAGCAAGACTGGGGCGTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa < 2:111666/89‑1 (MQ=255) gacaggTAGCAAGCCTGCGGCGTGCCCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:1047506/86‑1 (MQ=255) aaggggTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 1:1086787/87‑1 (MQ=255) aaaGGGTAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 1:1124507/90‑1 (MQ=255) aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:1040083/90‑1 (MQ=255) aggggTACAAAGACTGGGGCGTCCCCGCGCAAAAAATTTC‑CCCGCCGCTCCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt < 1:1499129/88‑1 (MQ=255) aaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt > 2:1494136/1‑90 (MQ=255) aaGGGTAGCAAGACGGCGGCGTGCCCGCGCAAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt < 1:1181420/90‑1 (MQ=255) aaGGGTAACAAGGCTGGGGGGGGCCCGCCCAAAAAAATTTCCCCCCCCCCCCTAATGCCCAATTTGCCCGTTGAATAGACGATGACGaaa < 1:750830/90‑1 (MQ=255) gTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat > 2:695627/1‑90 (MQ=255) aGCAAGACTGCGGCGTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa < 1:1103277/90‑1 (MQ=255) aCTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGt < 2:146181/58‑1 (MQ=255) aCTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGt > 1:146181/1‑58 (MQ=255) tGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa > 1:440249/1‑90 (MQ=255) cgcgCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGcc < 2:1094184/42‑1 (MQ=255) cgcgCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGcc > 1:1094184/1‑42 (MQ=255) cgcgCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGc < 2:302967/90‑1 (MQ=255) cgcTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg < 1:759354/90‑1 (MQ=255) | GTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCG > NZ_CP009273/3171021‑3171192 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |