Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,331,202 | 0 | A | C | 52.4% | 20.9 / 26.3 | 21 | G566G (GGT→GGG) | gyrA | DNA topoisomerase (ATP‑hydrolyzing) subunit A |
Reads supporting (aligned to +/- strand): ref base A (8/2); new base C (0/11); total (8/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACA > NZ_CP009273/2331124‑2331286 | gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc > 2:413434/1‑90 (MQ=255) agatatGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGAAAAATTCTTTTTTTTTAATGCGGGCGGCATATTTACCTTTCCCGCCacgac < 1:23868/88‑1 (MQ=255) atatGGTCGTGGGTGTTCGCCCCCAGCAGTGGGTCAAAAAAGTCTTTTTCTTTAATCCGTGCGGCATATTTCCCTTCCCCGCCACGACGc < 1:597115/90‑1 (MQ=255) atatGGTCGTGAGTTTTCCCCACCAGCAGTCGGTCGAAAAAGTCTTTTTTTTTAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGc < 1:625245/90‑1 (MQ=255) atGGTCGTGATTGTCCCCCACCAGCGGCGGGTCAAAAAATTTTTTTTTTTAAATCCGTGCGGCATTTTTCCCTTCCCCGCCACGACGCTg < 2:460316/90‑1 (MQ=255) gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc > 1:75449/1‑90 (MQ=255) gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc > 2:367510/1‑90 (MQ=255) gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc > 2:271428/1‑90 (MQ=255) gAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTcct > 2:23868/1‑88 (MQ=255) ccccccaGCAGTGGGTCAAAAAAGTCTTTTTTTTTAATCCGTGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGaa < 1:367510/86‑1 (MQ=255) ccaccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa > 1:618779/1‑90 (MQ=255) ccccaGCAGTGGGTCAAAAAAGTTTTTTTTTTTAATACGGGGGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAg < 1:679804/88‑1 (MQ=255) tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAc > 2:587860/1‑44 (MQ=255) tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAc < 1:587860/44‑1 (MQ=255) gaaaaaTTCTTTTTTTTTAATGCGTGGGGCATTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc < 2:75449/87‑1 (MQ=255) cttcttcttTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt < 1:181021/90‑1 (MQ=255) ttAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag > 2:122314/1‑90 (MQ=255) aaTCCGGGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag < 1:611491/90‑1 (MQ=255) agttttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac < 1:653262/87‑1 (MQ=255) gttttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcaca < 1:413434/88‑1 (MQ=255) gATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcaca < 1:122314/90‑1 (MQ=255) | GCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACA > NZ_CP009273/2331124‑2331286 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |