| 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.5% | 95.3 / 17.8 | 54 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (29/2); new base C (0/22); total (29/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.97e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC > NC_000913/2984295‑2984589 | tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt > 1:315058/1‑151 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg > 1:854554/1‑151 (MQ=255) cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑ccc > 2:320512/1‑151 (MQ=255) cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑ccc > 2:131484/1‑151 (MQ=255) cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑cc < 1:314196/150‑1 (MQ=255) aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACCCCCTGG‑CCCCCGCCGa < 1:724250/151‑1 (MQ=255) aTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGc > 1:391577/1‑151 (MQ=255) aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGtt > 1:460364/1‑151 (MQ=255) gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTg > 2:398274/1‑151 (MQ=255) gggAGCCCAGAAAATGTATGCCGACTTGCCAGCCAT‑ATTAATCAAAAATATAAATGACGGCAGATATTTATCGGAAAGAAGGAAAG‑GAAGACCTATTCTTACCGTTAAAGGGTCAAT‑CCCCCCCCTGGCCCCCCGCCGATACAAAGCGTTg < 2:377040/151‑1 (MQ=255) gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGcc > 2:170169/1‑149 (MQ=255) gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGcc > 2:929251/1‑149 (MQ=255) ccAGAAATGTTATGTCCCCTTACCCCCCATATTTTTTA‑AACATATTAATGCCGGCAGATTTTTATGGATAAGGAGAAAAT‑GCAGATCTTTTCGCCACTTAAAATTGTAAATCCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAg < 2:502201/149‑1 (MQ=255) ccAGAAAATGTATGTCGACTTCCCAGCCAACATTATTA‑AACAGTTAAATGGCGGCAGATATTAGACGG‑AAGGGGAAAATGGCGGGGCTAATGGTCCCGTCCAATGGACAGGCCCACCCCCGGG‑CCTCCGCCGATACAAAGCGTTGCCAGTc < 2:803477/151‑1 (MQ=255) aaaaaaTGTATGTCGACTGGCCACCCATCATTATTC‑AACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAAT‑GCAGAGTTTTTCGCCACGTTAAAGGGTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCt < 1:265155/149‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 2:568705/1‑151 (MQ=255) aaaaTGTATGTCGACTTGCCACCCATCATTATTC‑AACAAATTAATGACGGAAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGTTTTTCGTCACGTTAAATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCtt < 1:93633/150‑1 (MQ=255) aaTGTATGTCGACTTCCCAGCCATTATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGAAAGGAGAAAAT‑GCAGATTTTTTCGTCACGTTAAATGGTAAATGCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc < 1:159561/151‑1 (MQ=255) cgccTTGCCCTCCATCATTATTC‑AACAGATAATTGACGGCATATATTTATGGGTAGGGAGAAAAT‑GCAGAGTTTTTCGTCCCGTTCAAGGGTAAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt < 1:921496/145‑1 (MQ=255) ccagccaCCATTATTA‑AACAAATTAATGACGGAAGATATTTACCGGAAAGGAGGAAAG‑GCAGAGTTATTCGTCATTTTCAATGGTCAATCCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc < 1:391300/149‑1 (MQ=255) cagaccactttattC‑AAAAAATTAATGCCGGCAGATTGTTTGCGAGAGGGGGAAAAT‑GCAGACTTTTTCGTCACTTTAAAGGGCCAATGCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc < 2:148299/142‑1 (MQ=255) ccaTAATTATTA‑AACAGATTAAGGCCGGCAAATATTTATGGGTAAGGGGAAAAT‑GAAGGGTTTTTCGTCAGGTTAAAGGGTCAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa < 2:727313/146‑1 (MQ=255) atcatTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc > 2:931703/1‑151 (MQ=255) tcatTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc > 1:592620/1‑150 (MQ=255) tattC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa > 1:689033/1‑150 (MQ=255) attC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc > 1:862681/1‑151 (MQ=255) ttC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 1:939997/1‑151 (MQ=255) ttC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 1:518365/1‑151 (MQ=255) tacaaaaaaTTAATGACGGCAAATTTTTATGGGTAAGGAGAAAAT‑GCAGAGTTTTTCGTCACTTTCAATGGTAAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc < 1:407068/142‑1 (MQ=255) aaCAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAAG‑GCAGAGCTTTTCGTCACTTTAAATGGTCAATGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 1:672552/149‑1 (MQ=255) aaCAGATAAAGGACGGCATATATTTATCGGAAAGGGGAAAAT‑GCAGACCTATTCGTCACTTTAAAGGGTCATGCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa < 2:718374/146‑1 (MQ=255) aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga > 1:859491/1‑150 (MQ=255) gagggaaGATATTTATCTGAGAGGAGAAAAT‑GCAGATTTATTCCACGCGTTAAATTTACAAGGCCCCCCCCTCC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa < 2:856482/144‑1 (MQ=255) ggtaGATTTTTATGGGTAGGGGGAAAAT‑GAAGAGTTTTTCGTCGCGTAAATTGGTGAACGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc < 1:145439/139‑1 (MQ=255) ggaaGATATTTAGGGGTAAGGGGAAAAG‑GCGGGGCTGTTCGTCGGGTTCAGTGGTCAGGGCCCCGCCGGGG‑CCGCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa < 2:241378/144‑1 (MQ=255) ggCAGATTTTTATCGGTAAGGAGGAAAT‑GCAGAGTTATTCGTCACTTTCAATGGTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg < 1:757235/151‑1 (MQ=255) cAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 1:721411/1‑151 (MQ=255) tatTTATCGGTAGGAAGAAATT‑GAAGAGCAGTTCCTCTCGTACATTGGTTACGGCCCCCCCTTTG‑CCGCCGCCGATACAAATCGTTACCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTct < 2:180572/145‑2 (MQ=255) tatTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc > 2:554657/1‑151 (MQ=255) aTCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa > 1:83721/1‑149 (MQ=255) aTCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa > 2:518591/1‑149 (MQ=255) gaaaGGAGAAAAC‑GCAGGGCTTTTCTTCACGTTAAAGGTTCAAGCCCACACCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa < 2:3024/143‑1 (MQ=255) aaGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGTGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 2:727721/1‑150 (MQ=255) ggaAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg > 1:288326/1‑151 (MQ=255) ggaAAT‑GAAGATCTTTTCGTCACGTTAAATGGTCAATCCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg < 1:459033/151‑1 (MQ=255) gCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc > 1:390256/1‑151 (MQ=255) ttttCGTCACGTTAAATGGTCAATCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg < 1:421405/146‑1 (MQ=255) tattcttcgttttCAAGGGCAAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc < 2:491300/138‑1 (MQ=255) tATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 1:10225/1‑150 (MQ=255) ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc > 1:111941/1‑151 (MQ=255) gTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg > 2:743741/1‑151 (MQ=255) cACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 1:413865/1‑151 (MQ=255) gTTCAATGGTCAATGCCACACCCTGG‑CCCCAGCCGATACAAAGCGTTGCCAGTCCTGTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGACTATCGGGTGACCGAGTGCGATCGCGccac > 2:545731/1‑150 (MQ=255) ggTCAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac < 2:391418/150‑1 (MQ=255) acacCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc > 1:90347/1‑151 (MQ=255) | TATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTC‑AACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAAT‑GCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC > NC_000913/2984295‑2984589 |
| 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: ‑ |