Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,819,796 | T→C | 8.9% | intergenic (+33/‑442) | rpoS → / → PP_1624 | RNA polymerase sigma S factor/group II intron‑encoding maturase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,819,796 | 0 | T | C | 8.9% | 127.0 / 3.2 | 45 | intergenic (+33/‑442) | rpoS/PP_1624 | RNA polymerase sigma S factor/group II intron‑encoding maturase |
Reads supporting (aligned to +/- strand): ref base T (19/22); new base C (2/2); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01 |
GGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGG > NC_002947/1819659‑1819943 | ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGtttt > 6:169760/1‑149 (MQ=255) gCCTTACCCGTGAACGGGTGCGGCAAATCCATGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGttttt < 7:87173/149‑1 (MQ=255) cTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTtgt < 6:119047/149‑1 (MQ=255) cTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTtgt < 6:274230/149‑1 (MQ=255) ttACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTtgtg > 7:250295/1‑149 (MQ=255) ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGAGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTgcg > 4:312810/1‑149 (MQ=255) gggTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCAt < 2:16813/149‑1 (MQ=255) ggCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTATCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCggg < 5:65264/120‑1 (MQ=255) ggCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAACTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCggg > 6:65264/1‑120 (MQ=255) gCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGt > 6:260794/1‑149 (MQ=255) tCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTgcgc < 1:76978/149‑1 (MQ=255) gAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCC‑GTGATCACCGGGGCGtttttt > 1:264643/1‑114 (MQ=39) gAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCC‑GAGATCACCGGGGCGtttttt < 2:264643/115‑2 (MQ=38) aaGGGCTCAAGCGCTTGCGCGAGATCCTTGAAATGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGAGCTGTAACAAAAAACGCCCAGGTGATC‑TCGGGTCGTTTTTTTGTGCGCCAGGCATgg > 7:232172/1‑130 (MQ=255) tCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGc > 7:90720/1‑149 (MQ=255) aaGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt < 5:260794/149‑1 (MQ=255) aaGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt > 4:149526/1‑149 (MQ=255) gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg < 3:149526/149‑1 (MQ=255) gcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt < 5:315672/149‑1 (MQ=255) ttGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt > 3:109532/1‑149 (MQ=255) gcgcGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg < 5:169760/148‑1 (MQ=255) cTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt < 7:30231/148‑1 (MQ=255) aaaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt < 6:34472/148‑1 (MQ=255) aaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCACGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTTCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa > 4:221706/1‑149 (MQ=255) aaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc > 3:67851/1‑149 (MQ=255) tGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa < 4:67851/149‑1 (MQ=255) gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg < 8:90720/148‑1 (MQ=255) gAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc < 4:109532/149‑1 (MQ=255) aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca > 4:35811/1‑148 (MQ=255) aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca > 2:68210/1‑148 (MQ=255) tGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGAGCAACCCGCTTGGTTGTGGAGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacaca > 3:248843/1‑148 (MQ=255) cAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt < 4:107813/104‑1 (MQ=255) cAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt > 3:107813/1‑104 (MQ=255) gTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa < 3:312810/149‑1 (MQ=255) gcgcTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa < 1:68210/148‑1 (MQ=255) gcgcTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg > 7:290423/1‑149 (MQ=255) cgcTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt < 3:35811/149‑1 (MQ=255) tAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa < 8:250295/148‑1 (MQ=255) aaaaaCGCCCC‑GAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg < 1:286747/94‑1 (MQ=11) aaaaaCGCCCC‑GAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg > 2:286747/1‑94 (MQ=11) cccGGTGATC‑TCGTGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt < 3:115216/131‑1 (MQ=37) cccGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt > 4:115216/1‑131 (MQ=255) cccGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAggg > 1:321870/1‑149 (MQ=255) aTC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac > 2:61589/1‑109 (MQ=38) aTC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac < 1:61589/109‑1 (MQ=38) c‑tcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCggg > 2:101884/1‑149 (MQ=34) c‑tcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg < 2:321870/148‑1 (MQ=34) c‑tcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg < 4:248843/148‑1 (MQ=34) c‑tcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTAGCggg > 7:284361/1‑149 (MQ=21) | GGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGG > NC_002947/1819659‑1819943 |
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: ‑ |