Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 6,994 | 0 | A | C | 55.0% | 12.5 / 17.7 | 20 | G322G (GGT→GGG) | yaaJ | sodium:alanine symporter family protein |
Reads supporting (aligned to +/- strand): ref base A (8/1); new base C (0/11); total (8/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGCACCCCAGGAACCCA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGC > NZ_CP009273/6911‑7080 | tCAGCACCCCAGGAACCCA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGcc > 2:238725/1‑90 (MQ=255) tCAGCACCCCAGGAACCCA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGcc > 2:238724/1‑90 (MQ=255) cAGCACCCCAGGAACCCA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCa > 1:564741/1‑90 (MQ=255) ccccccAGAACCCCA‑TTAGCCCCCGCACCCCCTTCGGAAAAAGCTAAATCCCTTCCAGCGGCATGTAGGTTGTGCCGTTCCCCGCCAGTa < 1:206048/88‑1 (MQ=255) cccccCGGAACCCATTTAGCCCCCCCTCCCCCTTGGGAAAAAGCTAAATCCCTCCCGGGGGCATGTAGGTTGTGCGGTTCCCCGCCAGTa < 1:41361/89‑1 (MQ=255) cAGGACCCCA‑TTAGCCCCCGCATCGCCTTCTGGATAAGCTAAATCCCTCCCAGCGGCATGTGGGTTGTGCCGTTACCCGCCAGTAATATc < 2:232540/90‑1 (MQ=255) ccA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATATAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATgg > 1:303346/1‑90 (MQ=255) aGCCCCCGCATCGCCTTCGGGAAAAGCTGAATCCCTCCCGGGGGCATGTGGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCtt < 1:354135/90‑1 (MQ=255) cGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGcctgt > 1:524624/1‑87 (MQ=255) acccccTTCTGAAAAAGCTAAATCCCTCCCAGGGGCATGTGGTTGGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAg < 1:42918/86‑1 (MQ=255) aTCGCTTTCTGAAAAAGCTAAATCCCTCCCAGCGGCATGTGGGTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAg < 2:133251/90‑1 (MQ=255) ttCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGcc > 1:120660/1‑90 (MQ=255) ccAGCGGAATTTAGGTTGTCCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCaa < 2:11908/90‑1 (MQ=255) gggaTGTGGTTTGTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc < 2:556224/87‑1 (MQ=255) aTGTAGGTTGTGCCGTTACCCGCCAGGAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGGTAAAAATGCCAATCATCTgg > 1:279088/1‑90 (MQ=255) gggTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACaa < 1:225920/89‑1 (MQ=255) gggTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACaa < 2:640339/89‑1 (MQ=255) aGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACaa < 2:102781/90‑1 (MQ=255) ttGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATcc > 1:149037/1‑90 (MQ=255) gTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTgcgc < 1:104703/90‑1 (MQ=255) | TCAGCACCCCAGGAACCCA‑TTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGC > NZ_CP009273/6911‑7080 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |