Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,571,380 | 0 | A | C | 53.8% | 46.5 / 65.8 | 52 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (17/7); new base C (0/28); total (17/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.58e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATA > NZ_CP009273/1571295‑1571465 | tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 1:584898/90‑1 (MQ=255) aCCAGTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCAt < 1:1305562/90‑1 (MQ=255) gTCCGAACGCGTAGTGCGATACCCCAATTC‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTCCCCCGCCATTGaa < 2:1194673/90‑1 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:1750573/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:550484/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:456222/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:394384/1‑90 (MQ=255) ccGACCGCGTAGTGCGATAACCCAATTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGc < 1:1406944/90‑1 (MQ=255) ccGAACGCGTAGTGCGATACCCCAATTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTCCCCCGCCATTGAAGc < 2:95590/90‑1 (MQ=255) ccGAACGCGTAGTGCGATACCCCAATTC‑CCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGc < 2:413908/90‑1 (MQ=255) cgcgTGGTGCGATACCCCATTTC‑CCCCTTCCTATCATCAGCGTTAAAGCTTTTGATTTTTGCATAAATTCCCCCGCCATTGAAGCGACCa < 2:468317/90‑1 (MQ=255) gCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAc > 2:1897390/1‑90 (MQ=255) cccaaaTTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGcac < 1:932412/86‑1 (MQ=255) cccAATTTCCCCCTTCCTACCATAGGCGTAAAAGCGTTTGCTTTTGGCATCCTTCCCCCCGCCATTGAAGCGACCAAATTCAACCGGcac < 2:1111374/90‑1 (MQ=255) cAATTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa < 1:1231172/90‑1 (MQ=255) cAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa < 1:450247/90‑1 (MQ=255) aTTT‑CCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 2:394384/90‑1 (MQ=255) aTTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:185075/90‑1 (MQ=255) aTTT‑CCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:537102/90‑1 (MQ=255) aTTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc > 1:1587500/1‑90 (MQ=255) aTTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:550484/90‑1 (MQ=255) ttt‑CCCCTTGCTACAATCCGTGTTGAAGATTTTAGTTTTTAGAGCACTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACt < 2:967692/90‑1 (MQ=255) ttt‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACt < 1:881844/90‑1 (MQ=255) tt‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTg > 1:1432840/1‑90 (MQ=255) cACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTc > 1:1180901/1‑90 (MQ=255) aCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCa > 1:280984/1‑90 (MQ=255) ccTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCAt > 1:1332441/1‑90 (MQ=255) ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 1:1750573/90‑1 (MQ=255) ttGCTACCAAAAGCTTTGAAGCTTTTGATCTTTACATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 2:41090/90‑1 (MQ=255) atcatcAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt < 2:462539/90‑1 (MQ=255) atcatcAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt > 2:960730/1‑90 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 2:796206/1‑90 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 2:563252/1‑90 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 1:1734647/1‑90 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 1:1128504/1‑90 (MQ=255) agtttttGATTTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa < 2:1128504/86‑1 (MQ=255) gTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc > 1:1487134/1‑90 (MQ=255) gTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc < 1:796206/90‑1 (MQ=255) ttGTTCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:1387496/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 1:529914/90‑1 (MQ=255) ttttGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:311997/89‑1 (MQ=255) ttttACATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:1864379/89‑1 (MQ=255) cTTTGCATCATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:1107625/90‑1 (MQ=255) cTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 1:703867/90‑1 (MQ=255) cTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:456222/90‑1 (MQ=255) tttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:351867/90‑1 (MQ=255) tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:546937/90‑1 (MQ=255) ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa > 2:996235/1‑90 (MQ=255) cAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 1:790575/90‑1 (MQ=255) cAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 2:280984/90‑1 (MQ=255) aaTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata < 2:726478/90‑1 (MQ=255) aaTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata < 2:232031/90‑1 (MQ=255) | TGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTT‑CACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATA > NZ_CP009273/1571295‑1571465 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |