| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,351,234 | 0 | C | A | 21.4% | 58.4 / 12.5 | 28 | intergenic (‑195/+173) | fabI/ycjD | enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein |
| Reads supporting (aligned to +/- strand): ref base C (10/12); new base A (2/4); total (12/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGC > NC_000913/1351092‑1351334 | tcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggt > 5:187037/1‑148 (MQ=255) gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggtta < 6:187037/149‑4 (MQ=255) ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttag < 8:136453/149‑5 (MQ=255) ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttag < 6:71610/149‑5 (MQ=255) aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGAcc < 2:46765/127‑1 (MQ=255) ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttagg > 8:116401/1‑117 (MQ=255) ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggttagg < 7:116401/122‑6 (MQ=255) cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCAACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATa > 4:65586/1‑148 (MQ=255) tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg < 2:41460/92‑1 (MQ=255) tgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg > 1:41460/1‑92 (MQ=255) gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTc < 6:160609/137‑1 (MQ=255) gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTc > 5:160609/1‑137 (MQ=255) taGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCa > 2:161150/1‑149 (MQ=255) gTTGCCAGGTCCGACAGGAGCAGGCTGCGGGAGTGGGTGCGCGTTGCCCCTCACCCTAACCCGCTCCCCAGAGGGGCGAGGGTACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg < 3:166767/149‑1 (MQ=38) gACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaac > 2:150137/1‑108 (MQ=255) gACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaac < 1:150137/108‑1 (MQ=255) ccGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGAc > 2:62888/1‑149 (MQ=255) gAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGt < 1:161150/149‑1 (MQ=255) cGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg < 3:65586/149‑1 (MQ=255) ggCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc < 1:62888/149‑1 (MQ=255) gggCGCTTTTCCCCTCACCCT‑ACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATgg > 2:67318/1‑119 (MQ=255) gggCGCTTTTCCCCTCACCCT‑ACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATgg < 1:67318/119‑1 (MQ=255) tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGcc < 2:208340/89‑1 (MQ=255) tCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGcc > 1:208340/1‑89 (MQ=255) ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCaaa > 3:62926/1‑44 (MQ=255) ccTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCaaa < 4:62926/44‑1 (MQ=255) ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAg > 7:51962/1‑109 (MQ=255) ctctCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAg < 8:51962/109‑1 (MQ=255) | TCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGC > NC_000913/1351092‑1351334 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |