Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 56,392 | A→C | Y240D (TAC→GAC) | lptD ← | LPS assembly protein LptD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 56,392 | 0 | A | C | 93.1% | 75.5 / ‑0.4 | 29 | Y240D (TAC→GAC) | lptD | LPS assembly protein LptD |
Reads supporting (aligned to +/- strand): ref base A (0/2); new base C (14/13); total (14/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-02 |
GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAAA > NZ_CP009273/56306‑56478 | gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt < 1:722117/90‑1 (MQ=255) acgacgATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggt < 2:708650/90‑1 (MQ=255) acgacgATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggt < 2:284888/90‑1 (MQ=255) acgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGCtggtggt < 1:333696/90‑1 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTCGTtggtggtggg > 2:206224/1‑89 (MQ=255) aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTCGTtggggggggg > 1:72925/1‑89 (MQ=255) gCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGgtgt < 1:382018/90‑1 (MQ=255) tataATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATCGTTTGGGGTGGGGTGCt > 1:472772/1‑90 (MQ=255) aaTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTgg < 1:462592/90‑1 (MQ=255) ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATCGTTGGTGGGGGTGTACTTGGCGTTc > 1:261567/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATCGGTGGGGGGGGGGTACTTGGGGGTCg > 1:10273/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTCGTTGGTGGGGGTGTACTTGGCGTTCg > 2:215915/1‑90 (MQ=255) gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATCGTTGGGGGGGGTGTACTTGGCGTTCg > 1:198442/1‑90 (MQ=255) aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGGGGGGGGGTACTTTGCGTTCGGGGTCAAGAAAcc > 2:449267/1‑90 (MQ=255) ttGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc < 2:72925/90‑1 (MQ=255) gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGGGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt > 2:109764/1‑90 (MQ=255) gcgcGATGTTCCAGTAATATGGCAGGTAGAACACAAAATCGTTGGTGGGGGTGTACTTTGGGTTTGGGGTCAAGAAAACAGAGAGAAGtt > 1:593533/1‑90 (MQ=255) aaTATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa < 1:667761/90‑1 (MQ=255) aaTATGGCAGGTAGAACTCAAAATCGGTGGTGGTGGGGTACTTGGGGTTTGGGGTCAAGAAAACAGGGCGACGTTTTTCACCCCCCCGcg > 1:207208/1‑89 (MQ=255) atatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa < 1:984/90‑1 (MQ=255) tatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAAc < 2:697205/90‑1 (MQ=255) tGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTg < 1:191499/90‑1 (MQ=255) gCAGGTAGAACTCAAAGTCGTTGGTGGTGGGGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa > 2:304102/1‑90 (MQ=255) aGGTAGAACTCAAAGGCGTTGGGGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa > 1:105637/1‑90 (MQ=255) gAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAggg < 2:636197/90‑1 (MQ=255) aCTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGc < 1:318089/90‑1 (MQ=255) aCTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGc < 1:304102/90‑1 (MQ=255) cAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc > 1:330763/1‑54 (MQ=255) cAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc < 2:330763/54‑1 (MQ=255) agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATaaa < 1:611495/90‑1 (MQ=255) | GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAAA > NZ_CP009273/56306‑56478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |