| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,340,089 | 0 | G | C | 25.8% | 48.7 / 12.9 | 31 | intergenic (+74/‑257) | 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 (12/11); new base C (3/5); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCG > NC_002947/4339952‑4340230 | gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATTGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTcgcg < 2:253870/149‑1 (MQ=25)gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTcgcg < 5:387910/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTCGCgg < 1:382382/149‑1 (MQ=25) tGGCGGGCTATGTGCTGGTGGTGGTGCAGAGCATCCCGTTCGACCCTGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCggg < 6:145224/149‑1 (MQ=255) gggggCTATGAGCTTGATGTTGTGCAGATCATCCCGTTCGACCCGGGGTTTGTTTGCCTCACAACACACAAGACTTTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCggg < 1:121307/146‑1 (MQ=25) gCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCggg < 3:174119/147‑1 (MQ=255) tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTTACCGCGAATTC‑GAGGGTG‑‑ACACCACCGCCGCATTAGCGGGGTAACCCg > 1:93635/1‑149 (MQ=255) tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGGTTAACCGCGCATGG‑GGGGGGG‑‑ACCCCCCCGACCCCATCGCGGGTTATCCCg > 6:313066/1‑149 (MQ=255) cgCTGGAGGTGGTCCAGAGCATCCCGTTCGACCCGGGG‑CTGCATGGCTGAGAACACACAATCAAATCATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGCGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGcgc < 1:296173/147‑3 (MQ=38) aGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTc < 5:126565/148‑1 (MQ=255) aGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTc < 1:251289/148‑1 (MQ=255) gagCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCCAAGGC‑GACGGTG‑‑ACACCACCCACGCCATAGCCGGTAAAACCGCACCAACACGGCCCgtggtg > 4:47981/1‑149 (MQ=37) tCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTTACCGCGCATTCGGCGGTTGCAACACCGCCGAGGCAAACGCGGGAAAACCCCCACAAAAAGGGCCATTGGTggt > 1:322277/1‑147 (MQ=255) tCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTCGCgg > 2:239873/1‑115 (MQ=18) tCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTCGCgg < 1:239873/115‑1 (MQ=18) cccGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑CCACACCCGCCGCATTCGCgg > 6:183640/1‑105 (MQ=38) cccGGGG‑ATGGAGGGCTGAGACCACACAAGCAATTGATTTAGCGGGGCCCTGTCGGA‑GCGGGTGTACCCGCGAATGC‑GCGGGGG‑‑ACACCACCGACGCATTCGCgg < 5:183640/105‑1 (MQ=255) ccGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTcgc < 3:372467/102‑1 (MQ=17) ccGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GTCGGTG‑‑GTGTCACCCTCGCATTcgc > 4:372467/1‑102 (MQ=17) cGGGG‑ATGGATGGCTGACAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGC‑GGGGGTTTACCCGCTAAGGC‑GAGGGTG‑‑ACACCACCGACGCATTCGGGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAg < 3:47981/149‑1 (MQ=38) gCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGGCGGGTTTAACCGCGAATGC‑GAGGGTG‑‑ACACCACCGCCGCAATAGCGGGTAAACCCGCGCCTACAGGGCCCTTGGTGGGTTTTACAGGGCGGAGTGATGGCCggg > 1:75095/1‑148 (MQ=37) a‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa > 2:211473/1‑149 (MQ=255) ccGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGGCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 7:171295/1‑149 (MQ=255) cgcgAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc > 6:379639/1‑149 (MQ=255) gcgAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 2:239726/1‑149 (MQ=255) gAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt > 3:66551/1‑149 (MQ=255) gAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt < 2:322277/149‑1 (MQ=255) aTGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa < 2:75095/149‑1 (MQ=255) gC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCtttt < 4:245337/64‑1 (MQ=255) gC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCtttt > 3:245337/1‑64 (MQ=255) caccacCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCg > 6:88640/1‑149 (MQ=255) | GCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGG‑ATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGA‑GCGGGTTTACCCGCGAATGC‑GAGGGTG‑‑ACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCG > NC_002947/4339952‑4340230 |
| 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: ‑ |