Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 41.7% | 50.3 / 14.3 | 36 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (17/4); new base C (0/15); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.21e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984294‑2984591 | atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc > 2:283429/1‑151 (MQ=255) cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑cc < 1:183463/150‑1 (MQ=255) gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccg > 1:158158/1‑151 (MQ=255) gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccg < 1:741023/151‑1 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccgccg > 2:444555/1‑151 (MQ=255) aaTGAAAAAATGAGGGAGCCGAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGAAATTTAGCGGTAAGGAGGAAATGCAGAGCTATTAGTCACGTTCAATGGACAAGGCCACACCCTGG‑CCCccgccg < 1:736225/151‑1 (MQ=255) aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGa < 2:614573/151‑1 (MQ=255) aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGtt > 2:577459/1‑151 (MQ=255) gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGcc > 1:293804/1‑149 (MQ=255) cccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGt > 2:372615/1‑151 (MQ=255) ccaaaaaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGGAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCa < 2:17168/144‑1 (MQ=255) caaGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AAGGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATCCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGt < 2:802283/148‑1 (MQ=255) aaaaTGTATGTCGACTTGCCACCCATCATTATTAAACAGATA‑AATGACGGC‑AGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCtt < 2:956535/150‑1 (MQ=255) tatgccaccTTCCCAGACATTATTATTAAACAAATT‑AAGGACGAC‑ATATATTTACCGGAAAGGAGGAAATGCAGAGCTTTTCGCCTCTTTAAAGTTTCAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg < 1:765670/140‑1 (MQ=255) tggctatttcccacccaTCATTATTCAAAAAATA‑AATG‑CGGCAAAATATTTATCGGTAGGGAGAAAATGCGAGGCTTTTCGCCCCTTTAAATGGTAAATGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg < 2:302640/138‑1 (MQ=255) cGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc > 1:764228/1‑151 (MQ=255) cGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc > 1:950227/1‑151 (MQ=255) cTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt > 1:381984/1‑151 (MQ=255) ttGCCACCCATCATTA‑TCCACAAATTAAATGACGGA‑AAATATTAATCGGAAAGGAGAAAATGCAGATCTATTCGTCATGTTAAATGGTCAACCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt < 1:191202/149‑1 (MQ=255) ccaTCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt > 2:739582/1‑151 (MQ=255) ccaTCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt > 2:318770/1‑151 (MQ=255) ataaaataaaagACGGA‑AGATATTTACGGGAAGGGAGGAAATGCAGTGCTGTTCGCCACGTTAAATTTCAAACGCCACCCTCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 1:925684/138‑1 (MQ=255) aaGGAAGAAAGGCAGCGTTATTCTCCGCTTTCAATGCTCAACGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg < 1:220061/147‑1 (MQ=255) aaaaGGCAGAGTTTTTCGCCACGTTAAATGGTCAACCCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt < 1:53730/150‑1 (MQ=255) aaatgaaaggCTTTTCGTCACTTCAAATGGTAAACGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg < 2:780476/142‑1 (MQ=255) aaaTGCAGGGTTTTTCGTCA‑TTTTAATGGGTAAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa < 1:656894/144‑1 (MQ=255) cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 2:142525/1‑151 (MQ=255) cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 1:902474/1‑151 (MQ=255) aTTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt > 1:111613/1‑151 (MQ=255) ttCGCCACTTTAAGTGGTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa < 1:319833/149‑1 (MQ=255) tCGTCACTTTCAATGGTCAATCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc < 2:98110/150‑1 (MQ=255) cACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 1:577633/1‑151 (MQ=255) cACGTTCAATGGACAGCCCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc < 2:111648/151‑1 (MQ=255) aCGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 1:597028/1‑150 (MQ=255) ttCAAGGGTCAATCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc < 1:445871/150‑1 (MQ=255) aTGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac > 2:975077/1‑151 (MQ=255) cacCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACg > 2:63040/1‑151 (MQ=255) acCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc > 1:874767/1‑151 (MQ=255) | ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATT‑AATGACGGC‑AGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984294‑2984591 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |