| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,643,967 | C→G | 18.2% | G113A (GGC→GCC) | ydfU ← | Qin prophage; uncharacterized protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,643,967 | 0 | C | G | 18.2% | 47.1 / 4.0 | 22 | G113A (GGC→GCC) | ydfU | Qin prophage; uncharacterized protein |
| Reads supporting (aligned to +/- strand): ref base C (8/10); new base G (2/2); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01 | |||||||||||
ATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGC > NC_000913/1643824‑1644113 | gttCCCGCGTCTGCGTACCGCGTATTGCGTGACGGATGCCGTCAATCACCCATGCTGTCAGATTTTGTTGCGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGATa < 2:400661/148‑1 (MQ=255) tCCCGCTCCTGCGTACCGCTTATTGCGTGACGAATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTAcc > 1:341545/1‑149 (MQ=255) tCCTGCGTACCGCGTAGTGGGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACg < 2:449281/149‑1 (MQ=255) cccTGCGTACCGCTTATGGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACg < 2:369776/148‑1 (MQ=255) ccGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGAc > 1:606429/1‑149 (MQ=255) ttATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGCGATCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCAATGCGCCGGTACACAAACGACGGATGACTGtt > 2:108071/1‑149 (MQ=255) gACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGa < 1:559318/149‑1 (MQ=255) caCGGATGACGTCAATCCCCCATGCGGTCAGATTTTGGTGAGCAAGCCGCTCCAGTGATTCCCATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGa < 1:581637/148‑1 (MQ=255) cAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAAAACTAGTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCa < 1:136482/149‑1 (MQ=255) aCCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAAAGACTGTTTCAGTGTGATGGTAATCGCCATTAgg < 1:443846/149‑1 (MQ=255) gttgAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGGCGGTGAAATGACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACa > 1:678818/1‑149 (MQ=255) gttgAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGGCGGTAAATTAATTGGGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACa > 1:678859/1‑149 (MQ=255) ttgAGCGAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTTCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACAt < 2:35557/149‑1 (MQ=255) agcagtTCCCGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGGAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTaa < 2:647056/149‑1 (MQ=255) cTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCaa < 1:27113/149‑1 (MQ=255) tGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGAACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCGATAGGCCCCTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGcccaa > 2:452145/1‑146 (MQ=255) aGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCGAGTCGGACAATGCGCCAACGCCGCCCGCAGCACGACTCACCcgtccg > 1:328190/1‑149 (MQ=255) acaacaCCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGGTa > 2:388777/1‑149 (MQ=255) acCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTaaa > 1:27576/1‑149 (MQ=255) acCATTGCGCCGGTACCATAACGGTGAATGACTGGTTCAGTGTAATGGTAATCGCGATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCGGCAGCACGAATCACCCGTACGTTACTaaa > 1:664491/1‑149 (MQ=255) cttGCGGCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTaaaaaa < 2:678859/148‑1 (MQ=255) aTTGCGGCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTaaaaaa < 2:678818/149‑1 (MQ=255) cgcCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGAAGTAACATGACGTAATAGCCAGTCGGACAATGCGCCACCGCCGCCAGCAGCACGAATCACCCGGTCGTTACTAAAAAACGgc > 1:3535/1‑149 (MQ=255) | ATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGC > NC_000913/1643824‑1644113 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |