| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,171,384 | 0 | A | C | 32.3% | 41.5 / 26.4 | 31 | K285Q (AAG→CAG) | insF5 | IS3 element protein InsF |
| Reads supporting (aligned to +/- strand): ref base A (11/10); new base C (4/6); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.04e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCT > NC_000913/2171190‑2171629 | cAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACCATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGGAGGATCAgagtggtgttgt > 2:81137/1‑240 (MQ=2)cAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttga < 1:54530/251‑12 (MQ=2) aaaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGAGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGCACAATTTGAAAACCAGAAACTCGCTTAGGCCTGTGTCCATATTACGGGGGTAGGCTCAgagaggtgttgat > 1:75219/1‑239 (MQ=11) tgctgcTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGGGTCCATATTACGTGGGTAGGATCAgagtggtgttgatgcaggc < 2:144932/250‑19 (MQ=11) gctgctACGATAATGCCTGCGTGGAAAGCTGCTTTCATTCGCTGAAAGGGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatgcaggcc < 1:30432/250‑20 (MQ=11) ggctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGGGGATCAtttt < 1:268958/250‑1 (MQ=38) gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatgcaggccagtg > 1:19009/1‑228 (MQ=11) gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCAtttt < 1:12105/250‑1 (MQ=40) ctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTAGCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtg < 1:261396/234‑8 (MQ=2) aCGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGAc < 2:300620/250‑1 (MQ=255) gATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACACTTTGAAAACAAGAACCTCGCTTAGGCCTGTGACCATATTACGTGGGTAGGATCAGTCCTCCCTTCGTGGATCATTTTGACTgc > 2:263924/1‑251 (MQ=255) aaTGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTgccgc > 1:97471/1‑251 (MQ=255) gTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTCCGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTAt > 1:298933/1‑251 (MQ=255) tcttTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGGGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCTCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGAGCATGTTGACGGCCGCCTGCGCTAGTCTGttttt < 2:32055/250‑1 (MQ=255) tGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGAGCAACGGAGTTTAATTATATCGCATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGACGCCCGCGCTATTCTGTTTTTTGGCGCGCTGGCg > 1:254939/1‑251 (MQ=255) gTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTg < 1:409049/251‑1 (MQ=255) cATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCttt > 1:325992/1‑251 (MQ=255) aGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGGTCAgagtggtgttgatgcagg > 2:149461/1‑154 (MQ=11) aGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtggtgttgatgcagg < 1:149461/171‑18 (MQ=11) gCCGGGAAATAATGCGGGCAACGGTGTGTAATTATATCGAATGTGATTACAAGCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtc < 1:249135/251‑1 (MQ=255) ataatGCGGGCAACGGTGTTTAATTATATCGAGTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCCTCCCTTCGTGGATCATTTTGACTGCAGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCGGGGTCCCGGTTTc > 2:37787/1‑251 (MQ=255) aaCGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtg > 1:80480/1‑134 (MQ=2) aaCGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAgagtg < 2:80480/138‑5 (MQ=2) cGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGTGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTACTGCCAGAc < 2:351173/251‑1 (MQ=255) ggTGTTTAATTATATCGAATGTGATTACAATCGGGGGCGGCGGCACAGTTGGTGTGGCGGACTCAGTCCGGAACAATGTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTTCGTGGGTAGGATCAGTCATCCCTTCGTGGATCAGTTTGACTGCCGCCTGCGCTATTCTGTTTGTTGGCGCGGGGGGGCTCTGAAATCGCTCTTTAAACCCGACTTCCTGCGGCTCAGTTGCACTATTACCAGAcc < 2:249891/251‑2 (MQ=18) tgtTTAATTAAATCGAATGAGATTACAATCGGTGGCGGAGGCACAGTTGGTGTGACTGACTAAGTCCGGAACAAATTGAAAACACGCACCTCGCATCGGACTGTGTACATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGGTCATTTTGAATGCCGCCTGCACCATTATGTTTTTTGGCGTGCTGGCGTGCTTGAATCGCTCTTTAAAACCGACTTCCTGGGTCACGGTTTTACTATTACCCGACatt > 2:381061/1‑251 (MQ=255) ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGTCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtctc < 2:203985/222‑1 (MQ=255) ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtctc > 1:203985/1‑222 (MQ=255) ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAAAAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGGGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCTCCGCACCGCCATGa > 1:300075/1‑250 (MQ=255) ccagtCGGGGGCGCCGGCCCAGTTGGTGTGGCGGCCTCGGTCCGAAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCTTATTACGTGGGTAGGATCATTCTTCCCTTCGTTGATCATTTTGACTGCCGCCGCCGCATTTCTTTTTTTGGGCGCGCTGTCGAGCTGAAATCCCCCTTTAAACCCGACTTCCTGGGTTTCGTTTTCACTTTTACCAGACATTATTTTCGCCGCTCCGGCATGaaa < 2:152553/246‑1 (MQ=18) aaaaCAAGAACCTCGATTAGGCCTGTGTCCATATTAAGTGGGTAGGATCAGTCATCCCTTCGTGGAACATTTTGACTGCCGCCTGCGCTATTATGTTTTTTGGCGCGCTGGCGAGCTGACATCGCTCTTTAAACACGACTTCCTGGGTCTCGGTTTCACTATTACCCGACATTATTTTCGCCGCACCGCCCTGAAAGCGCGTAAACACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCt > 2:387092/1‑251 (MQ=255) | CAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCT > NC_000913/2171190‑2171629 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |