Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,753,266 | 0 | A | C | 24.4% | 106.1 / 15.5 | 41 | H491P (CAC→CCC) | recN | DNA repair protein RecN |
Reads supporting (aligned to +/- strand): ref base A (30/1); new base C (0/10); total (30/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.81e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTG > NC_000913/2753212‑2753328 | tCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACaa > 1:2449654/1‑65 (MQ=255) gCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTc > 1:326266/1‑65 (MQ=255) gCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGAGTGTTACCCACCTGCCACAAGTc > 1:1113501/1‑65 (MQ=255) ctgctgCTTCAACGCGGCGAATCAACTCAGGGGTTGTTTTTTCCCCCCCTGCCACAAGTCGCggg < 1:231897/65‑1 (MQ=255) ctgctgCCTCAACTCGGCGAATCAACTCAGGTGGTGTGTTTTACCCCCCTGCCACAAGTCGCggg < 1:1961173/65‑1 (MQ=255) gctgctTCAACTCGGCGAATCACCTCAGGTGATTTGTTTTCCCCCCCTGCCACAAGTCGCGGGAt < 1:1220538/65‑1 (MQ=255) gTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCCCAAGTCGCGGGATGTGGt > 1:1352049/1‑65 (MQ=255) gTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGt > 1:2698130/1‑65 (MQ=255) gTCAACTCGGAGAATCAACTCAGGTGATGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGt < 1:2744682/65‑1 (MQ=255) gTCAACCCGGCAAAGCAACTCAGGTCATGTCTGTTACCCACCTGCCACAAGTCGCGGGATGTGGt < 1:1169755/65‑1 (MQ=255) cagcTCGGCAAAAAAACTCAGGTGGTGTTTGTTACCCGCCTGCCACAAGTCGCGGGATGTGGtca < 1:3011056/62‑1 (MQ=255) cAACGCGGCGAATCAACTCAGGTGTTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtca < 1:1459130/65‑1 (MQ=255) cgcGGCGAATCAACTCGGGTGGTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtcatca < 1:2500900/63‑1 (MQ=255) cgcGGCGAATCAACGCAGGTGAGTTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGtcatca < 1:968250/63‑1 (MQ=255) cGGCGAATCAACTCAGGTGATGTGTTTTCCCCACCTGCCACAAGTCGCGGGATGTGGTCATCAac < 1:546423/65‑1 (MQ=255) gcaaaTCAGCTCAGGTGTTTTTTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAacac < 1:2201408/62‑1 (MQ=255) tCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTAttt > 1:2870846/1‑65 (MQ=255) tCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCa > 1:1329235/1‑65 (MQ=255) tGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCaaa > 1:2504699/1‑65 (MQ=255) gATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAg > 1:2355803/1‑65 (MQ=255) aTGTGTTTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGa < 1:437857/65‑1 (MQ=255) aTGTGTGTTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGa < 1:1900293/65‑1 (MQ=255) tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa > 1:411899/1‑65 (MQ=255) tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa > 1:1157999/1‑65 (MQ=255) tgtgtgTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGaa > 1:294053/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:1521880/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:44470/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:3750494/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:1256590/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:1390969/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:702743/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:2749502/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:1902208/1‑65 (MQ=255) tgtgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAc > 1:2134268/1‑65 (MQ=255) gtgtTCCCCCCCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAAcc < 1:249018/65‑1 (MQ=255) tgtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCg > 1:1929937/1‑65 (MQ=255) gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa > 1:1564072/1‑65 (MQ=255) gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa > 1:1762364/1‑65 (MQ=255) gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa > 1:644519/1‑65 (MQ=255) gtTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGa > 1:2128741/1‑65 (MQ=255) aCCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATgg > 1:849636/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg > 1:1091704/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg > 1:3851169/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg > 1:3600288/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg > 1:3005351/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg > 1:1684948/1‑65 (MQ=255) ccACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTg < 1:1010764/65‑1 (MQ=255) | TCGGCAAACTGCTGCGTCAACTCGGCGAATCAACTCAGGTGATGTGTGTTACCCACCTGCCACAAGTCGCGGGATGTGGTCATCAACACTATTTTGTCAGCAAAGAAACCGATGGTG > NC_000913/2753212‑2753328 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |