Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 42.3% | 66.2 / 32.6 | 52 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (25/5); new base C (0/22); total (25/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.07e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGA > REL606/2296255‑2296398 | tCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccac < 2:175866/76‑1 (MQ=255) cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc < 2:548495/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca > 2:1326030/1‑76 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca > 2:1326023/1‑76 (MQ=255) tATTCATGACTTTCTTTTGCCACTCCTCGCGTGACACCGCTTTTACCCTTCCTTTGCAGCGCATAGCTTCCACCAg > 1:1734775/1‑76 (MQ=255) tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg < 1:936104/76‑1 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc < 2:1107442/76‑1 (MQ=255) aTGAATTTTTTTTGCCACTACTAGCGTGCCCCCTTTTTTCCCCTTCATTTGCAGCGCATAGTTCCCCCCAGccact < 2:1700729/76‑1 (MQ=255) aTGAATTTTTTTTGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCCGCCCAAAGATTCCCCCAGccact < 2:380286/76‑1 (MQ=255) tGAATTTCTTTTGCCCCTACTAGCGTGCCCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGccacta < 2:1347662/76‑1 (MQ=255) aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacc > 2:547404/1‑76 (MQ=255) tCTTTTGCCACTACTAGCGTGCCCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGccactaccacta < 1:169894/76‑1 (MQ=255) tttGCCCCTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATa < 1:103818/76‑1 (MQ=255) tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa > 1:1134281/1‑76 (MQ=255) tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa > 1:1825589/1‑76 (MQ=255) tGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGt > 2:1004947/1‑76 (MQ=255) gCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTg < 2:173388/76‑1 (MQ=255) ccACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgc < 2:186525/76‑1 (MQ=255) ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgc > 1:600919/1‑76 (MQ=255) tactaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc < 2:1854963/76‑1 (MQ=255) tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc > 1:443973/1‑76 (MQ=255) actaGCGTGCCCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCa < 1:530932/76‑1 (MQ=255) ctaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAg > 1:490260/1‑76 (MQ=255) cgaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAg < 2:101583/74‑1 (MQ=255) tggCGTCACCCCGTTTTTCCCCTTCATTGGCAGCGCAGAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGt < 2:103443/74‑1 (MQ=255) taGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGt < 2:1580085/76‑1 (MQ=255) taGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGt < 2:518975/76‑1 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 1:1730026/1‑76 (MQ=255) aGCGTGACCCCCCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa < 1:1807659/76‑1 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 2:1754187/1‑76 (MQ=255) gCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:641987/76‑1 (MQ=255) cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa > 1:1566543/1‑76 (MQ=255) gTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt < 2:145490/76‑1 (MQ=255) gACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa < 2:155630/76‑1 (MQ=255) cccccGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa < 1:504317/73‑1 (MQ=255) accccGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa < 1:1545281/73‑1 (MQ=255) cTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa > 1:1244456/1‑76 (MQ=255) ttttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa > 1:613885/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 2:634626/1‑76 (MQ=255) ttccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt < 1:1291093/73‑1 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt > 1:629617/1‑76 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt > 2:450754/1‑76 (MQ=255) aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa > 1:1417512/1‑76 (MQ=255) ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt > 1:1765776/1‑76 (MQ=255) cTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc > 1:422243/1‑76 (MQ=255) cTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 2:1190605/76‑1 (MQ=255) cTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc < 2:1689442/76‑1 (MQ=255) ttCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 1:242064/76‑1 (MQ=255) ttCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa > 1:1850141/1‑76 (MQ=255) ttCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 1:916205/76‑1 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 1:213283/1‑76 (MQ=255) tCATTTGCAGCGCATAGATTACACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 2:711814/1‑76 (MQ=255) tCATTTGCAGCGCATAGATCNNCCCAGNCACNACCACTNATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg < 1:1379539/76‑1 (MQ=255) cttttGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt < 2:415382/74‑1 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 1:445530/1‑76 (MQ=255) gCAGCGAAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt < 1:790219/76‑1 (MQ=255) cAGCGCAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa < 1:625622/76‑1 (MQ=255) aGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATAc > 2:1391390/1‑76 (MQ=255) gcgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa < 2:140883/76‑1 (MQ=255) gcgcATAAATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa < 2:696051/76‑1 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:1355924/1‑76 (MQ=255) gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa > 1:507226/1‑76 (MQ=255) | TCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGA > REL606/2296255‑2296398 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |