| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,549,566 | 0 | T | C | 39.4% | 44.5 / 15.8 | 33 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
| Reads supporting (aligned to +/- strand): ref base T (20/0); new base C (0/13); total (20/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.74e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT > NZ_CP009273/3549478‑3549647 | ttttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTc > 2:205397/1‑90 (MQ=255) tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCa > 1:18941/1‑90 (MQ=255) cAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGa > 2:444456/1‑90 (MQ=255) aCGGTATTCCCCGGCCCCGGGTCGCGCGGCAGGTTATAAATATTTTTTTCATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGc < 1:382315/90‑1 (MQ=255) ccGCCCCCGGGCCGCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt < 2:151814/90‑1 (MQ=255) gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgga > 2:341118/1‑89 (MQ=255) gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc > 1:207738/1‑90 (MQ=255) gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc > 1:209449/1‑90 (MQ=255) gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc > 2:479077/1‑90 (MQ=255) gTCGCGCGGCGGTTTATAAATTTTCTGTTCATG‑AGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac < 1:217273/90‑1 (MQ=255) cgcgggGCGGTTTATAATTTTTCTTTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:48407/90‑1 (MQ=255) cgcgggGCAGGTTAAAATTTTTCTGCTCA‑GAGGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:520835/90‑1 (MQ=255) cAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa > 2:80110/1‑90 (MQ=255) cAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa > 1:375603/1‑90 (MQ=255) ttATAAATTTTCTGTCCA‑GGAGGGAAAACC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa < 2:375603/90‑1 (MQ=255) tGAATTTTCTGTTAAGGCGGGAAAAC‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt < 1:479077/90‑1 (MQ=255) ttCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa > 2:63683/1‑90 (MQ=255) ttCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa > 1:474458/1‑90 (MQ=255) tCTTTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc < 1:532592/90‑1 (MQ=255) tCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc < 1:126250/90‑1 (MQ=255) tCTGTTCA‑GGAGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGGACTTCTCCAc < 2:207738/90‑1 (MQ=255) tGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc > 1:221443/1‑90 (MQ=255) tGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc > 1:86340/1‑90 (MQ=255) gTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa > 1:271276/1‑90 (MQ=255) gTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa > 2:484404/1‑90 (MQ=255) tCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATc < 1:444456/90‑1 (MQ=255) cATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt > 1:503385/1‑90 (MQ=255) aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca > 1:323330/1‑90 (MQ=255) aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca > 1:190492/1‑90 (MQ=255) aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca > 2:76944/1‑90 (MQ=255) ggaaaaa‑c‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa < 2:156878/87‑1 (MQ=255) ggAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt > 2:378567/1‑90 (MQ=255) gcaaaa‑c‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt < 2:323330/88‑1 (MQ=255) | TTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT > NZ_CP009273/3549478‑3549647 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |