Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,561,205 | A→T | 9.9% | intergenic (+307/‑163) | yffL → / → yffM | CPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,561,205 | 0 | A | T | 9.9% | 107.9 / 5.0 | 44 | intergenic (+307/‑163) | yffL/yffM | CPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM |
Reads supporting (aligned to +/- strand): ref base A (14/24); new base T (2/2); total (16/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01 |
CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCA > NC_000913/2561068‑2561346 | cTAGGCGACGGGATTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg < 1:834164/143‑1 (MQ=255) cTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg < 1:661929/143‑1 (MQ=255) cTAGGCGACGGGAGTTTTGTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg < 1:404868/143‑1 (MQ=255) gCGACGTGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaa < 2:628092/143‑1 (MQ=255) tggcgtCCCGGTCGGTAATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactcac < 2:626052/138‑1 (MQ=255) cccGGGCGGTCATAGAGAGGCGCCTAGAGTGAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACaa < 1:801256/143‑1 (MQ=255) cccgtcggtcATCGGGAGTCGCCTAGAGTGACACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc < 2:587607/140‑1 (MQ=255) cATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc < 1:275514/143‑1 (MQ=255) cATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc < 1:818546/143‑1 (MQ=255) agagagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca < 1:315209/141‑1 (MQ=255) aTAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca < 1:477814/143‑1 (MQ=255) gTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACcaac > 2:629026/1‑143 (MQ=255) gTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGATCGGCTGATTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGATGGAAACGACCATTCACTCACGCACAACCTGAACAGCAAaacaac > 2:749939/1‑143 (MQ=255) gagagGTAAACAGTGGAAGTATTTCGCCGCCTGATAGGTCGGCGGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACca < 2:793386/142‑1 (MQ=255) agTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGAACACATTTCCATCGAGCTGTAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca < 1:846121/143‑1 (MQ=255) ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCCGACCCACAAACCATCATg > 1:117593/1‑143 (MQ=255) ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg > 1:492940/1‑143 (MQ=255) ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg > 1:9856/1‑143 (MQ=255) ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCAAGAAACGCCAATTCACTCAAGCACAACCTGAACAGCAGACCAACAAACCATCATg > 1:403509/1‑143 (MQ=255) aCAGAAGAATTATTTCGCCGCCTGAGAGGTCGGCTGGTTGCGAGCATTAAGCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAt < 2:384388/143‑1 (MQ=255) gccgccGGATAGGTCGGCTGGTTGCGAGCATTAACCAGGCTGATTAGCGCAATCATTTCCAGCCAGCTGGAAACGCCCATTCACTCACGCACAAGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcg < 1:177110/143‑1 (MQ=255) cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATACCTCAAAcgg > 2:806203/1‑143 (MQ=255) tGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATATCACAGGGt > 1:907082/1‑143 (MQ=255) gTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGGGCGCGGTATATCGCAAAgggcgc > 1:12769/1‑140 (MQ=255) gCATTCACCAGGCTGCGTAGCGCAATCACTTCCCGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATGTCTTa < 2:504409/143‑1 (MQ=255) cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGTAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAt > 2:584339/1‑143 (MQ=255) atgaaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt < 2:817032/140‑1 (MQ=255) aaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggt < 1:433855/143‑1 (MQ=255) cTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGCGGGATTTCTTATTCGTGGGccgt > 2:496821/1‑141 (MQ=255) tGATTAGCGCAATTATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCACACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc < 1:583606/143‑1 (MQ=255) gATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCAAGGGGAGGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCg < 1:893684/143‑1 (MQ=255) tcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg < 1:377462/142‑1 (MQ=255) gcgcAGTCATTTCCAGGAAGCTGGAAACGCCCATTCACTCACGCACATCCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg < 2:385649/143‑1 (MQ=255) ggcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATATCTTATTGGTGGGCCTTCGGAAGGCTTCg > 2:539191/5‑143 (MQ=255) tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCCCAAAGGGAGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa > 1:892820/1‑143 (MQ=255) ccATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa < 2:430755/142‑1 (MQ=255) gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa < 2:537635/115‑1 (MQ=255) gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa > 1:537635/2‑116 (MQ=255) aGCAAGGTGCAAACGCCCATTCACTACCACACAACATGAACAGCAGACCAACAAACCATCATGGCATTCGTTGATGGGGAGCGCGGTATATCTCAAATGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAgatcggaa < 2:100161/143‑8 (MQ=255) gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc < 2:350571/106‑1 (MQ=255) gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc > 1:350571/1‑106 (MQ=255) cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa < 1:611844/143‑1 (MQ=255) cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTCCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa < 1:83943/143‑1 (MQ=255) cAAGATGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGCATATCTCAAAGGGTGCTTGCTTAATGGTGGGCCGTCGCAAGGCTGCGAGCAATTGccc > 2:627675/1‑142 (MQ=255) | CTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCA > NC_000913/2561068‑2561346 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |