Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,171,271 | 0 | G | C | 35.6% | 78.5 / 15.4 | 55 | V150V (GTG→GTC) | CJ019_02060 | IS1182 family transposase ISPpu16 |
Reads supporting (aligned to +/- strand): ref base G (24/5); new base C (0/16); total (26/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.15e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATG > ACB143/2171162‑2171397 | actgTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCACCCTTTTCTCTGAGTCTATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTAAACAATCGCGACCGTAACTTCACCCa < 6:200202/133‑1 (MQ=255) ggTAAACCAATTCGAAGCTTCTCCCTTCTGTGCGTGACCCTATGCCTTAA‑CATCAACCTTTTATATTAGTCGATAATTCCAATCGACGGCAGCACTTTCAAAGCAGTCAACAAACGCGACCGTAACTTCACCCAgg < 2:296508/136‑1 (MQ=255) ccATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGt > 5:217840/1‑136 (MQ=255) cATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGt > 3:27231/1‑135 (MQ=255) gTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCCCCCAGGCCAAGGTGAAGGc > 8:81398/1‑136 (MQ=255) gTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGc > 1:282236/1‑136 (MQ=255) gTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGc > 3:200252/1‑136 (MQ=255) gTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CAGCAACCTTTTCTCTGAGTCGTTCATTGCCATAGACGGCAGCAAGTTCAAAGCAGTGACCCATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGc > 8:306832/1‑136 (MQ=255) tctcCCGTCAGTTCGGCAGCCGCTGCCGTAA‑CCTCAACCTTTTCTCTTAGTCGATCATTCCCAGCCACGGCAGAAAGTTCAACGCCGTCAACAAACGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCAt < 6:50341/136‑1 (MQ=255) cTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑ACTCAACCTTTTCTCTGAGTCGATCATTGCCATAGACGGCATCAACTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGCTGAAGGCCCGCATg > 5:162077/1‑136 (MQ=255) gCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATAGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATg > 1:261304/1‑134 (MQ=255) gCCGTCAGTTCGTGATCCTATGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATg > 8:159249/1‑134 (MQ=255) gatgTTCGTGTTCCTCTCCTGTAA‑CATCAACGTTTTATATGTTTCTTTCATTGCAATCGAGGGACTCATTTTCAACGCTGTTAACAATCGCGACCGTAACTTCACCCAGGCCAAGTTGAAGGCCCGCATGCAGCaa < 4:245189/133‑1 (MQ=255) ttCGTGATCCTCTGCCGAAA‑CCTCAAACTTTTCTCTGCGTCGATCATTGCCAGCGACGGCAGAAAGTTCAAAGCCGTCAACAAACGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCg < 3:196325/136‑1 (MQ=255) tCGTTACCCTCTG‑CGTAACCCACAAACTTTTATCTTAGGCGATCATTGCCATCGAGGGCAGAAAGTAAAAACCATTTAAAAAGCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCg < 2:118729/135‑1 (MQ=255) tCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATc > 5:36484/1‑134 (MQ=255) gtttaTCGTCTGCCTAAC‑CCTCAACCTTTTCTTTCAGGCCATAATTCCACTCGACGGCAGAATGTTCAAACCAGTCAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCga < 7:294835/131‑1 (MQ=255) cGTGATCCTCTGCCTTAA‑CCTCAACATTTTATCTGAGTCGATATTTGCCATCGACGGCAGCAAGTTCAAAGCAGTTAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCga < 4:210950/135‑1 (MQ=255) cGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACACAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgag > 7:191209/1‑136 (MQ=255) tcgcctGCCTAAA‑CCACAACTTTTTCTCTGAGTCGATATTGGCCATCCAGGGCATCAAGTTAAAACCAGTCAAAAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcaga < 6:100131/131‑1 (MQ=255) tCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCTATCATTGCCATCGAGGGCAGCAAGTTCAAAGCAGTCAACAACCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcaga < 6:154870/135‑1 (MQ=255) tCCGCTGCCGTAC‑CATCAACTTTTTCTGTTCGTATTTCATTCCCAGCGCGGGCTGCAAGTTCAAAGCAGTTAACAACCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcaga < 7:95546/135‑1 (MQ=255) ctgtGCCGAAA‑CCTAAACCTTTTCTCTGTGTCTATCATTGCCATCGACGGCAGCAAGTTCAAGGCAGTCAAAAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcagagc < 7:164632/132‑1 (MQ=255) ctGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt > 3:254068/1‑136 (MQ=255) gcggTAC‑CCTCAACCTTTTCTCTGAGTCCAGCTTTCCCAGCGACGCCTCCAAGTTAAAAGCCGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 6:336123/133‑1 (MQ=255) gcctaaa‑CCTCAACCTTTTCTGTGAGTCGAACATTCCCACCGACGGCTGCAAGTGCAAAGCTGTCAACAAGCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 3:326261/131‑1 (MQ=255) gCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt > 1:58261/1‑134 (MQ=255) gCCCTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 8:47555/136‑1 (MQ=255) cccGAAA‑CCACACCATTTTCTCGGAGTCTATCATTCCCATCCACGGCACCAAGATAAAACCCGCAAAAAAACGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 7:145752/135‑1 (MQ=255) cgttaa‑cctaaccTTTTTCTCTGAGGCTATCTTTGCAATCGCCGGCATAATTTACAAAGCAGTCAAAAAGCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 2:282236/132‑1 (MQ=255) cGTAA‑CCTAAACCTTTACTCTGAGTCTATCATTGCCATCGACGGCAGCAATTTCAAACAAGTCAACAAACGCGACCATAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa < 2:134170/134‑1 (MQ=255) tAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAgcg > 5:91317/1‑135 (MQ=255) tAA‑CCTCAACCTTTGCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAgcg > 1:180118/1‑135 (MQ=255) acccTCAACCTTTTCTGTGCGTCGATCATGGCACTGGAGGGCAGCAAGTTCACAGCTGTCACCAAACGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCt < 7:145148/134‑1 (MQ=255) tCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACAGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCt > 5:183544/1‑136 (MQ=255) aaaCCTTTTATCTGCGTCTATCTTTCCCAGCCCCGGCACACAGTTAAACGCTCTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCtt < 8:61791/135‑1 (MQ=255) aCCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTa > 2:36342/1‑135 (MQ=255) cTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAg > 6:291823/1‑136 (MQ=255) cTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAg > 1:323721/1‑136 (MQ=255) ctaTCATTCCCAGCGAGGGAAGCAAGTACAAACCATTCAAAAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGc < 5:134305/134‑1 (MQ=255) cGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg > 1:323694/1‑135 (MQ=255) aTTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGAcc > 4:215637/1‑136 (MQ=255) aTTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGAcc > 7:26797/1‑136 (MQ=255) ttcccATCGACGGCAGAAAGTTCAAAGCAGTTAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCg < 2:284403/133‑1 (MQ=255) cATCGACGGCAACAAGTTCAAAGCAGATAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAAACGACCAGAGCATAGAGCGCTACCTTACAGCGATGGATGCGGCAGACCGGGc > 4:242882/1‑135 (MQ=255) aTCCACGGCAGCAAGTCCAAAGCATTCAAAAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCaa < 1:75121/136‑1 (MQ=255) aCGGCAGCAAGTTCAAAGCAGTTAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAAACGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGcc > 5:200058/1‑136 (MQ=255) cGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGc < 7:49962/120‑1 (MQ=255) cGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGc > 8:49962/1‑120 (MQ=255) ggCACCAAGTACAAACCAGTAACAAACCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCgg < 2:186492/136‑1 (MQ=255) gcagcaAGTGCAAAGCTGCGAACATCCGCGACAGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGa < 2:74734/136‑1 (MQ=255) aaaGTTCAACGCAGTCAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTg < 7:81398/135‑1 (MQ=255) ttttCACAGCCGTCAACCACCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTg < 1:122386/132‑1 (MQ=255) gTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGc > 8:34157/1‑135 (MQ=255) caacccAGTCAACAACCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATg < 7:263310/131‑1 (MQ=255) | ACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAA‑CCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATG > ACB143/2171162‑2171397 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |