| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,984,441 | 0 | A | C | 47.6% | 60.0 / 12.9 | 42 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (21/1); new base C (0/20); total (21/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.17e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG > NC_000913/2984293‑2984587 | aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc > 2:677157/1‑151 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg > 1:236062/1‑151 (MQ=255) tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg > 2:271862/1‑150 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg > 1:242764/1‑151 (MQ=255) aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa > 1:376750/1‑150 (MQ=255) aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa > 2:379537/1‑149 (MQ=255) gggggAGCCCAAAAAATGTATGTCGACTTGCCACCAATCATTAATAAACAGATTAATGACGGCAGATTTTTATCGGAAA‑GGAGAAAATGCAGGGCTTTTCCCCACGTTAAAGGGCCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCg < 1:835900/148‑1 (MQ=255) gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCg > 2:652079/1‑150 (MQ=255) ggACCCCAAAAAATGTTTTTCGTCTTGCCCCCCATTATTATTCAACAGATTAATGACGGCAAATTTTTATGGGTAA‑GGAAAAAAGGCAGAGCTTTTCCCCCCGTTAAATGGTCAATCCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGc < 1:130595/151‑1 (MQ=255) ccaaaaaaaTGTTTGTCGTCTCCCCAGCCATCATTATCCAAAAGATAAATGACGGCAGATATTTATCGGAAA‑GGAGAAAAGGCAGACTTTTTCGTCCCGTTAAATGATAAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGt < 1:352814/146‑1 (MQ=255) gAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCtt > 1:388048/1‑151 (MQ=255) aaaaTTTATGCCGCCTGCCAACCCATAATTATTAAACAGATTAATGCGGAAAGATTTTTACCGGTAA‑GGGGAAAGCGCAGAGCTTTTCGTCCCGTTAAAGGGTAAAGCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt < 2:183367/150‑1 (MQ=255) atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg > 2:593985/1‑151 (MQ=255) tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATATTTCTCGGTAA‑GGAGAAAAGCCAGAGCTATTCGTCACTTTCAAGGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc < 1:322167/149‑1 (MQ=255) gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg > 2:367974/1‑151 (MQ=255) gACTTGCCCGCCACCATTTTTCAACAGATTAATGCCGGAAAATTTTTATCGGAAA‑GGGGAAAATGCAGACTTATTCTTCATGTTAAAGGGCAAATGCCGCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc < 2:651095/150‑1 (MQ=255) ttcccacccaTCTTTTTTAAAAAGATAAATACCGCCAATTATTTAGCGGAAA‑GGAGAAAATGCAGAGTTTTTCGTCCCTTTCATTGGTAATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg < 1:800958/145‑1 (MQ=255) cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc > 1:554616/1‑151 (MQ=255) gccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATTCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGCCCCCCGCCGATACAAAGCGTTACCAGTCCTTTGCGGTCATTACGTTTCACCAt > 1:621222/1‑151 (MQ=255) tcatTATTCAACAGATAAAGGACGGCATTTTTTTACCGGTAA‑GGAGGAAATGAAGGGCTTTTCTTCCCTTTCAAGGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa < 1:389210/146‑1 (MQ=255) attattCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 2:615205/1‑151 (MQ=255) ttattCAAAAGATTAAGGACGGCAGATTTTTATCGGTAA‑GGAGGAAATGAAGAGCTTTTCGTCCCGTTCAATGTTAAACCCCCCACCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt < 2:450704/146‑1 (MQ=255) tattCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa > 2:227602/1‑150 (MQ=255) ttaaaaaGATTAAGGAGGGCAGATATTTATGGGTAG‑GGAGAAAATGCAGAGCTATTCGTCACGTTCAAGGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa < 2:26965/142‑1 (MQ=255) aCAGATTAATGGCGGCAGATATTTATGGGTAA‑GGAGAAAGTGCAGAGTTTTTCTTCCCGTAAAATGATCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt < 2:359771/142‑1 (MQ=255) acaataagcccGGCAGATTTTTA‑CGTTAACGGAGAAAATGCAGTTCTTTTCTTCCCGTTCAATGTTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga < 1:24027/140‑1 (MQ=255) aTTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc > 1:288328/1‑151 (MQ=255) aaTGACGGAAGATATTTATCGGTAG‑GGAGAAAATGCAGATCTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc < 2:587201/149‑1 (MQ=255) aaTGACGGAAGATATTTATCGGAAG‑GGAGAAAATGCAGACTTTTTCGTCACTTTAAATGGTAAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc < 2:430429/149‑1 (MQ=255) aGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 1:476708/1‑150 (MQ=255) aGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 1:413493/1‑150 (MQ=255) tatTTATCGGAAA‑GGAGAAAATGCAGACTTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg < 1:265516/145‑1 (MQ=255) atgggaaa‑GGAGGAAAGGAAGGCTTTTTCGCCAGGTCAAATGCTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa < 2:367375/144‑1 (MQ=255) ggTAA‑GGAGAAAAGGCAGAGTTTTTCGTCACTTTAAATGGTCAAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa < 1:39714/147‑1 (MQ=255) ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg > 1:69709/1‑151 (MQ=255) aaaTGAAGAGTTTTTCGTCAGTTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg < 2:669256/151‑1 (MQ=255) aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg > 2:683312/1‑150 (MQ=255) agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 1:435735/1‑149 (MQ=255) tATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 1:30412/1‑150 (MQ=255) cACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc < 2:709518/148‑1 (MQ=255) atgggTCAAACCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat < 2:108600/146‑1 (MQ=255) ccACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg < 2:234650/151‑1 (MQ=255) | AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAA‑GGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG > NC_000913/2984293‑2984587 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |