| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,302,452 | 0 | A | C | 33.3% | 26.6 / 19.3 | 24 | L101R (CTG→CGG) | napA | nitrate reductase, periplasmic, large subunit |
| Reads supporting (aligned to +/- strand): ref base A (8/8); new base C (5/3); total (13/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.79e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAAC > NC_000913/2302309‑2302599 | cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGc > 2:2894/1‑149 (MQ=255) cGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGta < 1:2894/149‑2 (MQ=255) aCATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCaaa < 4:80395/149‑1 (MQ=255) ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg < 6:87443/139‑1 (MQ=255) ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg > 5:87443/1‑139 (MQ=255) aaTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg < 1:51729/134‑1 (MQ=255) aaTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg > 2:51729/1‑134 (MQ=255) cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAgg < 8:119550/121‑1 (MQ=255) cacaTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAgg > 7:119550/1‑121 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtat > 2:147127/1‑84 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtat < 1:147127/86‑3 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg < 6:23027/112‑1 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg > 5:23027/1‑112 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc < 6:123921/126‑1 (MQ=255) aTCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc > 5:123921/1‑126 (MQ=255) aaGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCtgctgc > 4:73221/1‑78 (MQ=255) aaGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCtgctgc < 3:73221/78‑1 (MQ=255) cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg < 4:71439/108‑1 (MQ=255) cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAg > 3:71439/1‑108 (MQ=255) gggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAAc > 2:33976/1‑127 (MQ=255) gggTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAAc < 1:33976/127‑1 (MQ=255) cttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGt > 2:114699/1‑149 (MQ=255) cGTTTTTCATACGCAGCCGCGGCTGCGTCAAACGGTCTTGACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCg > 8:69487/1‑149 (MQ=255) gaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcg > 6:142177/2‑149 (MQ=255) cagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAAc < 5:75058/149‑1 (MQ=255) | CAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAAC > NC_000913/2302309‑2302599 |
| 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: ‑ |