| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,316,729 | 0 | C | A | 55.0% | 19.1 / 29.2 | 20 | E159* (GAA→TAA) | trpE | anthranilate synthase component I |
| Reads supporting (aligned to +/- strand): ref base C (7/2); new base A (0/11); total (7/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.64e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATA > NZ_CP009273/1316653‑1316818 | accACCATCAGCTTTTAAGGAAAAAAAAAAAAAAAATAGGGGCATTTTTTTTCCGCTGACAGTTGCGGAAAATCTTAAAATCCCGCcaca < 1:883700/88‑1 (MQ=255) ccccAGCGTTCCAGCAAAAAAAAAAAAAAAATCAGGGCAGTTTTTTTCCGCTGACAGTTGCGGTAAATCTTCAAACCCCGCCACAAGGTc < 2:404094/87‑1 (MQ=255) caccAGTGTTCCAGAGAAATAAAAACAGAAATCAGGGTATTTATTTCCCCCAGACTGTTGAGGAAAATTTTAAAATCCCGCCACAAGGTc < 2:1020094/87‑1 (MQ=255) cATCGGCTTTCCAGCGAAAAAAAAAAAAAAACCAGGGAATTTTTTTCCCCCTGCCAGTTGGGGAAAATTTTAAAACCCCGCCACAAGGTc < 2:1241462/90‑1 (MQ=255) cATCGGCGTTTCAGCGAGAAAAAAACAGAAATCGGGGCATTTTTTTTCCGCTGACAGTTGCGGTAAATCTTAAAATCCCGCCACAAGGTc < 2:114217/90‑1 (MQ=255) ccAGCTTTTCAGGGAAAAAAAAAAAAAAATCGGGGCATTTTTTTCCCGCTGACAGTTGCGGAAAATCTTAAAACCCCGCCACAAGGTCAt < 1:154437/89‑1 (MQ=255) tttCAGCGAGATAAAAGCAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAgaga > 1:1241462/1‑90 (MQ=255) ttCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGaa > 2:1053200/1‑90 (MQ=255) aaaaacaaaaacCGGGGCATTTTTTTTCCGCTGACGGTTGGGGTAAATCTTAAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCga < 2:1058200/90‑1 (MQ=255) aaaaaaaaaaACCGGGGAATTTTTTTCCCGGTGACAGTTGGGGAAAATCTAAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCga < 1:662773/90‑1 (MQ=255) aaaaaCAAAAATCGGGGCATTTTTTTCCCGCTGACAGTTGGGGAAAATCTTAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCga < 2:963201/90‑1 (MQ=255) aaaaaaTCGGGGCATTTTTTTCCCGCTGACGGTGGCGAAAAATCTAAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAAc < 1:940047/87‑1 (MQ=255) gTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGttc > 1:668482/1‑90 (MQ=255) tCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc > 2:1274096/1‑90 (MQ=255) gCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGt > 2:1348479/1‑90 (MQ=255) aattCTTAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTc < 2:995577/87‑1 (MQ=255) aaaTCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTc > 2:111596/1‑90 (MQ=255) aaaTCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTc > 1:75511/1‑90 (MQ=255) aaaTCTAAAAACCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTc < 2:846154/90‑1 (MQ=255) ttttCAACCCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGc < 1:1499001/88‑1 (MQ=255) cAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATa > 2:746993/1‑90 (MQ=255) | ATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATA > NZ_CP009273/1316653‑1316818 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |