| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 273,581 | 0 | T | G | 36.2% | 96.8 / 11.1 | 58 | G136G (GGT→GGG) | yagU | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base T (12/25); new base G (21/0); total (33/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC > REL606/273510‑273655 | tttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt < 2:1372958/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 2:987530/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 2:739137/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 2:1640863/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 2:132246/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 1:897001/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCttt < 1:352007/76‑1 (MQ=255) ttGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGttt > 1:371062/1‑76 (MQ=255) tGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGcc > 1:421541/1‑76 (MQ=255) tCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa < 1:1095425/76‑1 (MQ=255) cGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGcccca > 1:817268/1‑74 (MQ=255) gTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACt > 2:1177658/1‑76 (MQ=255) tATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACttt > 1:161138/1‑76 (MQ=255) aTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtttt < 1:100971/76‑1 (MQ=255) ttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACttttt > 1:905201/1‑76 (MQ=255) tgtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTg < 2:976757/76‑1 (MQ=255) gtgtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGTCAGGGCTTCCTGGCAGGTGCTTTAGCCCAACTTTTTGt < 2:1482419/76‑1 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACttttttttt > 1:390298/1‑75 (MQ=255) gtgGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTc > 1:40181/1‑76 (MQ=255) ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCAACTTTTTGTTCat > 2:55853/1‑76 (MQ=255) ggTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTat > 1:1560942/1‑76 (MQ=255) cTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCCACTTTTTGTTCATATGAt > 1:121168/1‑76 (MQ=255) cTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAt > 1:1650350/1‑76 (MQ=255) tGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtt < 2:975793/76‑1 (MQ=255) tGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGGTTTTGCCCCACTTTTTGTTTATATGAtt > 1:448769/1‑76 (MQ=255) gAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGAttt > 2:690520/1‑76 (MQ=255) aaGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttt > 2:363814/1‑75 (MQ=255) gTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttctt < 2:1048031/76‑1 (MQ=255) gTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtttctt < 2:1081462/76‑1 (MQ=255) aTTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGAtttctttc > 2:440707/1‑76 (MQ=255) tttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATTAtttcttttc > 2:262454/1‑74 (MQ=255) ttCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATCATTTTTTTtcc > 2:1342598/1‑76 (MQ=255) tCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTTATATGATTTTTTTCCct > 1:580986/1‑76 (MQ=255) aaaaaTTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTTATATGATTTTTTTCCCtctc > 2:1345588/1‑76 (MQ=255) tAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGa < 1:172666/76‑1 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTGTTTATATGATTTCTTTTCCTCTCATGGGAc > 2:314620/1‑76 (MQ=255) aaaCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGAc > 2:763871/1‑76 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACt < 2:176474/76‑1 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACt < 2:1014479/76‑1 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGGGCTTTCGCCCCACTTTTTTTTTATATTGTTTCTTTCCCTCTCATTGGaac > 1:1429444/1‑74 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCCACTTTTTTTTCATATGATTTCTTTTCCTCTCATGGGaat > 1:1429457/1‑74 (MQ=255) aaCTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCCTCTCATGGGACt > 1:1429441/1‑76 (MQ=255) aCTCTGGCAGGGCTTACTGGCNGGTGCTTTAGCCCAACGTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTg < 2:270805/76‑1 (MQ=255) aCTCTGGCAGGGCTTACTGGCAGGGGCTTTTGCCCCACTTTTTTTTTATATGATTTCTTTTCCTCTCATGGGACTg > 1:722137/1‑76 (MQ=255) ctctGGCAGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGa > 2:1253814/1‑76 (MQ=255) ctGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTTTTTCCCTCTCATGGGACTGaca > 2:1522449/1‑76 (MQ=255) tGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGacac < 2:321054/76‑1 (MQ=255) tGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTTATATGATTTTTTTTCCTCTCATGGGGCTGacac > 1:850073/1‑76 (MQ=255) aGGGCTTACTGGCAGGGGCTTTTGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccacc > 2:1129440/1‑76 (MQ=255) aGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccacc > 2:1367492/1‑76 (MQ=255) aGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACaccacc > 2:588390/1‑76 (MQ=255) gggCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACct < 1:1050030/76‑1 (MQ=255) aCTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGATTTCTTTTCCTCTCATGGGACTGACACCACCTCTGttt > 2:1478029/1‑76 (MQ=255) tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa < 2:195654/76‑1 (MQ=255) tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGa > 1:806152/1‑76 (MQ=255) tGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCCCTCATGGGACTGACACCACCTCTGTTTGa > 2:288883/1‑76 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt < 2:124213/76‑1 (MQ=255) ggCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGAt < 2:1121622/76‑1 (MQ=255) gCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATc > 1:366686/1‑76 (MQ=255) cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt > 1:165103/1‑76 (MQ=255) cAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCATTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCt > 2:670688/1‑76 (MQ=255) aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt < 2:1601091/76‑1 (MQ=255) aGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCtt > 1:902298/1‑76 (MQ=255) gTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTcc < 2:1224663/76‑1 (MQ=255) | TTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTTCTTTCCCTCTCATGGGACTGACACCACCTCTGTTTGATCTTCC > REL606/273510‑273655 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |