Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 45.7% | 38.0 / 17.6 | 35 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (1/18); new base G (16/0); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAC > REL606/3492010‑3492156 | cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg > 1:1179678/1‑75 (MQ=255) ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg < 1:707869/76‑1 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggggg > 1:856347/1‑73 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGtata > 2:1504966/1‑76 (MQ=255) cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc < 1:803633/76‑1 (MQ=255) cGTCGCACAAAACAAAGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc < 2:58669/76‑1 (MQ=255) gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGNGAGTCAGTGGCTGGTGGAGTATACgca < 1:1524054/76‑1 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACgcag > 1:851524/1‑76 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc < 2:732142/76‑1 (MQ=255) gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca < 1:685567/76‑1 (MQ=255) acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGa < 1:606119/76‑1 (MQ=255) caAAACAAGGCATGNGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGan < 2:219001/76‑2 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCCGTc > 1:1327568/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGATGTc > 2:1208348/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc > 2:1369288/1‑76 (MQ=255) aaGGCATGNGGCTAGATGCNGATCAACAATCTCGTGTGAGTCNGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTc < 2:1290082/76‑1 (MQ=255) gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt < 2:869118/76‑1 (MQ=255) gCATGTGGCTAAATGCGGATCAACAATCTCGTGNGAGTCAGTGGCTGGGGGAGGATTCGCAGCAGAAGCTGTCGCt > 2:1292421/1‑76 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg < 2:1178847/76‑1 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg < 1:926229/76‑1 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGcctg > 2:1101975/1‑73 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTgg > 1:1235348/1‑76 (MQ=255) tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCCGGCGCTggg > 2:67820/1‑75 (MQ=255) gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCCGAAGCTGTCGCTGGcg > 1:348481/1‑76 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtc < 1:596014/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCCGgcggcggg > 2:608474/1‑74 (MQ=255) tAAAGGAGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg > 1:89277/1‑76 (MQ=255) naaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCCGNAGCTGGCGGCGTgg > 2:218377/2‑75 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTgg > 2:1253654/1‑75 (MQ=255) gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATg < 1:1489398/76‑1 (MQ=255) tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa < 2:1480087/76‑1 (MQ=255) tacaaTCTCGTGTNAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg < 1:196945/75‑1 (MQ=255) acaaTCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGTCGGGGTTTAATGTAAgg > 1:435884/1‑76 (MQ=255) caaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTAATGTAAggg > 1:1078729/1‑75 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg > 2:1066099/1‑76 (MQ=255) aTCTCGTGTGATTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGTTTAATGTAAGgggg > 1:841840/1‑76 (MQ=255) cGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAg < 1:436080/76‑1 (MQ=255) tgATTCGGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCa < 2:298527/76‑1 (MQ=255) gTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTCGCGTCGTGATTTAATGTAAGNTGGCCAGGTCAGATTTCa < 1:1468212/76‑1 (MQ=255) gCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCGGGGCAGATTTTAATc > 1:664096/1‑76 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg < 1:376649/76‑1 (MQ=255) gtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAc < 2:413679/76‑1 (MQ=255) | CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAC > REL606/3492010‑3492156 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |