| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,763,767 | 0 | T | . | 58.2% | 92.8 / 130.2 | 67 | noncoding (77/1542 nt) | rrsG | 16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base T (28/0); new base . (0/39); total (28/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.71e-19 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATC > minE/1763722‑1763817 | taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:663113/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:177187/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:190928/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:230091/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:318786/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:342943/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:350950/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:413511/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:422890/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:483907/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:503469/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:613234/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:613558/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:307597/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:694970/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:738556/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:743007/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:774262/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:800131/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:812711/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:816414/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:879245/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:884642/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:994722/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:9997/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1361268/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1218024/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1208455/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1198768/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1281600/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1283422/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1091134/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1415950/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1445006/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1035923/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1047380/58‑1 (MQ=2)taTTACTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1492900/58‑1 (MQ=2) cTCACCCGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1498990/54‑1 (MQ=2) cccGTTCGCCACTCGTCAGCG‑‑G‑AGCAAGCT‑C‑‑CCTGTTACCGTTCGACTTg < 1:1082571/50‑1 (MQ=2) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAt > 1:1232390/1‑54 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:931527/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1010567/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:668178/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1018204/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:706736/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:992780/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1140863/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:923589/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:589687/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1081173/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1062548/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:644289/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1276479/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:561341/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:522908/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:47642/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:452666/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1343788/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:358740/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:142291/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1005789/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:265246/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1449302/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:227485/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1504531/1‑55 (MQ=33) tgctTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATc > 1:1499520/1‑55 (MQ=33) tgctTCCTGCTACCGTTCGACTTGCATGTGTTAGGCCTgccgc > 1:149874/1‑43 (MQ=255) | CATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATC > minE/1763722‑1763817 |
| 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: ‑ |