| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,357,755 | 0 | A | T | 53.3% | 28.0 / 65.4 | 45 | P530P (CCA→CCT) | yhcD | fimbrial biogenesis outer membrane usher protein |
| Reads supporting (aligned to +/- strand): ref base A (14/7); new base T (0/24); total (14/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.97e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGA > NZ_CP009273/3357679‑3357844 | aCTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCaa < 1:1146582/90‑1 (MQ=255) tGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAac < 2:298972/90‑1 (MQ=255) gggCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAac < 1:144355/89‑1 (MQ=255) gcagcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAacac < 2:156977/90‑1 (MQ=255) cagcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCa < 1:921414/79‑1 (MQ=255) cagcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCa > 2:921414/1‑79 (MQ=255) gcaTCAATTTAGTGTTTCAAAAAAAAATAAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGc < 1:272742/90‑1 (MQ=255) aTCAACTTAGTGCTTAAAAAAATAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca < 1:1035021/90‑1 (MQ=255) aaCTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCacac < 2:698895/90‑1 (MQ=255) aGTGCTTCAAAAAATAATGAAGACGAGGGTGATCTTTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtat < 2:959074/90‑1 (MQ=255) gTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtata > 1:735968/1‑90 (MQ=255) tCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGc > 2:781252/1‑90 (MQ=255) tCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGc > 1:1177043/1‑90 (MQ=255) aaaaaaaaATGAAGACGAGGGTGATTTTTTCTCCGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCa < 2:875150/88‑1 (MQ=255) aaaaaaTAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCa < 2:607833/88‑1 (MQ=255) ataatGAAGACGATGGTGATTTTTTCTACGTCAACATCAGTATTCTTTTTGGCAGATCAAACCCCGCCACACTGTATACTCGCCATGACg < 2:688279/90‑1 (MQ=255) aaaatAAAGACGATGGTGATTTTTTCTCCCACAACCTCAGTATTCCTTTTGGAAAACAAAACCCCGCCACACTGTATACTCGCCATGACg < 1:781252/88‑1 (MQ=255) aaaatAAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGACCAAACCCCGCCACACTGTATACTCGCCATGACg < 2:1231848/88‑1 (MQ=255) taatGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAAATCAAACCCCGCCACACTGTATACTCGCCATGACGa < 1:301276/90‑1 (MQ=255) aatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGAt < 2:735968/90‑1 (MQ=255) tGAAGACGATGGTGATCTTTTCTACGTCAACATCATTTTTCCTTTTGAAAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCa < 1:1157649/90‑1 (MQ=255) aGACGATGGTATTCTGTTCTCCTTAAACATCAGTATCCTTTTTGGAAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATaa < 1:741448/90‑1 (MQ=255) agggTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc < 2:539377/88‑1 (MQ=255) agggTGATCTGTTCTACGTAAACATCAGTATTCTTTTTGGAAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc < 1:6997/88‑1 (MQ=255) aTGGTGATCTGTTCTACGTAAACACCATTTTTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc < 2:607791/90‑1 (MQ=255) ggTGTTCTTTTCTCCGTCAACATCATTATCCCTTTTGGAAGAACAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCAc < 1:81598/90‑1 (MQ=255) ggTGATCTTTTCTACGACAACATAAGTTTTCCTTTTGGCAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCAc < 1:1054107/90‑1 (MQ=255) tttCTTTTCTCCGTAAACATCAGTTTTCCTTTTGGAAGAAAAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt < 1:168962/88‑1 (MQ=255) tttCTGTTCTACCCCAACACCATTTTCCCTTTTGGCAGACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt < 2:1164623/88‑1 (MQ=255) gctttGTTCTCCTCCAACATCATTTTTCCTTTTGGCAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt < 1:611340/86‑1 (MQ=255) gATCTGTTCTACTCCAACACCATTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATTCTCGCCATGACGATCATAAAACCCACTAt < 1:214530/90‑1 (MQ=255) gATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt > 2:1054107/1‑90 (MQ=255) cTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGt > 2:514453/1‑90 (MQ=255) gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc > 2:425220/1‑90 (MQ=255) gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc > 1:680939/1‑90 (MQ=255) aCGTCAACATCATTTTCCCTTTTGGCAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg < 1:825749/90‑1 (MQ=255) aCGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg > 1:1101006/1‑90 (MQ=255) aCGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg > 2:447090/1‑90 (MQ=255) aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATg > 2:450089/1‑90 (MQ=255) aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATg > 1:188336/1‑90 (MQ=255) aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATg > 2:1198068/1‑90 (MQ=255) atttttCCTTTTGGCAGACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc < 1:1244073/86‑1 (MQ=255) tCCTTTTGGCAAACCACACCCCGCCACCCTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGtctc < 2:1114552/90‑1 (MQ=255) tCCTTTTGGAAGCACAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGtctc < 1:837655/90‑1 (MQ=255) tCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGtctc > 2:685608/1‑90 (MQ=255) aTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCaga > 1:923563/1‑90 (MQ=255) aTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCaga > 2:1109003/1‑90 (MQ=255) aTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCaga > 1:1215296/1‑90 (MQ=255) | ACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGA > NZ_CP009273/3357679‑3357844 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |