Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,032 | 0 | T | G | 68.4% | 3.1 / 10.1 | 19 | V288G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (13/0); total (15/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.87e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG > NZ_CP009273/594945‑595121 | tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg < 1:743024/90‑1 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg > 1:139999/1‑89 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg > 2:164804/1‑89 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 1:625628/90‑1 (MQ=255) cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg > 2:426535/1‑88 (MQ=255) gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggg > 2:365959/1‑87 (MQ=255) gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggggg > 2:740957/1‑87 (MQ=255) gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa < 2:391247/90‑1 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATcccg > 2:386988/1‑88 (MQ=255) gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGGTCCGGCCAAAAACCCCCggaa > 1:170702/1‑87 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGG‑CAAAAAAGCCCGAGAAggggt > 2:670965/1‑86 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCGG‑CAAAAAAGCCCGAGAAGGGAt > 2:604186/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGG‑CAAAAAAGCCCGAGAAGGGGTcgccg > 1:198378/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGCCCGG‑CAAAAAAGCCCGGGAAGGGATggccc > 1:346953/1‑89 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGCTCCGG‑CAAAAAAACCCGAGAAGTGATcgccc > 2:353687/1‑89 (MQ=255) aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCCCGATCCCG‑CCAAAAAACCCGAGAAGGGAGCGccgccgcg > 1:375074/1‑88 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGGGATCGCCGCCGTGa > 2:323644/1‑90 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc < 1:336034/90‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGGGATCCGG‑CAAAAAAACCCCAGGAGTGAACGCCCCCCGGAAgggc > 2:628947/1‑88 (MQ=255) aaCTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGG‑CAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGGAgggaa > 1:322938/1‑87 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 1:353687/90‑1 (MQ=255) agTGGCGGGCGGGGGGGGGATCCTGCGATCCGG‑CAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGAGAAAATGGGAACGCTGaaaaa > 2:651372/1‑89 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGa < 1:783884/60‑1 (MQ=255) cgggcgggGGGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGa > 2:783884/1‑60 (MQ=255) gggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg > 1:786954/1‑90 (MQ=255) gggTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg > 1:425281/1‑90 (MQ=255) ggtggtGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:330823/90‑1 (MQ=255) ggtggtGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:175017/90‑1 (MQ=255) ggtggtGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:697122/90‑1 (MQ=255) tggtgATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg < 2:623082/90‑1 (MQ=255) | TGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGG‑CAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG > NZ_CP009273/594945‑595121 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |