Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,340,076 | 0 | G | C | 30.4% | 25.2 / 14.3 | 23 | intergenic (+61/‑270) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base G (7/9); new base C (3/4); total (10/13) | |||||||||||
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.04e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT > NC_002947/4339928‑4340216 | ggggCCGCCGCTGAACTCCGCGGGGGGGGCGGGCTATGCGCTGGAGGGGGGGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACCCAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAg < 8:115205/149‑1 (MQ=255) gggCCGCCGCTGAAATACGGGGGGCTGGCGGGCTATTCGCTGGCGGTGGTGCAGAGCATCCCTTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAg < 8:120840/148‑1 (MQ=255) ggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAggg < 7:235330/149‑1 (MQ=255) ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAGGGTGaca < 1:91935/149‑1 (MQ=255) gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGCGGATGGATGGCTGAGAACACACAAGCAAGTGATTTAGCGGGGCCCTGTA‑GGAGCGGCTTTACCCCCGAAAGCGACGGTGACcccacc > 3:252547/1‑149 (MQ=255) gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCTCTGTA‑GGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg < 2:372468/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCATTTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg < 2:102783/149‑1 (MQ=255) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg < 2:171481/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg < 6:372548/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg < 1:15655/149‑1 (MQ=255) cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTCAGAACACACAAGCAATTGATTTAGCGGGGCCCTTTA‑GGAGCGGGTTTTACCCCCAATGGCTGGGTTACACCCCCCACGCCATACCGGGTAAACCCCCTACTACAGGGCCCGTgg > 7:365440/1‑149 (MQ=255) ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGGAGCGGGTTTAACCCCGAATGCGAGGGTGACACCCCCGACGCATTAGCGGGTAAACCCCCACCTACAGGGCCCTTGGGGGGTTTGACAGGGGGGAGTGAttg > 1:88081/1‑147 (MQ=255) caaCCCACAAGCAATTGCTTTAGCGGGGCCCTGGA‑GGAGCGGGGTTACCCGGTAATGGCCGGGTGACACCACCGACGCATGCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc < 4:252547/147‑1 (MQ=37) caaCACACAAGAAATTTATTTCGCGGGGCCCGGTA‑GGCGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGcccc < 2:88081/148‑1 (MQ=37) ggCCCTGTA‑GGCGCGGGTTTACCCGCTAATGCGAGGGTTACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATa < 6:358649/149‑1 (MQ=255) tGTA‑GGAGCGGGTTTACCCCCGAATGCGAGGGTGACACCCCCCCCGCATTAGCCGGTAAACCCCCTCCAACAGGGCCCATGGTGGGTTTTACAGGGCGGAGTGATGGCCGGGAGCCCCAGACCACGTGCCAAAAACTTGGAAAACCaa > 2:93578/1‑147 (MQ=19) cccGCTAATGCGAGGGTGACACCCCCGACGCATTCGCGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa < 8:102086/148‑1 (MQ=255) cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAt < 7:163114/127‑1 (MQ=255) cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAt > 8:163114/1‑127 (MQ=255) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTTGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 5:166128/122‑1 (MQ=18) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 5:67840/1‑149 (MQ=25) ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGATTTTGACAGGGAGGAGTGATGGCCGGTCGCCCCATACCTATTGCCAAAATATTg > 6:166128/1‑122 (MQ=17) ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 2:347073/1‑149 (MQ=255) cgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt > 8:249433/1‑149 (MQ=255) gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt > 7:72112/1‑143 (MQ=255) gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt < 8:72112/143‑1 (MQ=255) | GGGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTA‑GGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT > NC_002947/4339928‑4340216 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |