| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,826,926 | 0 | C | . | 31.8% | 48.7 / 21.1 | 22 | coding (1334/1344 nt) | astB | N‑succinylarginine dihydrolase |
| Reads supporting (aligned to +/- strand): ref base C (9/6); new base . (4/3); total (13/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATC > NC_000913/1826811‑1827065 | ggCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATa > 2:149036/1‑140 (MQ=255) cGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATaa > 2:210569/1‑140 (MQ=255) ccATCGCCCAGCCAGCGCCAGCGAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCa > 2:226428/1‑140 (MQ=255) ccagcgccagcgAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTc > 1:43423/1‑140 (MQ=255) gcgccagcgAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAAt < 1:226428/140‑1 (MQ=255) gAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAg < 2:43423/140‑1 (MQ=255) cGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTc > 2:311606/1‑140 (MQ=255) tttCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGccc > 1:271793/1‑140 (MQ=255) tcgctcgGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTcgc < 1:149036/140‑1 (MQ=255) ccGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCgg > 2:19970/1‑140 (MQ=255) cccGTTAAGGTCAGAGCAAGAAAATTATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATc > 1:211882/1‑140 (MQ=255) ttAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGgc < 1:175752/140‑1 (MQ=255) agCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCgg < 1:205713/129‑1 (MQ=255) agCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCgg > 2:205713/1‑129 (MQ=255) agCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAgg > 1:178050/1‑140 (MQ=255) gAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGgca > 1:118535/1‑126 (MQ=255) gAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGgca < 2:118535/126‑1 (MQ=255) aaaTTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATcgc < 1:19970/140‑1 (MQ=255) ttATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGAt > 1:306366/1‑133 (MQ=255) ttATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGAt < 2:306366/133‑1 (MQ=255) ttATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATcgcg > 2:270013/1‑137 (MQ=255) ttATCCATTG‑CCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATcgcg < 1:270013/137‑1 (MQ=255) cccccTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATc > 1:135371/1‑140 (MQ=255) | GGCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAACGCCGTTGATTTCTCGCTCGGTGATAACCGGTTTTTTACCCGTTAAGGTCAGAGCAAGAAAATTATCCATTGCCCCCTCCCTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATC > NC_000913/1826811‑1827065 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |