| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 20,600 | G→A | intergenic (‑92/+215) | insA1 ← / ← rpsT | IS1 family protein InsA/30S ribosomal subunit protein S20 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 20,600 | 0 | G | A | 90.0% | 22.0 / ‑5.8 | 10 | intergenic (‑92/+215) | insA1/rpsT | IS1 family protein InsA/30S ribosomal subunit protein S20 |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (5/4); total (6/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑CT‑‑‑‑‑‑TTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCC > NC_000913/20372‑20821 | gtgtAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCACAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGt > 1:327243/1‑249 (MQ=255) aaaaatatagtggtcttaagtttaatacctcaatgtgtatcacaatatccacattCTTTGTATGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATTTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgct < 1:358690‑M1/226‑88 (MQ=255) gtcttaagtttaatacctcaatgtgtatcacAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggcttt > 1:41731‑M1/14‑152 (MQ=255) tcttaagtttaatacctcaatgtgtatcaCAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttt > 2:133827‑M1/13‑151 (MQ=255) cttaagtttaatacctcaatgtgtatcACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggcttttt > 1:375565‑M1/12‑150 (MQ=255) taagtttaatacctcaatgtgtaTCACAATAGCCACATTCTGTGTAGGGG‑TTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgctggcttttt < 2:375095‑M1/240‑102 (MQ=255) ttaatacctcaatGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaa < 1:202983‑M1/246‑108 (MQ=255) aTCACAATATCCACATTCAATGAAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGACACGAAGAGTTCAAATAGTACAACGGCATGTCTGAAAGAACCTGCAACATTTGACTATCGGGCAACCGTACGTGTGTTCTCTACCGCCTCACAAAAGCAGCAGGCATAACCAAAACCGCTTGAGCGGGCTTATTCACCAAGCTACAGCAAATTGGCg > 2:304890/1‑250 (MQ=255) cacaATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGAGATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgat < 1:220783‑M1/250‑127 (MQ=255) tCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacacgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagcc > 1:27391‑M1/1‑117 (MQ=255) | GTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑CT‑‑‑‑‑‑TTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCC > NC_000913/20372‑20821 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |