| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 578,811 | 0 | T | G | 50.0% | ‑1.1 / 13.3 | 18 | Q28P (CAG→CCG) | borD | prophage lipoprotein BorD |
| Reads supporting (aligned to +/- strand): ref base T (6/3); new base G (5/4); total (11/7) | |||||||||||
| 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.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACT > NC_000913/578570‑579042 | gCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACCCGAGCTTCCGGCGCAGGCTAAATGCATAACGTACTAAAACCGAGCAATCCATTTCCGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGtttt > 1:44892/1‑251 (MQ=255) gtGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATGCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCggtttg < 1:146489/211‑1 (MQ=255) ctgctaGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTTTTTTAACAACACTTTCTGCGCCTCCACAAATTTTAGCTGCATCGACAGTTGTCTGCTGTTCAATTCCAGAAACGAAGAAATGATGGTTGATGTTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTACATGTCTGTTGAg < 1:181275/214‑1 (MQ=255) ccAGCGTAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATGCCAGAAACGAAGAAATGATGGGTTATGGGTTCCTTTTGTGCTACTGCAGTCTGTTTGTTTAGAGCAGTAAATGTCTGTTTTGCACATTCTGTAATAAGCAGGGCCAGCGCAGTAGCg < 1:127967/251‑1 (MQ=255) gCCTAAAGTAATAAAACCGAGCAATCCATTTCCGAATGTTTGCTGGGTTTATGTTTTCCCAACATTTTCTGCGCCGCTACAAATTTTCGCTGCACAGCCAGTTTTCTTCTGCCAAATTCCAGAAACGAAGAAATGATGGGTGCTGGCTTCCTTTGGTGATACTGCTGTCTGTTTGTTTTGATCAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAg > 1:223692/1‑240 (MQ=255) aacaTTTTNTNNNCNNCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCNGAAACGAAGAAATGANGGNTGATNGTTNCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTg > 2:180964/1‑138 (MQ=255) aacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAATAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTg < 1:180964/138‑1 (MQ=255) caTTTTCTGCGCCGCTAAAAATTTTAGCTGCATCGACAGTTGTCTTCTGCCCAATTCCAGAATCGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTTGGTTTGTTTTGAACAGTAAATGGCTGTGGAGCACATCCTGGAATAAGCAGGTCCAGCGTAGTATCGAGTAGCATTGGTTTCATGTTGTTAGGCCCGATGCTTTTTGAAGTTCGCAGAATCGTAGTTGTAGAAAATTaaacaaac < 1:246363/251‑1 (MQ=255) gccgccACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACNg < 1:199179/113‑1 (MQ=255) gccgccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaacaa > 1:146081/1‑239 (MQ=255) caAATTTTGGCTGCATCGACANTTTTCTTCTNNCCNATTCCNGNNNCNNNGANATGATGGGTGATGGTTTCCNTTNNTGCTACNGCTGCCGGTTTGTTTTGAACAGTaaa < 2:60303/110‑1 (MQ=255) caAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTaaa > 1:60303/1‑110 (MQ=255) gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAAc > 1:246183/1‑95 (MQ=255) gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGGTTCCTTCGGTGCTACTGCTGCCGGTTTGTTTTg > 1:237842/1‑92 (MQ=255) cTGCATCTACATTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCtgtttgtgt < 1:193078/89‑3 (MQ=255) gACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTCTGGTGCTACTGCTGTCTGTTTGTTTTGCACAGTAAATGTATGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTAtg > 1:128561/1‑198 (MQ=255) aagaaATGATGGGTGATGGTTTCCTTTGGGGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATNCCCNNNNNtt > 2:157813/1‑144 (MQ=255) cTTTGGTGCTACTGCTGCCGGTTTGGTTTGAACAGTAAAAGCCTGTTGAGCGCATCCTGTACTAAGCCGGGCCGGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGGTACTACCGATGCTTTTTGAAGTTCGCCGCATCGGATGTCGAGCAAATTAACCAAGCCCTACACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTGTGTCAGGTGCGATTTATCGTCATTGTATTCCCGGATTAACt > 1:132394/1‑251 (MQ=255) gTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATg < 1:1606/219‑1 (MQ=255) ntACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCtttt > 1:76989/2‑115 (MQ=255) | GCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACT > NC_000913/578570‑579042 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |