Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,224,025 | 0 | C | A | 39.0% | 96.2 / 10.2 | 41 | intergenic (+118/+254) | ycgI/minE | pseudogene/cell division topological specificity factor |
Reads supporting (aligned to +/- strand): ref base C (2/23); new base A (16/0); total (18/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAT > NC_000913/1223880‑1224175 | ggTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGa < 1:283894/151‑1 (MQ=255) gTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaa < 2:645937/151‑1 (MQ=255) ggTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCACCCTATAGGGGAaaaaaaaaaaa > 1:566351/1‑141 (MQ=255) gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAt < 1:244760/151‑1 (MQ=255) aaCTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCATGGGGGAAAAAATAaaaatgaaaatg > 2:226204/1‑150 (MQ=255) ggTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCgg < 2:627033/149‑1 (MQ=255) cGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGtt < 2:531709/150‑1 (MQ=255) tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatc < 1:134680/150‑1 (MQ=255) gACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatca < 2:675290/148‑1 (MQ=255) aTCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCGACCCATGGGGGAAAAAAGAAAAATAAAAATGCCGGGGTTTCATCcccc > 2:542376/1‑150 (MQ=255) aaCGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCAATCTATGGGGGAAAAAAAAATAATGAAAATTGGGGGGTATTATCTccc > 1:11787/1‑147 (MQ=255) aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATACTCTTGGGGAACAATGACTAATGAAAATGAGGGGGTATCATCTCCCAGtat > 1:12244/1‑148 (MQ=255) aaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATa < 2:384646/151‑1 (MQ=255) aaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATa < 2:187386/151‑1 (MQ=255) tGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATCCTATGGGGGAAAAAAAAATAAAGAAATTGGGGGGGTGTTATCTCCCAATAAATACATa > 1:255328/1‑148 (MQ=255) tCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAAt < 1:672289/150‑1 (MQ=255) cGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATTCCCTTGGGGAAAAAAAAATAAATAAAATGGGGGGGTATTATCTCCCCGGATAATCATACTAAAAATa > 1:195020/1‑150 (MQ=255) cGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATaa < 2:85015/151‑1 (MQ=255) aGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCATTCTTTTGGAGGAAAAAGAAAAAATAAAATGGGGGGGTATTATCTCCCCGTATATTCAaaataaaaataa > 1:10579/1‑145 (MQ=255) aGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTa < 1:59888/150‑1 (MQ=255) gTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATTCTATGGGAGAAAAAAGAAAAAAGAAATTGCGGGGGTATCATCTCCCAGTATATCCATACTAACAATAAAGTTAtt > 1:192111/1‑151 (MQ=255) tttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATGGGGGAAAAAAGAATAAAGAAATTTCGGGGGTATCATCTCCCCGTATATCCCTACTAAAAATAAGGTTAttt > 2:421168/1‑150 (MQ=255) ttCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTa < 2:205466/150‑1 (MQ=255) tCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACACATCCTATTGGAGAAAAAAGAATAATGAAAATTCGGGGGTATCATTTCCCAGTATATCCATACTAACAAAAAAGGTAttt > 1:354052/1‑148 (MQ=255) cGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACt < 1:53724/150‑1 (MQ=255) cGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCCATGGGAGAAAAAAGAATAATGAAAATGCGGGGGTATCATCTCCCAGTATATCCATAATAACAATAAGGGTATTTaat > 2:415150/1‑148 (MQ=255) cGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCCATGGGGGAAAAATGAATAATGAAAATGCGGGGTTATCATTTCCCCGTATATTCATACTAACAAAAAAGGTAttt > 1:327350/1‑147 (MQ=255) gTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTc < 2:204219/150‑1 (MQ=255) ttACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTTTGGGGGAAAAAAGAATAATGAAAATTCGGGGGTTTCATCTCCCCGTATATTCATACTAAAAATAAAGTTATTTACTcca > 2:85663/1‑148 (MQ=255) tACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAAc < 2:652731/150‑1 (MQ=255) cAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAAAACCATCCCATGGGGGAAAAAAGAATAAATAAAATGGGGGGTTATCATCTCCCAGAATATTACTACTAACACTAAAGGTAtttttt > 1:453998/1‑141 (MQ=255) gAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATaa < 1:241502/148‑1 (MQ=255) gggAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACAtttt < 2:271656/151‑1 (MQ=255) aaaaTTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGGAACCTCATCCTATGGGAGAAAAAAGAAAAATGAAAATGCGGGGGTATCATTTCCCCGTATAACCACACTAACAAAAAGGGTATTTTCTCAACCACGCATAAACAttttgtt > 1:186789/1‑151 (MQ=255) ggggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTtgtg < 1:495475/148‑1 (MQ=255) gggCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTg < 1:68558/151‑1 (MQ=255) cAAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCTATGGGGGGAAAAAGAATAATTAAAATTCGGGGTTATCATTTCCCAATATATCCACACTTAAAATAAGGTTATTTTTTCAAACAGGCACTAACAATTTTTTTTTTGCGTggg > 2:209904/1‑150 (MQ=255) aaaaaGCCCGCCAGTTACGGCGGGGAACCTCATACTATGGGGGAAAAAAGAAAAATTAAAATGCGGGGGTATTCTCTCCCCGGATAATCATAATAACAAAAAAGTTTTTTAATCCAACACGCATAAACAATTTTTTTTTTGGGTgg > 1:119553/1‑146 (MQ=255) ggggAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAggggg > 2:44412/2‑151 (MQ=255) cTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGc < 2:374708/150‑1 (MQ=255) gagaACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgagg < 1:666640/151‑1 (MQ=255) gaACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATgaggag < 2:192149/151‑1 (MQ=255) cAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAt < 1:139526/151‑1 (MQ=255) | GGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTAT > NC_000913/1223880‑1224175 |
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: ‑ |