Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 47.4% | 10.7 / 11.2 | 19 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (3/7); new base A (9/0); total (12/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380220‑2380478 | cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAAAAGGACatttttt > 2:471508/1‑137 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:223038/138‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:574733/138‑1 (MQ=255) gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAAAGAAATGCACATTTTATTTTTTCTCAACTGtttt > 2:223038/1‑137 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGACATGCACATTTTATTTTTCCTCCTCTGTTTTTAGATGc > 1:240948/1‑139 (MQ=255) ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAACGACAAGCACAATTTTTTTTTCCTCATCTGTTTTTATATGAGAAAc > 2:680742/1‑138 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataa < 1:216113/138‑1 (MQ=255) tataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGCACAATTTATTTTTCCTCATCTGTTTTTAAAAGGGGAAAATGAAAATCAACGGGGGAAAAATAAAAGGGAAt > 1:410491/1‑138 (MQ=255) aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAtt < 1:201290/138‑1 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGAACATTTTATTTTTCCTCATCTGTTTTTAGATGGGAAAAATGAAATTTAAAGGAGGAAAAATAAATTGGGATGATTGGc > 1:668792/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTAAGAAGCGGAACATGAAATTCAACGGAGGAAAAATAAATGGGGATGATTTGCt > 1:305989/1‑139 (MQ=255) caGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACAATTTATTTTTCCTCATCTGTTTTTAGATGCGGAACAAGAAAATCAACGGGGGAAAAATAAATGGGGATGATTTGCTATTTTTCAGa > 1:152481/1‑138 (MQ=255) caGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGGGACATGCACATTTTTTTTTTCCTCTTCTGTTTTTAGTTGGGAAAAATGAAATTTAACGGGGGAAAAAAAAATTGGGATGATTGGGTATCTTGctttt > 1:67227/1‑135 (MQ=255) aataatAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTATTTTCTCTCCTCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGGAAAAAAAATGGGGATGATTGGCTATTTTTCAGAAAATGa > 2:60503/1‑139 (MQ=255) aGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa < 2:773770/139‑1 (MQ=255) ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa < 1:471508/139‑1 (MQ=255) ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa < 2:61628/139‑1 (MQ=255) aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 1:99723/1‑139 (MQ=255) aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 2:440292/1‑139 (MQ=255) | CTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380220‑2380478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |