Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,198,949 | 0 | A | T | 29.6% | 36.2 / 17.1 | 29 | intergenic (‑250/+27) | ethA/PP_2806 | FAD‑containing monooxygenase EthA/TetR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (9/10); new base T (3/5); total (12/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCAGCAA > NC_002947/3198811‑3199098 | aTTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCTGCGTCGCTGCACGGc < 1:297353/150‑1 (MQ=255) aTTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGc < 1:585056/150‑1 (MQ=255) cAGGATGATTTGCACTTTTTACTTCGCCGATCAACGATCAGACTTGGTCAGTCTGGTAGCTATCGTCTTCGGGACGCTAGACGGGTGGGGCTGCGGTTTCTGGGGCTGTAGCGGGACTGGATTGCCGTGGTGCGTCGCTGCACGGCGatc < 2:885277/150‑3 (MQ=255) aGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATc < 1:555661/150‑1 (MQ=255) ccTGTTTCTTCGCGGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCa < 2:897909/150‑1 (MQ=255) ccTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCa < 2:563389/150‑1 (MQ=255) tGTTACTTCGCCGATCAACGATCAGACTTGGTCAGTCTGGTAGATATCGTCATCGGAACCTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCTGCGTCGCTGCACg < 1:622001/128‑1 (MQ=255) tGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCTGCGTCGCTGCACg < 1:621901/128‑1 (MQ=255) tGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCTGCGTCGCTGCACg > 2:622001/1‑128 (MQ=255) tGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCTGCGTCGCTGCACg > 2:621901/1‑128 (MQ=255) aCTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACt < 2:227614/150‑1 (MQ=255) cTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTg < 1:984837/150‑1 (MQ=255) tCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCGGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGc > 1:298564/1‑150 (MQ=255) gATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGccacca > 1:600139/1‑150 (MQ=255) aTCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgca < 2:914270/103‑1 (MQ=255) aTCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgca > 1:914270/1‑103 (MQ=255) tCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcg < 2:678861/150‑1 (MQ=255) tGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggt > 2:852739/1‑150 (MQ=255) ggAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCa > 2:1018904/1‑150 (MQ=255) gAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg < 2:45537/87‑1 (MQ=255) gAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg > 1:45537/1‑87 (MQ=255) cGAAGGGGGGCGCTGCGGCTTTCAGATCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtg < 1:287830/76‑1 (MQ=255) cGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtg > 2:287830/1‑76 (MQ=255) cTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGCCGCCGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGa < 2:534394/150‑1 (MQ=255) gCGGCTTTCGGCGCGGGCGGCGGCGGGGCTAGCCGTGCTGGGCCGCAGCACGGCGAATCCTGCGGGGCCGGAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGc < 2:315161/150‑1 (MQ=255) gCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAg > 1:871669/1‑150 (MQ=255) gaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCATTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa > 2:468458/1‑150 (MQ=255) gCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTc < 1:878020/143‑1 (MQ=255) gCCGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTc > 2:878020/1‑143 (MQ=255) aCGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCAGCaa < 1:1018904/150‑1 (MQ=255) | ATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCAGCAA > NC_002947/3198811‑3199098 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |