Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,349,621 | 0 | A | C | 26.7% | 82.9 / 32.0 | 45 | intergenic (‑149/‑66) | inaA/yfaH | putative lipopolysaccharide kinase InaA/putative uncharacterized protein YfaH |
Reads supporting (aligned to +/- strand): ref base A (33/0); new base C (0/12); total (33/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.48e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCTTTTGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGC > NC_000913/2349561‑2349685 | ttCTTTTGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAg > 1:4463074/1‑65 (MQ=255) ttCTTTTGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAg > 1:1004884/1‑65 (MQ=255) tCTTTTGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGc > 1:1178103/1‑65 (MQ=255) ttttGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGcca > 1:270445/1‑65 (MQ=255) gCCTCTACTAGCTTGCCCCCTCTTTTCCCCTTCATTGGCGGCGCA‑TATATCCCCCCAGccactac < 1:1307416/65‑1 (MQ=255) ctctactaGCTTGCCCCCTCTTTTCCCCTTCATTTGCAGCGAA‑TAGTTCCCCCCAGccactacca < 1:4531669/63‑1 (MQ=255) ccctactgCCTTGCCCCCGTTTTTCCCCTTCATTTGCGGCGCA‑TAGATCCCCCCAGccactacca < 1:4656849/63‑1 (MQ=255) cACTCCTACCTTACCCCCGTTTTTCCCCTTCTTTTGCAGCGAA‑TAGTTCCCCCCAGccactacca < 1:3290960/65‑1 (MQ=255) cACTACTAGCTTGCCCCCTTTTTTCCCCTTCTTTGGCAGCGCA‑TAAATCCCCCCAGccactacca < 1:87368/65‑1 (MQ=255) cACTACTAGCTTGCCCCCCCTTTTCCCCTTCATTGCCAGCGCA‑TAGATCCCCCCAGccactacca < 1:4662576/65‑1 (MQ=255) cACTACTAGCTTGACCCCGCTTTTCCCCTTCATTGGC‑GCGGAGAACATCCCCCCAGccactacca < 1:1314348/65‑1 (MQ=255) cACTACTAGCTTCCCCCCGTTTTTCCCCTTCATTTGCAGCGCA‑TAGATCCCCCCAGccactacca < 1:4608794/65‑1 (MQ=255) ctactaGCTTGCCCCCCCTTTTCCCCTTCATTGGCAGCGCA‑TAGATCCCCCCAGccactaccact < 1:152760/65‑1 (MQ=255) ctaGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAAt > 1:2899493/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCCCTAATAGTg > 1:1929526/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAggg > 1:3581876/1‑63 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:3148710/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:4098153/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:2975798/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:2962689/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:3007451/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:2190686/1‑65 (MQ=255) cTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTg > 1:1558933/1‑65 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTgcg > 1:1698298/1‑65 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTgcg > 1:2747950/1‑65 (MQ=255) cccacCGTTTTTCCCCTTAGTTAGCAGCGCA‑TAGTTCCCCCCAGCCACTACCACTAATAGTgcgc < 1:4630198/63‑1 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 1:2569492/1‑65 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 1:1903220/1‑65 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 1:4065986/1‑65 (MQ=255) acCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCca > 1:410901/1‑65 (MQ=255) acCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCca > 1:3400190/1‑65 (MQ=255) cttttttCCCCTTTATTGGCAGCGCT‑TATCTCCCCCCAGCCACTACCACTAATAGTGCGCCcagc < 1:958906/62‑1 (MQ=255) gCTTTTCCCTTGCATCTGCAGCAAA‑TATTTTCCCCCAGCCACTACCACTAATAGTGCGCCcagca < 1:672178/65‑1 (MQ=255) ttttCCCCTTATTTGGCGGCGAA‑GAGTTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGa < 1:2493835/65‑1 (MQ=255) ttttACTCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGa > 1:721252/1‑65 (MQ=255) ttttACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGa > 1:1286851/1‑65 (MQ=255) cccTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATaaaa > 1:4247898/1‑65 (MQ=255) ccTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAAt > 1:2206171/1‑65 (MQ=255) ccTTCATTGGCACCCCA‑TAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAAt < 1:1583324/65‑1 (MQ=255) cATTCGCAGCGCA‑TAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGa < 1:254581/65‑1 (MQ=255) tttGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAAt > 1:2837944/1‑65 (MQ=255) gCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTa > 1:2365668/1‑65 (MQ=255) gCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTa > 1:4595636/1‑65 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGa > 1:3540496/1‑65 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATa > 1:2787451/1‑65 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATa > 1:4232631/1‑65 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATa > 1:3850543/1‑65 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATAc > 1:527237/1‑65 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAAGTTAATCAGGCGGATAc > 1:4197085/1‑65 (MQ=255) accaGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGc > 1:2125588/1‑65 (MQ=255) | TTCTTTTGCCACTACTAGCTTGACACCGCTTTTACCCTTCATTTGCAGCGCA‑TAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGCAGATAAAATGTGAATTTAATCAGGCGGATACAGC > NC_000913/2349561‑2349685 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |