Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 58.3% | ‑3.4 / 32.2 | 36 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (7/8); new base A (21/0); total (28/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.13e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAA > NZ_CP009273/2380208‑2380463 | taaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACa > 1:1182889/1‑139 (MQ=255) cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttat < 2:1182889/139‑1 (MQ=255) ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACAACATGCACattttttttttt > 2:846037/1‑137 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 2:1006635/138‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:507308/138‑1 (MQ=255) cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGCACATTTTATTTTTCCTCATCTG‑tt > 1:1192950/1‑139 (MQ=255) cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTG‑tt > 1:257297/1‑139 (MQ=255) aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTTTTTTTCCTCATCTG‑TTTTTA‑GATGCG‑AAAAAt > 1:149108/1‑138 (MQ=255) agggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACggagg < 1:20392/137‑1 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTATTTTTCCTCCTCTG‑TTTTTA‑AATGGG‑AAAAATGAAAAT‑TAAAGGAGGaaaaaa > 2:805162/1‑137 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACATGGACATTTTTTTTTTTCTCCTCTG‑TTTTTA‑AATGCG‑AAAAATGAAATT‑TTACGGGGGGAAAATAAATtgg > 2:1203837/1‑139 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGGCATGCACATTTTATTTTTCCTCCTCTG‑TTTTTA‑GATGCG‑AAAAATGAAAAT‑TAAAGGGGGGAAAATAAATggg > 1:1179022/1‑139 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATggg < 2:1254768/139‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATACACATTTTATTTTTTCTCCTCTG‑TTTTTA‑GATGCG‑AAACATGAAATT‑CAACAGAGGAAAAATAAATgct > 2:20392/1‑137 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTATTTTTCCTCATCTT‑TTTTTA‑GATGCG‑AAACATGAAAAT‑CAACGGAGGAAAAAAAAATggg > 1:459097/1‑139 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTG‑TTTTTA‑GATGCG‑GAACATAAAATT‑TAACGGGGGGAAAAAAAATggg > 1:1058979/1‑138 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTG‑TTTTTA‑GAAGGG‑GAAAAAGAAAAT‑TAACGGGGGGAAAAAAAATggg > 1:187000/1‑138 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAAGATGGACCTTTTATTTTTCCCCCTCTG‑GTTTAA‑GACGCG‑GAAAAAGGAAAT‑AAAGGGGGGAAAAATAAATagg > 1:497039/1‑138 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATg < 1:27003/138‑1 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTG‑TTTTTA‑GAAGCG‑AAACAAGAAAAT‑TCACGGGGGGAAAAAAAATTGGGATGATTggg > 2:1061068/1‑137 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGACAAGCACATTTTTTTTTTTCTCATCTT‑TTTTTA‑GATGGG‑GAAAATGAAATT‑CAACGGAGGGAAAATAAATTGGGATGATTTGCt > 2:651045/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCATAAAAAACGACATGTACAATTTTTTTTTCCCCCTCTG‑TTTTTA‑GAAGCG‑GCACAAGGAATT‑CCACCGAGGAAAAAAAAATGGGAATGATTTGCt > 1:591846/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACCTTTTTTTTTTTCTCCTCTG‑TTTTTA‑GATGCG‑AAACATGAAATT‑TAACGGAGGGAAAATAAATTGGGAAGATTGGCt > 2:24545/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTATTTTTCCTCACTTG‑TTTTAA‑GAAGCG‑AAAAAAGAAATT‑CAACGGAGGAAAAATAAATGGGAAGGATTggg > 2:910881/1‑137 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTATTTTTCCCCATCTG‑TTTTAA‑GAAGCG‑AAAAATGAAAAT‑TAACGGAGGAAAAATAAATGGGGATGATTGGc > 1:264697/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGAATACCAAAAAAGAGATAAGAAATTTATTTTTTCTTATCTGTTTTTAAGGGTGCA‑AAATA‑AAAATT‑CAAGGGAGG‑AAAATAAATTGGGattatttgg > 2:553419/1‑135 (MQ=255) gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAACACATGCCCAATTTATTTTTCCCCATCCT‑GTTTAA‑GGTGCG‑AAAACTGCAATTCCAACGG‑GGGAAAATTAATGTGAATATTTGGCTAt > 2:1278198/1‑139 (MQ=255) tAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAAGCACCTTTTATTTTTCCTCATCTG‑TTTTAA‑GATGGG‑AAACATGAAATT‑CAACGGAGGAAAAAAAAATGGGGATGATTGGACATTTTTAAGAAAtt > 1:166710/1‑139 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTG‑TTTTTA‑GATGCG‑AAAAATGAAATT‑TAAAGGAGGGAAAATAAAAGGGAAAGATTGGGTATTTTTCAGAAAtttaat > 1:1189083/1‑132 (MQ=255) taaAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTATTTTTCCTCATCTG‑TTTTTA‑GATGCG‑AAACATGAAAAT‑TAACGGAGGAAAAATAAATGGGAATGATTTGCTATTTTTCAGAAAATGACt > 2:162722/1‑136 (MQ=255) taaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTG‑TTTTTA‑GATGCG‑AAAAAAGAAAAT‑TAACGGGGGAAAAAAAAATGGGAATGATTGGCTATTTTTTAGAAATTGACt > 2:674565/1‑136 (MQ=255) taaAAATTTCTCGAGGCGGGAGCGAAAAAAGCGGGATGGACATTTTTTTTTACCTCCACTG‑GTTTAA‑AAAGAGCAAAAATAAAATA‑AAAC‑GAAGAAAAATAAAATGGAAAAATTGGCTATTTTGGAGaaaattnnn > 1:412179/1‑130 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGAAGAAATTGACTCAGaa < 2:766142/139‑1 (MQ=255) ttCTCGAGGCGGGTCCCGAAAAAGCGACATGCACATTTTATTTTTCCTCCTCTG‑TTTTTA‑GATGCG‑AAACATGAAATT‑TAACAGAGGGAAAAAAAATGGGGATGATTGGGTATTTTTTAGAAAATGAATCAGaaaaaaa > 2:890255/1‑136 (MQ=255) aaaaaGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt < 2:1162820/111‑1 (MQ=255) aaaaaGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt > 1:1162820/1‑111 (MQ=255) gCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa > 1:676426/1‑123 (MQ=255) gCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa < 2:676426/123‑1 (MQ=255) | TAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTG‑TTTTAA‑GATGCG‑AAACATGAAATT‑CAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAA > NZ_CP009273/2380208‑2380463 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |