Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,815,801 | 0 | A | . | 55.9% | 102.2 / 149.7 | 59 | intergenic (‑33/+62) | pyrE/rph | orotate phosphoribosyltransferase/truncated RNase PH |
Reads supporting (aligned to +/- strand): ref base A (12/14); new base . (16/17); total (28/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCG > NC_000913/3815752‑3815847 | cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa < 1:2197199/50‑1 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa < 1:5850053/50‑1 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa < 1:5596221/50‑1 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa < 1:1501550/50‑1 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa > 1:3360447/1‑49 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa > 1:3470240/1‑49 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa > 1:5366160/1‑49 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa > 1:200245/1‑49 (MQ=255) gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa > 1:3967517/1‑49 (MQ=255) tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaa > 1:3945511/1‑49 (MQ=255) atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaa > 1:1598036/1‑49 (MQ=255) atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaa > 1:4087775/1‑49 (MQ=255) atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaa < 1:4556259/48‑1 (MQ=255) tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaa < 1:5632079/50‑1 (MQ=255) tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaa > 1:3162679/1‑49 (MQ=255) atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaag > 1:5333896/1‑49 (MQ=255) tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaag < 1:3485646/49‑2 (MQ=255) tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaaa < 1:14916/49‑1 (MQ=255) tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg < 1:1928223/50‑1 (MQ=255) tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg < 1:3739298/50‑1 (MQ=255) tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGa < 1:4875834/50‑1 (MQ=255) tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGa > 1:2104446/1‑50 (MQ=255) cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGAc < 1:1001208/49‑1 (MQ=255) cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt > 1:3564071/1‑49 (MQ=39) cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt > 1:102877/1‑49 (MQ=39) cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt > 1:1621152/1‑49 (MQ=39) tGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtcat > 1:5803440/1‑50 (MQ=255) ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcat < 1:1320118/49‑1 (MQ=255) ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcat > 1:4814592/1‑49 (MQ=255) ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtcatc < 1:2822232/49‑1 (MQ=39) ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAg > 1:3789351/1‑50 (MQ=255) ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGt > 1:645619/1‑50 (MQ=39) cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAg < 1:4584027/48‑1 (MQ=255) cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc > 1:2836329/1‑50 (MQ=255) cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc > 1:1512493/1‑50 (MQ=255) cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc > 1:4426550/1‑50 (MQ=255) cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc > 1:2415147/1‑50 (MQ=255) cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGt < 1:2884471/48‑1 (MQ=39) cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGt < 1:394392/48‑1 (MQ=39) cGCTCCTCATCTTACATTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTc < 1:5475638/49‑1 (MQ=38) gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc < 1:2227132/49‑1 (MQ=255) gCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGc > 1:1572503/1‑50 (MQ=39) gCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGc > 1:5178191/1‑50 (MQ=39) ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCtt < 1:1198928/50‑1 (MQ=255) ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCtt < 1:838867/49‑1 (MQ=39) ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCtt > 1:5291758/1‑49 (MQ=39) ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCtt > 1:3454819/1‑49 (MQ=39) ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCt < 1:1912780/48‑1 (MQ=39) tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCtt > 1:1334639/1‑49 (MQ=255) tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCtt > 1:3215391/1‑49 (MQ=255) cATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaa > 1:2614473/1‑49 (MQ=39) aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa < 1:2457693/50‑1 (MQ=255) aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaa > 1:674199/1‑49 (MQ=255) aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaaa < 1:2864252/50‑1 (MQ=39) aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaaa < 1:5424434/50‑1 (MQ=39) aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaa < 1:4245437/49‑1 (MQ=39) aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaa < 1:3874187/49‑1 (MQ=39) aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaa < 1:5211811/48‑1 (MQ=39) tCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaaaa > 1:3472669/1‑50 (MQ=39) ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt > 1:607459/1‑50 (MQ=255) aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc < 1:327838/49‑1 (MQ=255) aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt < 1:4807183/47‑1 (MQ=39) aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc < 1:4119023/48‑1 (MQ=39) aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg > 1:4964794/1‑50 (MQ=39) cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa < 1:5918793/49‑1 (MQ=255) cTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGt < 1:1946491/50‑1 (MQ=39) ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGt > 1:5579337/1‑50 (MQ=255) ttttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt > 1:650984/1‑50 (MQ=39) tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt > 1:1815490/1‑50 (MQ=255) tttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt > 1:2455007/1‑50 (MQ=39) tttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt > 1:628988/1‑50 (MQ=39) tCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg < 1:5198185/49‑1 (MQ=39) cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa < 1:1657416/49‑1 (MQ=255) cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg < 1:5192785/50‑1 (MQ=255) cAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa < 1:3441307/48‑1 (MQ=39) aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc < 1:1971223/50‑1 (MQ=255) gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg < 1:4000858/50‑1 (MQ=255) | CGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCG > NC_000913/3815752‑3815847 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |