Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,296,329 | 0 | A | C | 33.9% | 111.5 / 25.3 | 59 | intergenic (‑149/‑66) | inaA/yfaH | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (36/3); new base C (0/20); total (36/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.34e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGT > REL606/2296254‑2296401 | gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca > 1:1434458/1‑76 (MQ=255) cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc < 1:1726776/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca < 2:201621/76‑1 (MQ=255) gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca > 1:938856/1‑76 (MQ=255) ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc > 2:1314874/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca > 2:17246/1‑76 (MQ=255) tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca < 1:235615/76‑1 (MQ=255) aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact > 2:1730656/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 1:1407134/1‑76 (MQ=255) gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac > 1:1407145/1‑76 (MQ=255) tttttttCCCACTCCTAGCGTGCCCCCGTTTTTCCCTTTCATTTGCAGCGCATAGATCCCCCCAGccactaccact < 1:842961/73‑1 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1687556/1‑76 (MQ=255) ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact > 1:1400693/1‑76 (MQ=255) cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa > 1:946350/1‑76 (MQ=255) tttGCCACTACTAGCGTGCCCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATa < 2:1243679/76‑1 (MQ=255) tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCCTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa > 1:994263/1‑76 (MQ=255) tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa > 2:1357913/1‑76 (MQ=255) tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa > 1:518263/1‑76 (MQ=255) ttcccACTACTAGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAg < 2:1139672/73‑1 (MQ=255) ccACTACTAGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgc < 1:245770/76‑1 (MQ=255) ctactaGGGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGcc < 1:1264036/76‑1 (MQ=255) ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 1:116075/1‑76 (MQ=255) ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc > 2:525764/1‑76 (MQ=255) tcctaGCGTCACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc < 1:348658/74‑1 (MQ=255) tactaGCGTGACCCCGCTTTTACCCTTCATTGGCAGCGCAAAAATCCCCCCAGCCACTACCACTAATAGTGCGccc < 1:1288494/76‑1 (MQ=255) tactaGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc < 1:1606601/76‑1 (MQ=255) ccaaGCGTGCCCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCa < 2:1178754/73‑1 (MQ=255) actgGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCa < 2:1449438/76‑1 (MQ=255) actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa > 1:1360670/1‑76 (MQ=255) ctaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAg > 2:806673/1‑76 (MQ=255) taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt > 2:411228/1‑76 (MQ=255) aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa > 1:1439308/1‑76 (MQ=255) gCGTGACCCCCTTTTTCCCCTTCATTTGCAGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAg < 1:781543/76‑1 (MQ=255) gggACACCGTTTTTCCCCTTCATTTCCACCGCATATATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGAt < 2:773493/74‑1 (MQ=255) tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa > 1:1905518/1‑76 (MQ=255) gACCCCCCTTTTCCCCTTCATTGGCGGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa < 2:1543639/76‑1 (MQ=255) acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa > 2:795766/1‑76 (MQ=255) cccGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAt < 2:675960/75‑1 (MQ=255) cccGCTTTTCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAt < 2:1754331/75‑1 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:1519127/1‑76 (MQ=255) tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt > 1:1495908/1‑76 (MQ=255) ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGAAGATAAAATGTGAAtt > 1:831251/1‑76 (MQ=255) tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 2:1288893/74‑1 (MQ=255) tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt < 1:1317849/74‑1 (MQ=255) tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACAAATAGTGCGCCCAGTAGATAAAATGTGAAttt > 1:1735134/1‑76 (MQ=255) cccTTCAGTTGCAGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa < 2:497845/76‑1 (MQ=255) ttCATTTGCAGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 2:962728/76‑1 (MQ=255) ttCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa < 2:962795/76‑1 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 1:491514/1‑76 (MQ=255) tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg > 2:493370/1‑76 (MQ=255) naTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt > 1:1283133/2‑76 (MQ=255) cATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt > 1:85719/1‑76 (MQ=255) tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa > 1:158911/1‑76 (MQ=255) cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa > 1:1584951/1‑76 (MQ=255) gcgcATAGATCCCCCAAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa < 1:958644/76‑1 (MQ=255) cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg > 1:1220089/1‑76 (MQ=255) aTAGTTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa < 2:521504/76‑1 (MQ=255) aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa > 1:1693303/1‑76 (MQ=255) aTAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa < 1:1929141/76‑1 (MQ=255) tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt > 2:930411/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 1:1085392/1‑76 (MQ=255) aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg > 2:1317185/1‑76 (MQ=255) aGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg < 2:956013/76‑1 (MQ=255) gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa < 1:1652658/76‑1 (MQ=255) aTTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGaa > 2:1533494/1‑76 (MQ=255) ttCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAg > 1:1264003/1‑76 (MQ=255) tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt > 1:163489/1‑76 (MQ=255) | GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGT > REL606/2296254‑2296401 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |