Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,359,867 | 0 | G | A | 40.9% | 66.0 / 105.3 | 88 | noncoding (1228/2905 nt) | rrlA | 23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base G (26/26); new base A (30/6); total (56/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.59e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC > minE/2359816‑2359935 | ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca > 1:1907301/1‑70 (MQ=255) tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:2361158/1‑71 (MQ=255) tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:2280583/1‑71 (MQ=255) tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:3135465/1‑71 (MQ=255) tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:2968886/1‑71 (MQ=255) ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc < 1:1967007/70‑1 (MQ=35) ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:3150403/1‑70 (MQ=255) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt > 1:1939222/1‑70 (MQ=35) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt > 1:2470773/1‑70 (MQ=35) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt > 1:2733858/1‑70 (MQ=35) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt > 1:3061162/1‑70 (MQ=255) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt > 1:3037095/1‑70 (MQ=255) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt > 1:1749878/1‑70 (MQ=255) tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCACCATACCTCACAGTACACCTTCACAGGCt > 1:1969175/1‑70 (MQ=11) aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt < 1:332117/70‑1 (MQ=35) aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa > 1:151163/1‑71 (MQ=34) aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt < 1:1649295/70‑1 (MQ=255) tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa < 1:794686/70‑1 (MQ=35) tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa < 1:3181243/70‑1 (MQ=35) gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg < 1:539391/71‑1 (MQ=34) tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc > 1:2370869/1‑67 (MQ=35) tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc > 1:420282/1‑67 (MQ=35) tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa > 1:335653/1‑71 (MQ=34) tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa > 1:1868856/1‑71 (MQ=34) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:2298462/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:2237071/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:1121859/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:1931644/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:1399411/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa < 1:1204994/70‑1 (MQ=35) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa < 1:1242278/70‑1 (MQ=255) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa < 1:383993/70‑1 (MQ=255) aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa < 1:428843/70‑1 (MQ=255) cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg < 1:1911558/71‑1 (MQ=34) cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg < 1:1684853/71‑1 (MQ=34) aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt < 1:1304075/70‑1 (MQ=35) aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc > 1:145065/1‑71 (MQ=34) aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc > 1:147692/1‑71 (MQ=34) aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc > 1:57369/1‑71 (MQ=255) tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc < 1:691510/71‑1 (MQ=34) tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGATc < 1:1877005/70‑1 (MQ=25) cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc > 1:758414/1‑67 (MQ=35) aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt < 1:714709/55‑1 (MQ=35) aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt < 1:3133983/55‑1 (MQ=35) aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt < 1:1618436/70‑1 (MQ=255) aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt < 1:605396/70‑1 (MQ=255) gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg > 1:1570900/1‑51 (MQ=11) cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:3132319/1‑46 (MQ=32) cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc > 1:2322363/1‑55 (MQ=35) cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc > 1:941583/1‑71 (MQ=34) cATTCGCACTTCTGATACCTCCAGCATGCCT‑ACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc > 1:1792660/1‑70 (MQ=255) cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc > 1:1694804/1‑71 (MQ=255) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:3088941/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:1371566/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:1437037/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:1360791/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:2842535/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca > 1:200207/1‑45 (MQ=32) aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa > 1:3080901/1‑58 (MQ=11) aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa > 1:731266/1‑58 (MQ=11) aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa > 1:2928787/1‑58 (MQ=11) aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa > 1:2463169/1‑58 (MQ=11) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa > 1:1925760/1‑71 (MQ=255) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa > 1:1237256/1‑71 (MQ=255) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa > 1:822661/1‑71 (MQ=255) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa > 1:23665/1‑71 (MQ=255) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca > 1:2073456/1‑70 (MQ=255) tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca > 1:853895/1‑70 (MQ=255) cGCACTTCTGATACCTCCAGCATACCTCACAGTACCCCTTCACAGGc > 1:555927/1‑47 (MQ=11) cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:1869252/1‑47 (MQ=11) cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc > 1:839819/1‑47 (MQ=11) cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa > 1:947785/1‑70 (MQ=255) gCACTTCTGATACCTCGAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac < 1:1265470/70‑1 (MQ=255) gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac < 1:2240640/70‑1 (MQ=35) gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac > 1:1080682/1‑70 (MQ=255) gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac > 1:2137436/1‑70 (MQ=255) gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac > 1:3401/1‑70 (MQ=255) cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa > 1:256485/1‑71 (MQ=34) cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa < 1:2844971/71‑1 (MQ=34) aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac < 1:924353/68‑1 (MQ=35) ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca > 1:985126/1‑70 (MQ=37) cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac < 1:1832262/69‑1 (MQ=39) aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt > 1:2876018/1‑71 (MQ=255) cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg < 1:738200/71‑1 (MQ=255) gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc < 1:2290805/71‑1 (MQ=255) aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg > 1:1210600/1‑70 (MQ=255) aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg > 1:2734822/1‑70 (MQ=255) aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc < 1:2549879/71‑1 (MQ=255) gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc < 1:2171599/38‑1 (MQ=40) gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc < 1:3169475/38‑1 (MQ=40) | CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC > minE/2359816‑2359935 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |