| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,287,044 | 0 | G | A | 27.3% | 67.7 / 23.7 | 45 | intergenic (+518/+22) | narI/rttR | nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript |
| Reads supporting (aligned to +/- strand): ref base G (20/12); new base A (4/8); total (24/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGC > NC_000913/1286911‑1287169 | aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAg > 1:64997/1‑139 (MQ=255) gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAgg < 2:483071/139‑1 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAggg > 2:272466/1‑139 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAggg < 1:282829/139‑1 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAggg > 1:523093/1‑139 (MQ=255) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCggg < 2:976501/139‑1 (MQ=11) cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCggg > 1:205381/1‑139 (MQ=11) gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCgggg < 2:205381/139‑2 (MQ=11) tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACt > 1:618189/1‑139 (MQ=255) cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGA < 1:584753‑M1/139‑1 (MQ=255) ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGAT < 1:937948‑M1/139‑1 (MQ=255) tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGT > 2:1037297‑M1/1‑139 (MQ=255) gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGG < 1:272466‑M1/139‑1 (MQ=255) gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGG > 1:22372‑M1/1‑139 (MQ=255) gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGT < 1:91629‑M1/135‑1 (MQ=255) gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGT > 2:91629‑M1/1‑135 (MQ=255) tGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGA < 1:613282‑M1/139‑1 (MQ=255) tGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGA < 1:212829‑M1/139‑1 (MQ=255) cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGAGGAAGGATTACt < 2:649562/139‑1 (MQ=11) gCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcg < 2:789750/139‑1 (MQ=11) ctggcGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc > 2:868377/5‑139 (MQ=11) cGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGA‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg < 1:462840/139‑1 (MQ=14) gAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggca > 1:76093‑M1/1‑114 (MQ=255) ccGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagat < 2:618189‑M1/139‑29 (MQ=255) atcGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctg < 1:1037297‑M1/138‑39 (MQ=255) aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGA‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg > 1:194544/1‑139 (MQ=14) aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccc > 1:862137‑M1/1‑96 (MQ=255) gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGA‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCa > 2:745980/1‑139 (MQ=11) tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggc < 2:523093‑M1/139‑50 (MQ=255) ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAgg > 1:699700/1‑48 (MQ=255) ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAgg < 2:699700/48‑1 (MQ=255) cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcggga < 2:413559‑M1/139‑59 (MQ=255) cGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGA‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc > 2:282829/1‑139 (MQ=14) cAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgc < 2:194544/139‑2 (MQ=11) atgtttgatTGCGGATTCGCATGAG‑ATTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggt < 2:309033‑M1/136‑73 (MQ=255) tGCGGATTCGTTGGGA‑AGTTCAGAGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt > 2:706/1‑121 (MQ=14) tGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattgaaat > 1:289876‑M1/1‑58 (MQ=255) tGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaat > 1:92735‑M1/1‑58 (MQ=255) tGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaac > 1:248177‑M1/1‑58 (MQ=255) tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCttaatcg > 1:388712/1‑134 (MQ=11) tGCGGATTCG‑TTGGGAAGTTCAGAGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt < 1:706/121‑1 (MQ=14) cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcaggtttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaattt > 1:64734‑M1/1‑56 (MQ=255) cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaattt > 1:905387‑M1/1‑56 (MQ=255) cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaattt > 1:259126‑M1/1‑56 (MQ=255) cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggtaattcaaattt > 2:243065‑M1/1‑56 (MQ=255) gattcgttgggaagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCACCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGc > 2:937948/13‑139 (MQ=11) | AGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGC > NC_000913/1286911‑1287169 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |