Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,337,723 | 0 | A | C | 45.8% | 73.5 / 11.0 | 49 | G566G (GGT→GGG) | gyrA | DNA gyrase (type II topoisomerase), subunit A |
Reads supporting (aligned to +/- strand): ref base A (23/3); new base C (0/22); total (23/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.48e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC > NC_000913/2337576‑2337871 | gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt > 2:296228/1‑151 (MQ=255) gcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCttt > 2:344792/1‑151 (MQ=255) gcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTcc < 1:214537/151‑1 (MQ=255) ccACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGc < 1:587603/151‑1 (MQ=255) aCGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc > 2:1032969/1‑150 (MQ=255) tGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCacgacg > 1:941451/1‑151 (MQ=255) tGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAtt > 1:741244/1‑151 (MQ=255) aacttttcccgaaaTAGCCGCGCCCAGGGCTGGAAAAGCACAGAATAGGGTCGTGAGTGTTCCCCACCAGCATTCGGTTTATAAATCCTTCTTCTTAAATTCGTGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTc < 1:581367/136‑1 (MQ=255) tttCATCAAATAGACGCGCCCGCGGCGGGAAAAGCAAAAAATATGGTCGTGATTTTTCGCCCCCAGCAGTGGGTCAATAAAGTCTTCTTTTTTAATCCGGGCGGAAAATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa < 2:849125/149‑1 (MQ=255) aTCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCTCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt > 1:718636/1‑151 (MQ=255) aaTAGCCGCGACCCCGGTGGGAAAAGCAAAAAATATGTTCGGGAGTTTCCCCCCCCAGCAGTCGGTCAAAAAATTTTTCTTTTTTAATACGGGCGGCAATTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCg < 1:147018/144‑1 (MQ=255) gccgggcccggggggggAGAACCTCAGAATTTGGTGGTATGTGTGCCCCGCCAACACTGGTTCGATAAGGTCTTCTTATTTAATGTGGGCGTGACTTTGCCCTTGTCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa < 1:328983/136‑1 (MQ=255) cgcccAGGGCGGGAGAAGCACAGAATATGGTCGTGAGTTTCCCCCACCAGCAGTCGGTCAAAAAAGTCTTCTTTTTTAATGCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc < 1:1107931/146‑1 (MQ=255) ccgcggcggGAGAAGCAAAAAATATGGTCGTGAGTGTCCCCCACCACCAGCCGGTCGATAAATTCTTTTTCTTTAATCCGTGCGGAAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg < 2:741223/144‑1 (MQ=255) ccccggcggGAGAAGCACAGAATATGGTCGTGAGTTTGCCCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc < 1:253198/143‑1 (MQ=255) ccACGGCGGAAGAAGCAAAAAATTGGTTCGTAGGTGTCCCCCACCAGCGGCCGGTGAAAAAATTCTTTTTTTTTAATGCGTGCGGCAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt < 2:783627/148‑1 (MQ=255) cACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc > 1:714542/1‑151 (MQ=255) ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg > 1:419523/1‑151 (MQ=255) ggAGAAGAACAAATTATGTTGGTGAGTTTCCCCCCCCAGCGGTCGGTCAAAAAGTTCTTTTTTTTTAACACGTGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTgg < 2:263900/148‑1 (MQ=255) aaGCACAGAATATGGTCGGGAGTGTTCGCCCCCAGCAGCCGGTCAATAAAGTTTTCTTTTTTAATACGGGGGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagaga < 1:286401/151‑1 (MQ=255) gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCAc > 2:423470/1‑151 (MQ=255) aaTATGGTCGTGAGTGTCCGCCCCCAGCATTCGGTCGATAAAGTCTTCTTCTTTATTACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGt < 2:887008/147‑1 (MQ=255) tatGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAc > 2:392681/1‑151 (MQ=255) tGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt > 1:844208/1‑151 (MQ=255) tGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt > 1:421248/1‑151 (MQ=255) gAGTGTTCGCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc < 1:437391/148‑1 (MQ=255) gtgtTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt > 2:528181/1‑151 (MQ=255) gtgtCCCCCCCCAGCAGTCGGTCGATAAAGTTTTCTTTTTTAATACGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt < 2:979227/148‑1 (MQ=255) gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg > 2:573694/1‑150 (MQ=255) cccccaagcggTCGGTCAAAAAATTCTTCTTTTTTAATACGTGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGa < 2:468948/145‑1 (MQ=255) ccaccaGCAGTCGGCGGAAAAATTCTTCTTTTTTATTCCGGGGGGAAATTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTcc < 2:984017/143‑1 (MQ=255) ccaGCAGTCGGTCAATAAAGTCTTCTTTTTTAATACGGGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCa < 1:296238/150‑1 (MQ=255) cagTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCt > 1:146665/1‑151 (MQ=255) cggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc > 2:225051/1‑150 (MQ=255) ggtcgATAAATTCTTCTTTTTTAATGCGTGCGGAAGTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAg < 2:300744/151‑1 (MQ=255) ggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc < 2:248010/149‑1 (MQ=255) cgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt > 2:380534/1‑151 (MQ=255) tAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt > 2:647439/1‑151 (MQ=255) aattcttcttctTTAATGCGGGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt < 2:20118/148‑1 (MQ=255) aGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc > 1:1127262/1‑151 (MQ=255) ttcttTAATACGTGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc < 1:403388/147‑1 (MQ=255) ttcttAAATACGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg < 1:544600/148‑1 (MQ=255) tttaaagtgTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg < 2:4121/137‑1 (MQ=255) aaTACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCg > 1:1077727/1‑151 (MQ=255) tgggTGGGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg < 1:344803/147‑1 (MQ=255) tACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg > 1:855367/1‑150 (MQ=255) tGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAtt > 1:906533/1‑151 (MQ=255) gcgcaaaaTTTTCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAta < 1:1032995/143‑2 (MQ=255) tACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac > 2:427551/1‑150 (MQ=255) | GCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC > NC_000913/2337576‑2337871 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |