Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 926,174 | A→T | 22.4% | intergenic (‑203/+51) | serW ← / ← infA | tRNA‑Ser/translation initiation factor IF‑1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 926,174 | 0 | A | T | 22.4% | 131.3 / 37.7 | 58 | intergenic (‑203/+51) | serW/infA | tRNA‑Ser/translation initiation factor IF‑1 |
Reads supporting (aligned to +/- strand): ref base A (23/22); new base T (7/6); total (30/28) | |||||||||||
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.53e-01 |
CATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATG > NC_000913/926074‑926306 | cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa < 1:888964/102‑1 (MQ=255) cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa > 2:888964/1‑102 (MQ=255) cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt < 3:520949/120‑1 (MQ=255) cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt > 4:520949/1‑120 (MQ=255) tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttcttt > 3:216545/1‑119 (MQ=255) tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttcttt < 4:216545/119‑1 (MQ=255) gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa > 3:575161/1‑141 (MQ=255) gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa < 4:575161/141‑1 (MQ=255) cAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa > 4:815076/1‑91 (MQ=255) cAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa < 3:815076/91‑1 (MQ=255) aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc < 3:928450/91‑1 (MQ=255) aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc > 4:928450/1‑91 (MQ=255) tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa > 3:915469/1‑102 (MQ=255) tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa < 4:915469/102‑1 (MQ=255) cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc > 1:992982/1‑114 (MQ=255) cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc < 2:992982/114‑1 (MQ=255) ccAAAGTACTTCATACATATCACAAGATAAAAAGGCCGGTTaa < 2:676590/43‑1 (MQ=255) ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa > 1:676590/1‑43 (MQ=255) aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt > 4:343852/1‑148 (MQ=255) aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt < 3:343852/148‑1 (MQ=255) tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa < 3:902519/93‑1 (MQ=255) tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa > 4:902519/1‑93 (MQ=255) cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc > 3:327175/1‑129 (MQ=255) cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc < 4:327175/129‑1 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg < 4:159990/62‑1 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg > 3:106917/1‑62 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg > 3:159990/1‑62 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg < 4:106917/62‑1 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt > 1:270672/1‑100 (MQ=255) caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt < 2:270672/100‑1 (MQ=255) aaaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAgg > 1:775260/1‑148 (MQ=255) aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGt > 4:740667/1‑36 (MQ=255) aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGt < 3:740667/36‑1 (MQ=255) aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg > 3:353003/1‑114 (MQ=255) aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg < 4:353003/114‑1 (MQ=255) aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg < 2:775260/148‑1 (MQ=255) aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg > 3:528917/1‑148 (MQ=255) tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg > 3:1095631/1‑87 (MQ=255) tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg < 4:1095631/87‑1 (MQ=255) tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt > 4:531771/1‑133 (MQ=255) tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt < 3:531771/133‑1 (MQ=255) taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCATTTCAACAGTCACTTTg > 1:82041/2‑131 (MQ=255) taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTg < 2:82041/130‑1 (MQ=255) aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg > 4:1090976/1‑109 (MQ=255) aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg < 3:1090976/109‑1 (MQ=255) aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc > 1:1105657/1‑53 (MQ=25) aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc < 2:1105657/53‑1 (MQ=21) aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc < 1:1046275/87‑1 (MQ=39) aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc > 2:744247/1‑87 (MQ=39) aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc < 1:744247/87‑1 (MQ=39) aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc > 2:1046275/1‑87 (MQ=39) aggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc < 2:190810/113‑1 (MQ=255) aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAg > 3:291480/1‑76 (MQ=38) aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc > 1:190810/1‑117 (MQ=255) ggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCa < 4:190118/111‑1 (MQ=255) ggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCa > 3:190118/4‑114 (MQ=255) gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa < 1:58648/58‑1 (MQ=255) gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa > 2:58648/1‑58 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa > 3:619580/1‑63 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa < 4:619580/63‑1 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg > 3:997656/1‑98 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg < 4:997656/98‑1 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTc < 4:308967/129‑1 (MQ=255) ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTc > 3:308967/1‑129 (MQ=255) | CATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATG > NC_000913/926074‑926306 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |