| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,337,723 | 0 | A | C | 56.8% | 21.5 / 43.0 | 37 | G566G (GGT→GGG) | gyrA | DNA gyrase (type II topoisomerase), subunit A |
| Reads supporting (aligned to +/- strand): ref base A (11/5); new base C (0/21); total (11/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.11e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-06 | |||||||||||
CGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC > NC_000913/2337586‑2337871 | cGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCCCGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCa < 1:63674/150‑1 (MQ=255) ttCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTgc > 2:422658/1‑151 (MQ=255) aTAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc > 1:128971/1‑150 (MQ=255) aTAAACTTTAACGAAATACACGCGCCAACGGCGGAAGAGGCAGAGAATAGGGTCGTGGGTATTCGCCGCCAGCTGCCGCCAGAAAAAGTGTTTTTCTTAAAAGCGGGCGGTACATTTCCCTTCCCCGCCACGACGCTGCGCTTCGt < 1:54377/146‑1 (MQ=255) tttCATCGAGTAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGc > 2:232348/1‑151 (MQ=255) cATCGAAAAGACGCGCCCAGGGTGGGAGAAGCACAGAATTGGGTCGGGAGTGTTCGCCACCAGCGGTCGTTCGATAAAGTCTTCTTCTTTATTCCGGGCGGAAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg < 1:259053/150‑1 (MQ=255) tCGAATAGCGGCGACCCCGGCAGGAGAACCACAGAATATGGTCGTGAGTGTTCCCCCCCAGCTGTCGTGCGATAAAGTCTTTTTTTTTAATCCGTGCGGCAGGTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGc < 1:45972/146‑1 (MQ=255) aaTAGCCGCGACCACGGCTGGAGAAGCACAGAATATGTTCGTGATTGTTCGCCACCAGCAGTCGGTCAAAAAAGTCTTCTTTTTTAATACGTGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt < 2:288609/147‑1 (MQ=255) cgACCAGGGGTGGAGAAGAAAAAAATATGGTCGGGATTTTTCGCCACCAGCAGTCGGTCGAAAAAGTCTTTTTTTTTAATCCGGGCGGAATATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg < 2:286351/150‑1 (MQ=255) ccgcggtgggAGAAGCACAAAATTTGGTCGTGATTGTTCCCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATGCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAg < 1:382662/142‑1 (MQ=255) cccgggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCCCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTa < 2:297687/145‑1 (MQ=255) ggCGGGAAAAGCACAAAATTTGGTGGTGAGTGTTCCCCCCCAGCATTCGTTGGATAATGTCTTCTTTTTTAACACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt < 1:329848/150‑1 (MQ=255) agccagaGAATAGGGTCGTGATTGTCCCCCACCAGAAGTGGGTCGATAAACTCTTCTTCTTAAATCCGTGCGGGAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGccc < 2:577038/135‑1 (MQ=255) atagGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTTTTTTTTAATACGTGCGGAAGATTTACTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCAc < 1:431761/149‑1 (MQ=255) gttCGTGAGTGTTCGCCCCCAGCAGTCGTTCGATAAAGTCTTCTTCTTTAATACGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAcc < 1:349966/147‑1 (MQ=255) tCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACAt > 2:247137/1‑151 (MQ=255) aGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTcc > 2:443888/1‑151 (MQ=255) ttttCGCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATCCGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt < 2:341636/148‑1 (MQ=255) cGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg > 2:462564/1‑151 (MQ=255) cccccaGCAGTCGGTCGATAAAGTCTTCTTTTTTAATCCGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGa < 2:236621/147‑1 (MQ=255) ccaGCAGTCGGTCAATAAAGTCTTCTTTTTTAATACGTGCGGCATATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAg < 1:47883/151‑1 (MQ=255) ggtcgATAAGGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCa < 2:293759/150‑1 (MQ=255) cgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt > 2:313451/1‑151 (MQ=255) aaaaaGTCTTCTTCTTAAATCCGTGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg < 1:74699/148‑1 (MQ=255) aTAAAGTCTTCTTCTTTAATCCGGGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTg < 1:567330/150‑1 (MQ=255) aaaTTCTTCTTCTTTAATGCGTGCGGCAAATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg < 2:118053/151‑1 (MQ=255) tcttcttctTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt < 2:192182/150‑1 (MQ=255) ttttttAATGCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTg < 1:497796/143‑1 (MQ=255) tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt > 2:251353/1‑151 (MQ=255) tttAATACGTGCGGCAGTTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGtt < 1:577940/150‑1 (MQ=255) aaTCCGTGCGCCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCCGTTgg < 2:245501/149‑1 (MQ=255) aaTACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc < 2:557684/150‑1 (MQ=255) ataTTTCCTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAttt < 1:361593/144‑1 (MQ=255) aGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt > 2:291151/1‑150 (MQ=255) aTTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAATCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGGGATTTCAGTacg > 1:156796/1‑151 (MQ=255) tttCCCTTGCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTa < 1:569584/148‑1 (MQ=255) ttACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac > 1:229182/1‑151 (MQ=255) | CGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAC > NC_000913/2337586‑2337871 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |