Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,296,583 | T→C | 28.6% | intergenic (‑336/‑230) | nhaA‑I ← / → PP_1133 | Na(+)/H(+) antiporter NhaA/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,296,583 | 0 | T | C | 28.6% | 35.5 / 16.0 | 28 | intergenic (‑336/‑230) | nhaA‑I/PP_1133 | Na(+)/H(+) antiporter NhaA/transposase |
Reads supporting (aligned to +/- strand): ref base T (10/10); new base C (3/5); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01 |
CTTTCGCGGGTAAACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGGGTTCGTGTGCGCGGGACAGAGATCGTAGCCAGCACTGGTGTGTGTCTGCACACCCGCGTAGGAGAGTGATCCGCACGGTATCGGCGTGAAGT > NC_002947/1296435‑1296712 | cTTTCGCGGGTAAACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGt > 4:230421/1‑149 (MQ=255) tCGCGGGTAAACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTg < 2:128582/148‑1 (MQ=255) gagcgGGTTTACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTgg > 3:321581/3‑149 (MQ=255) gcgGGTAAACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGc < 4:84606/148‑1 (MQ=255) ccGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGc > 3:119342/1‑148 (MQ=255) ccGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGc > 3:34888/1‑148 (MQ=255) cTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAgcc < 3:230421/147‑1 (MQ=255) cAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtt < 4:34888/149‑1 (MQ=255) cAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCtt < 4:19799/149‑1 (MQ=255) gggACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTc > 7:57244‑M1/1‑148 (MQ=255) acaaTTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCtt < 4:271069/146‑1 (MQ=37) ccATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttt < 8:57244‑M1/149‑6 (MQ=255) cATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTg > 3:217095/1‑149 (MQ=38) aTTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGAGTACGGCGTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt > 5:20829/1‑149 (MQ=255) ttGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGt > 8:89005/1‑148 (MQ=35) gtgtAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaaaatgtc < 5:202357‑M1/149‑16 (MQ=255) tgtAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTcttttt > 5:243224‑M1/1‑133 (MQ=255) tgtAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTcttttt < 6:243224‑M1/139‑7 (MQ=255) tAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaaaatgtca < 7:237482‑M1/148‑17 (MQ=255) gaTCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaaaa > 8:248060‑M1/1‑126 (MQ=255) gaTCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaaaa < 7:248060‑M1/137‑12 (MQ=255) ttCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaa < 8:186138‑M1/123‑10 (MQ=255) ttCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTctttttcaa > 7:186138‑M1/1‑114 (MQ=255) agGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTgctttttcaaaatgtcagcgtcgctgcggaactgaa > 4:342762‑M1/1‑33 (MQ=255) agGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTgctttttcaaaatgtcagcgtcgctgcggaactgaa < 3:342762‑M1/69‑37 (MQ=255) gCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcatgtttgctgcg > 2:359984‑M1/1‑31 (MQ=255) gCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcatgtttgctgcg < 1:359984‑M1/90‑60 (MQ=255) gTTGGCCCGGTAGGGCCAGCCGCCTTTTctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcatgtttgctgcgcgacagcgcagcccggaagggctgggtgatctcctacagggccc < 1:262856‑M1/131‑104 (MQ=255) gTTGGCCCGGTAGGGCCAGCCGCCTTTTctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcatgtttgctgcgcgacagcgcagcccggaagggctgggtgatctcctacagggccc > 2:262856‑M1/1‑28 (MQ=255) | CTTTCGCGGGTAAACCCGCTCCTACAGGGGACCATTGTTCCGTGTAGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCGCAGCAAACATGAGGTGTACGGCTTCAGTTCCGCAGCGACGCTGACATTTTGAAAAAG‑AAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTGCTTTGTTAGGGTTCGTGTGCGCGGGACAGAGATCGTAGCCAGCACTGGTGTGTGTCTGCACACCCGCGTAGGAGAGTGATCCGCACGGTATCGGCGTGAAGT > NC_002947/1296435‑1296712 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |