| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_hsa_tpiA | 3,657,601 | 0 | G | C | 25.0% | 37.3 / 16.1 | 25 | intergenic (+34/+385) | hdeD/arrS | acid‑resistance membrane protein/antisense sRNA ArrS, function unknown |
| Reads supporting (aligned to +/- strand): ref base G (7/11); new base C (3/3); total (11/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAT > NC_000913_3_hsa_tpiA/3657455‑3657747 | gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt < 1:394538/149‑1 (MQ=255) tCTGTAACGCTGCTTTCCACACTGGTGGGAATTGAACTGGTATTGAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt < 2:406126/149‑1 (MQ=255) aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATGTCGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct < 1:249055/149‑1 (MQ=255) cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa < 1:193511/135‑1 (MQ=255) cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa > 2:193511/1‑135 (MQ=255) tGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGCTCATAGTCACTTGATATAACTATGCTCCGATTAATActatct > 2:337764/1‑149 (MQ=255) aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTTATGTGACTATGATCCGATTCATACTCTCACCCCTACCAACTGTTGTaa > 2:6246/1‑148 (MQ=255) cgAACTTATATTGAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc < 1:417883/148‑1 (MQ=255) gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt > 1:422683/1‑105 (MQ=255) gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGCATCATAGTCACTTGATGTGACTATGAt < 2:422683/105‑1 (MQ=255) ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt > 1:54893/1‑137 (MQ=255) ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt < 2:54893/137‑1 (MQ=255) cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCAGAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAc < 2:263498/149‑1 (MQ=255) atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATc < 2:331682/149‑1 (MQ=255) ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCt > 1:56749/1‑149 (MQ=255) cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTTAACCTACTGAATTTCATCTAAtt > 1:112875/1‑149 (MQ=255) cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttt > 1:260939/1‑149 (MQ=255) ccGGACCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGCATTTCATCTAATTTTTTTCACTCTa > 2:300019/1‑149 (MQ=255) aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa > 2:48950/1‑61 (MQ=21) aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa < 1:48950/61‑1 (MQ=21) cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTaa > 1:188953/1‑129 (MQ=255) cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTaa < 2:188953/129‑1 (MQ=255) atcatCACTTGATGTGACGTTTATCTGATTAGCACCCTCTCCGTTACGAAGCGTTGTAGATCAATTTCCCACTATCATTGAAATAATTACCTTCTAGTGATTATTTCACTCTACTGAATTTCATCTAATTTTTTTCTCTCTATGGCaaa < 2:34860/145‑1 (MQ=255) tatgcACTTGGTGTGACTATGGTCCGATTAATACTTTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAt < 2:225778/145‑1 (MQ=255) ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTAACAATTTCAACCCACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAt < 2:265123/149‑1 (MQ=255) | GTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAT > NC_000913_3_hsa_tpiA/3657455‑3657747 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |