Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 993,913 | 0 | G | C | 50.0% | 47.5 / 13.2 | 34 | G44G (GGC→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base G (16/1); new base C (0/17); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG > NC_000913/993768‑994050 | gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCccac > 2:489065/1‑151 (MQ=255) aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCccacca < 2:269845/151‑1 (MQ=255) cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCcac > 2:341936/1‑150 (MQ=255) cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCcac > 1:823746/1‑150 (MQ=255) catTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCgg > 1:649652/1‑150 (MQ=255) ttggcttgtttccgtGAATTT‑CACCCGGGGGTGTGGGCCCCGCTAACCCCTCG‑CTGGGG‑TTGGGCTTTTCCAGCCCGCCCAGCGGGCGCAGCAGGTTTTTTTGC‑CCACCCCCG‑CTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa < 2:17867/135‑1 (MQ=255) cgcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat > 2:93586/1‑151 (MQ=255) gcgTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg > 2:167522/1‑151 (MQ=255) tctcctttattct‑CAGCCGGGGGGGGGGTCCCCGCTAACACCTCGCTTGCGG‑TT‑GGCTTTCCAAGCCCGGCCAGCGGGCGCAGCAGGTTACTTTGC‑CCCCCCCCG‑CGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTaaa < 2:243849/138‑1 (MQ=255) tCCTTAATTTCCGGCCGGGGGTGTGGCCCCCTCAAACCCATCGCCTGGGG‑TG‑GGCTTTCCCCCCCCCGCCAGCGGGCGCAGCGGGTTATTTTCC‑CCACCACCG‑CGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGt < 1:619908/145‑1 (MQ=255) tCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACACACTGACCTGCCGGAATATGTAAATcc > 2:807527/1‑150 (MQ=255) tCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc > 2:423426/1‑150 (MQ=255) ccTAAATTC‑CGGCCGGGGGTGGGGGCCCCGCAAACACACCGCCGGCGT‑TG‑GGCTTCCCCAGCCCCCCCAGCAGGCGCAGCAGGTTTTTTTCC‑CCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg < 1:407791/151‑1 (MQ=255) cTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt > 1:149548/1‑151 (MQ=255) tt‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc > 2:267540/1‑151 (MQ=255) ccaggGGGGGGGTGCCCGCAAACCCATC‑CCTGGGGTTT‑GGCTTTCCCAGCCCTCCCAGCGGGCGCAGGGGGTTTTTTTTC‑CCCCCCCCG‑CTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa < 2:738013/147‑1 (MQ=255) cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg > 1:647798/1‑151 (MQ=255) gcGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat > 1:157670/1‑151 (MQ=255) gtggtcCCCGCAAACACCCCCCCGGGGG‑TG‑GGCGTTCCCACCCCGGCCGGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt < 2:799230/149‑1 (MQ=255) ggcccccGCAAACCCATCGCCTGGGG‑TT‑GGCGTTCCCGGCCCCCCCAGCGGGCGCGGGGGGGTTCTTTTC‑CCACCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc < 1:691398/147‑1 (MQ=255) gccccccCTAACCCACCGCCTGGGT‑TT‑GGTGTTCCAAGCCCGCCCAGCAGGCGCAGCAGGTTACTTTTG‑CCCCCACCG‑CGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt < 2:457366/144‑1 (MQ=255) ccccatAACCCACCGCCTGCGT‑TG‑GGCTTTCCCACCCCGGCCAGCGGGCGCAGCGGGGTTTTTTCC‑CCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg < 2:94178/145‑1 (MQ=255) tAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt > 2:788438/1‑151 (MQ=255) cacaCCGCCTGGGG‑TG‑GGTTTTCCCAGCCCGCCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCG‑CTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt < 2:242563/150‑1 (MQ=255) caccgccggggttgggCTTTCCCAGCCCTCCCAGCAGGCGCAGCGGGGTATTTTTG‑CCCCCCCCG‑CTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg < 2:14495/135‑1 (MQ=255) ttcccagcccgcccAGCAGGCGCAGCA‑GGTTTTTTTC‑CCCCCCCCGCCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt < 1:191056/118‑1 (MQ=255) ttcccaccccTCCCAGCGGGCGCAGCGGGGTCCTTTTC‑CCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca < 2:477561/139‑1 (MQ=255) cccagcccccccAGCGGGCGCAGCGGGTTTTTTTCC‑CCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg < 1:560457/141‑1 (MQ=255) ccggCCAGGAGGCGCGGCAGGGT‑CTTTTGCCCCCCCCCG‑CGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca < 1:39274/139‑1 (MQ=255) agGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc > 2:32591/1‑149 (MQ=255) agGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc > 1:238124/1‑149 (MQ=255) tccccccccccg‑ccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt < 1:715115/137‑1 (MQ=255) cccccaccG‑GTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg < 1:238410/147‑1 (MQ=255) caccacCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg > 2:404734/1‑151 (MQ=255) | GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGG‑TT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCG‑CTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG > NC_000913/993768‑994050 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |