Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_pgi | 2,383,959 | 0 | A | T | 26.1% | 20.5 / 15.0 | 23 | intergenic (+35/+36) | elaD/yfbK | protease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein |
Reads supporting (aligned to +/- strand): ref base A (10/7); new base T (2/4); total (12/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCT > NC_000913_3_pae_pgi/2383811‑2384107 | aGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCACAAGCTGAAGCAGCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAata > 1:528107/1‑149 (MQ=255) aaaaTTGTAAATTAGAAAAAGATGGCGATCAAAATAACAATAAGTATAATAGAAACTTATAAATAAAACACAAACAAACAAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAAATAAGTCGACCACTCAACAActat < 2:605518/148‑1 (MQ=255) gTTATTTAGATAAATATGGCGATACAAATGGCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAAACCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGAc < 2:298601/149‑1 (MQ=255) ttttttAGATAAAGAGTGAAATGCAAAAATAAATCATTAGTATAGAAACTTAGAAATAAAACACTCACAACCCAAGAATAGAGCATCCGGCAAAAGAGTGAGTGAACATCATTGTTAATTAATGCGACCAGTCAGCAAATATGCTGACt < 2:357436/146‑1 (MQ=255) ttATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt < 2:389408/149‑1 (MQ=255) tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTg < 1:459363/149‑1 (MQ=255) aTAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc > 1:574516/1‑149 (MQ=255) tatGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTGAACATCAAAGTTAATGAATACGACCAGTCATCAAATATGCTGACTAGTGGAACATAgg > 1:150645/1‑148 (MQ=255) gCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCCTCCGGCAAAAGAGTCAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCCAACATAGtttcc > 2:405539/1‑147 (MQ=255) aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTCACATCATTGTTAATTAATTCGCCCAGTCAGCAAATATGCTAACTGGTCGAACATAGTTTGCTAAACAGTCATCac > 2:599409/1‑149 (MQ=255) ttAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcaca > 2:391724/1‑149 (MQ=255) acacTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAg > 2:376242/1‑149 (MQ=255) ccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCCCACCATCCGCCAGTTCAATCAGGc > 2:205107/1‑149 (MQ=255) aaGAGTGAGTCAACATCATTGTTAATTAATTCGACCAGTCAGAAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAca < 2:232843/149‑2 (MQ=255) aCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc < 1:341501/149‑1 (MQ=255) aCATCATTGTAAAGTAAGTCGACCAGACAGCAAATATGCCGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc < 1:150254/149‑1 (MQ=255) atTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca < 1:265434/86‑1 (MQ=255) atTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca > 2:265434/1‑86 (MQ=255) ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACGGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGc > 2:365070/1‑149 (MQ=255) ttCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCCTCACTGACTGATGTCAGTCACACCAGCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGCATATTCACAATTTGCCTGCTCAg > 2:293993/1‑149 (MQ=255) tCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCACCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGc < 1:599409/149‑1 (MQ=255) ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc > 2:271962/1‑149 (MQ=255) ccAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc < 1:365079/149‑1 (MQ=255) aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatca > 2:542751/1‑50 (MQ=21) aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatca < 1:542751/50‑1 (MQ=21) atGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATctgct < 1:523119/149‑1 (MQ=255) | AGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCACTGTTTGATCTGCT > NC_000913_3_pae_pgi/2383811‑2384107 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |