Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,346,703 | 0 | C | T | 44.1% | 36.5 / 13.8 | 35 | L208L (CTG→CTA) | fumB | anaerobic class I fumarate hydratase (fumarase B) |
Reads supporting (aligned to +/- strand): ref base C (3/16); new base T (15/0); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.48e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTG > NC_000913/4346563‑4346850 | tCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGA‑GGAAGGTTTTTAGTTTGCCgg > 2:543681/1‑151 (MQ=255) aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGGGGGCGCAGTACCGAGGGGACGCATTTTTTTGACGAGGGAGGTTTTTAAGTTGGCgggg > 1:448285/1‑151 (MQ=255) aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGa < 1:525649/151‑1 (MQ=255) cGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGGAGTACCGCGGGGACGCATTTTTTCGACGA‑GGGAGGTTTTTCGTTTGGCggggg > 1:166737/1‑148 (MQ=255) aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCTACGA‑GGAAGGTTTTTAGTTTGCCggggg > 2:54609/1‑145 (MQ=255) aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCCAGGGGAAGCATTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCggggg > 2:539788/1‑144 (MQ=255) aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGGGGGGAAGCATTTTTTTGGCGA‑GGGAGGTTTTTAGTTTTCCGGGGGGCAGCAggg > 1:473450/1‑149 (MQ=255) gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGCAGGGCTGCAGGGACGAGGGGGCGCCTTTTTTCGAGGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCgggggctt > 2:322470/1‑146 (MQ=255) cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGGGGGGACGGCTTTTTTTGGCGA‑GGGAGGTTTTTAGTTTTCCGGGGGTcagacggg > 1:643120/1‑140 (MQ=255) cTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTATCGAGGGGACGCCATTTTTCGACGA‑GGGAGGTTTTCAGTTTGCCGGGGGTCAGCAGGGGttt > 2:684663/1‑149 (MQ=255) tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCCTTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCGGGGGGCAGCCGGGCTTTGGtt > 2:91077/1‑151 (MQ=255) cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGA‑GGGAGGTTTTTAGTTTTCCGGGAGTCAGCAGGGGTTTTGttt > 2:443060/1‑150 (MQ=255) aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt < 1:441710/151‑1 (MQ=255) aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt < 2:553710/151‑1 (MQ=255) aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg < 2:88699/151‑1 (MQ=255) tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg < 1:535738/149‑1 (MQ=255) aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAGTACCGAGGGGAAGCATTTTTTCCACGA‑GGAAGGTTTTTCGTTTGCCGGGGGGCAGCAGGGCCTTGGTTTTCTGGTag > 2:193176/1‑150 (MQ=255) aCGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTGCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTaga < 1:366511/151‑1 (MQ=255) ggatttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGTTGCAGCACCGAGGGGGGGGACTTTTTTGATGA‑GGGAGTTTTTTAAATTGCCGGGGGTGCGCCGGGGTTTTGTTTTCTg > 2:661702/4‑144 (MQ=255) ttttCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCGGCAGTAAAGAGGGTACGCATTTTTTTGGCGA‑GGGAGGTTTTTCGTTTGCCGGGGGGCCGCAGGGCTTTGGTTTCCTGGgagaga > 1:284947/1‑149 (MQ=255) tttCATGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCATGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTCGTTCCCTGTTAGAGATa < 1:88745/150‑1 (MQ=255) aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGGCTTTTCTCGACGA‑GGGAGATTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTTGTTTTCTGGgagagaaatggtt > 2:70845/1‑144 (MQ=255) ttggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCCGGGTGCAGTACCGAGGGGACGCATTTTCTCGCCGA‑GGGAGATTTTTAGGTTGCCGGGAGTCCGCAGGGGTTTGGTTTCCTTGGAGAGATAAGttt > 1:411366/1‑148 (MQ=255) ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCaga < 2:200344/151‑1 (MQ=255) gCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGCGGGTGCAGTACCGAGGGTACGCACTTTCTCGACGA‑GGAAGGTTTTCAGTTTGCCGGGGGTGAGCAGGGGGTTTGGTTCCTGGGAGAGATACGTTTTTTTGGCagag > 1:112320/1‑150 (MQ=255) gACGTACCGCCAAACACAAACGCGATATGGTACGGCGGGCAGGGCGCCGAAACCGGGGGACGCATTTTTTCGAAGA‑GGGAGTTTTTTTGGTTGTCGGGGGGTAGCAGGGGTTTGGGTTCCTGGTAGGGATACGttttt > 1:564209/1‑137 (MQ=255) gCCAATCACAAACGCGATATGGTACGGCGGGGAGGCTGCAGGACCGCGGGTACCCATTTTTTCGACGA‑GGAAGTTTTTTAGGTTGCCGGGGGGCGGCAGGGCTTTTGTTTCCTGGTAGAGATACGGTTTGTTGGGAGAGCCGCCGCCctt > 2:230753/1‑150 (MQ=255) caAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacgca < 1:573055/151‑1 (MQ=255) atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt < 1:459174/151‑1 (MQ=255) atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAtt < 1:304212/151‑1 (MQ=255) cggcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCa < 1:641176/151‑1 (MQ=255) aGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg < 1:579006/151‑1 (MQ=255) gTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGt < 2:442706/151‑1 (MQ=255) ggAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGt < 1:598788/151‑1 (MQ=255) ttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGgtg < 2:473469/150‑1 (MQ=255) | TCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGA‑GGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTG > NC_000913/4346563‑4346850 |
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: ‑ |