| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,606,094 | 1 | . | G | 20.0% | 42.1 / 15.1 | 20 | noncoding (72/87 nt) | leuV | tRNA‑Leu |
| Reads supporting (aligned to +/- strand): ref base . (8/8); new base G (1/3); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.91e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATC > NC_000913/4605845‑4606303 | gccaggATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑a < 1:59357/247‑1 (MQ=255) gTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGacac > 1:193598/1‑251 (MQ=255) aTCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg < 1:240885/251‑1 (MQ=255) gATGGGAGTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg < 2:66749/251‑1 (MQ=255) tAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcacactc > 2:83463/1‑247 (MQ=255) gTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAt < 1:218451/251‑1 (MQ=255) ccGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt > 1:64168/1‑251 (MQ=255) cGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAAAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtttttt < 2:185838/251‑1 (MQ=255) aTCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc < 2:100138/251‑1 (MQ=255) gTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc < 1:30313/251‑1 (MQ=255) ccATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATGGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGAAATCGCCCCGTTGGGTGCGAGGGGGGGGAATTGAAccccc > 1:28269/1‑250 (MQ=255) gTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCTTCTTACTTTATATGATATCGCATCGTGGGGTGCGAGGCGTGGGACTTTAACCCTCACGTc > 2:35082/1‑250 (MQ=37) ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCt < 1:185291/250‑1 (MQ=255) aaTGAATTTGAATCAGGGATGAAAAATCAATGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTAcc < 1:58499/251‑1 (MQ=255) aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACACCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAATCTTACTTTTTTTGATATCGCCTCGTTGGGTGCGCGGGGAGGGACTTTAACCTCTACGTCAGTAAGAACACCAAC‑ACATGGGGCTAGCGCGTCTACCAATTcctccc > 2:230668/1‑247 (MQ=25) gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc > 1:3660/1‑122 (MQ=255) gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc < 2:3660/122‑1 (MQ=255) aaaTCAAGGAAGAAACAAGAAAGGAAGGAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg < 1:159546/251‑1 (MQ=255) aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTCTTACTTTTTTTGATATCCCCTCGTTTCGGGCGAGTGGGGGGGCTTGCACCCCCACGTCCGTAAGAACAC‑ACCGACCTAGAGCTAGTTCGTGTACCAATTCCGACACCTACGccacgtgccaggaatgtt > 2:164678/1‑234 (MQ=18) gTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCACTTCCGCCACCTTCGCACCGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGCGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc > 2:256428/1‑217 (MQ=255) gTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGGTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc < 1:256428/217‑1 (MQ=255) | GCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAAC‑ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATC > NC_000913/4605845‑4606303 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |