| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,740,467 | 0 | T | G | 28.0% | 51.0 / 32.6 | 50 | intergenic (+32/+39) | hutF/PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
| Reads supporting (aligned to +/- strand): ref base T (14/22); new base G (8/6); total (22/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.43e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG > NC_002947/5740320‑5740613 | gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAGGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑Gt < 2:359621/148‑1 (MQ=255) cgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GGGCTTGCAcgcg < 6:289094/149‑1 (MQ=255) cgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCACCCCTTCGC‑GTGCAAGccccc > 5:446936/1‑146 (MQ=255) gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCcgcg > 8:316170/1‑148 (MQ=255) aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGaa < 6:388714/149‑1 (MQ=255) aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGC‑GGGCTTGCACGCGaa < 5:293196/149‑1 (MQ=39) tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAg < 1:426851/149‑1 (MQ=39) aTGGTGAATGGGCAATGGGTGGTGAGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGACAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTGCAGACCATAGC‑GTGCAACCCAGCGAAGGGtt > 2:71529/1‑149 (MQ=255) ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCTTAGC‑GTGCCAACCCCCGAAGGGttt > 7:182112/1‑147 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCGTAGC‑GGGCAAGCCCCCGAAGGGTTgca > 6:355005/1‑147 (MQ=255) aaggggCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTgcaac < 3:290307/145‑1 (MQ=255) gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTAGCGGGGTTGCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAgggct < 8:320070/124‑3 (MQ=37) gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GGGCTTGCACGCGAAg > 6:163807/1‑120 (MQ=38) gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GGGCTTGCACGCGAAg < 5:163807/120‑1 (MQ=38) gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAgggct > 7:320070/1‑122 (MQ=37) cccggggCGGCATGCCCAGGAAGATGCACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTg < 4:278911/146‑1 (MQ=255) ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTCGTTTCAGCCCATAGCGGTGCAA‑ACCGCGCAAGGTTGGAACG‑CAACCCCAAAacacaacggg > 8:72047/1‑143 (MQ=255) ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTg > 6:446267/1‑149 (MQ=255) cccAGGAAGATTACAGCAGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGGTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTAttt < 6:305090/149‑1 (MQ=255) ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGtt < 7:209168/108‑1 (MQ=255) ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGtt > 8:209168/1‑108 (MQ=255) cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc < 6:335384/149‑1 (MQ=255) aaaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGAAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc < 5:355005/148‑1 (MQ=255) aGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGGGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGTCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgc < 1:315661/149‑1 (MQ=255) gggcgggcGTTTGCCATGGTGTTGCCCCAGCTTTTGGCCTCAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc < 8:182112/149‑1 (MQ=255) gggcgggcGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc < 7:316170/149‑1 (MQ=255) gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc < 5:206267/132‑1 (MQ=255) gTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGGAAGGCCGCGAAGGGTTTCAAAG‑CAACCCCAAAAAAATCAGTGGATTTTTAATATTTGCTTATAGGAAGCCCGCCAGAtg > 2:279760/1‑146 (MQ=255) aaGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCATAGC‑GTGGAAGCCCGCGAAGGGTTTGAACG‑CAAACCCAAAAAAAACAGTGGATTTTTTCTATTTTCTGATAGGACGCCCCCCAGATAAAGCGGCTg > 1:126686/1‑149 (MQ=255) gcgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGC‑GTGCAAGCCCCCGCACGGTTTCAACC‑CAACCCCAAAAAACTCAGTGGAATTTTAATATTTTGTGATAGGACGCCCGCCAGAAAAAGCGGCTGGG‑GTCAt > 2:39354/1‑148 (MQ=255) gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc < 5:185224/137‑1 (MQ=255) caaGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAACCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATAccc < 7:72047/147‑1 (MQ=255) cAGCTTTTGGGCTGTGCTTACCGGGTTTGCGTGGCAGCCGTTCGC‑GTGCAAGCCCGCGAAGGGTTGCATCT‑CAACCCCAAAAAACTCACTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATAcccc < 7:280009/149‑1 (MQ=255) aGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCACCCCATAGC‑CTGCAACCCCGCCAAGGGTTGCAACC‑CCACCCCAAAAAACAAACTGCAATTTTAAGAATTTCTGATCCGGCGCGCGCCAACGAAAGCGGCGGGG‑GTCATAacc > 8:309252/1‑147 (MQ=255) gggCTGAGCTCACCGGGGTGGCGTTCCAGCCCTTCGC‑GTGCAAGCCCGCTATGGGTTGAAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGcc < 5:32510/148‑1 (MQ=255) gcaCCGGGGTTGCGGTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGAAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGa < 8:312417/147‑1 (MQ=255) tACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTgcaacg‑caac < 1:148435/57‑1 (MQ=255) ccGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGCT‑GTCATACCCCTGCTGCCGCGCCCTGGAca < 1:39354/148‑1 (MQ=255) ggggTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGAcagca < 2:126686/149‑1 (MQ=255) gcgtttcATCCCTTAGC‑TTGCAATCCCGCGAATGGTTGCAACG‑CAACCCCAAAAAACTCACTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTCGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt < 1:71529/148‑1 (MQ=255) gttgcAGCCCTTCGC‑GTGCAAGCCCGCGAAGGTTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa < 1:279760/149‑1 (MQ=255) tgcAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg > 3:334147/1‑149 (MQ=255) caacccTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt > 7:185019/4‑149 (MQ=38) aacccTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCgcgc > 7:375557/3‑125 (MQ=37) aacccTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCgcgc < 8:375557/123‑1 (MQ=37) acccTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCctgc > 2:387795/2‑116 (MQ=37) cTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg > 1:27367/1‑149 (MQ=39) ttCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCct > 5:386956/1‑110 (MQ=255) ttCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCct < 6:386956/110‑1 (MQ=255) ttCGC‑GGGCTTGCACGCGAAGGGCTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 1:12787/1‑149 (MQ=39) c‑gTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg > 3:386944/1‑149 (MQ=255) | GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACG‑CAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGT‑GTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG > NC_002947/5740320‑5740613 |
| 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: ‑ |