Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,388,923 | 0 | A | C | 53.3% | 30.5 / 36.5 | 30 | G151G (GGT→GGG) | nuoJ | NADH‑quinone oxidoreductase subunit J |
Reads supporting (aligned to +/- strand): ref base A (9/5); new base C (0/16); total (9/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCGTTTTTCTTTTCGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTG > NZ_CP009273/2388834‑2389008 | tccGTTTTTCTTTTCGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGa > 1:542092/1‑90 (MQ=255) tttcGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcag > 1:129574/1‑90 (MQ=255) tcTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGcc < 2:479532/90‑1 (MQ=255) cGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTcc > 2:492074/1‑90 (MQ=255) tCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCa > 1:317182/1‑90 (MQ=255) aGCACTTCCCCCCCACGCTCTTCACGCCCGACGGGGAAGGCCACAACCAGCCCTGCGAGCAGCGGCAAAAAACCCAGTTCCACCGCCAGt < 1:286606/90‑1 (MQ=255) tccACCCCCACCCTCTTCCCGACCGACGTGGAAGCCCAAAACCAGCCCTGCGAGCAGCAAAAAAAAAGCCAGTCCCACCGCCAGTACGTa < 1:125561/88‑1 (MQ=255) tCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcagcag < 1:308917/59‑1 (MQ=255) tCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcagcag > 2:308917/1‑59 (MQ=255) tCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTaa > 2:935739/1‑90 (MQ=255) ccGCCCCCTCTTCACGCCCGACGTGGAAGCCCACACCCAGCCCGGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGcc < 1:121075/90‑1 (MQ=255) ccGCCCCCTCTTCACGCCCGACGGGGAAGGCCACACCCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGcc < 1:20119/90‑1 (MQ=255) cccGCTCTTCACGACCGACGTGGAAGGCCACAACCAGCCCGGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGa < 2:283909/88‑1 (MQ=255) ccctctCCACGCCCGACGGGAAAGCCCAAAACCAGCCCCGCGAGCAGCAGAATAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAAc < 1:872078/88‑1 (MQ=255) cGCTTTCCACGCCCGAGGGGAAAGGCCAAACCCAGCCCGGCGAGCAGCAGAATAGAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAAc < 1:567789/90‑1 (MQ=255) cGCTCTTCCCGCCCGAGGGGGAAGGCCACAACAAGCCCTGGGAGCAGCAGAATAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAAc < 2:542092/90‑1 (MQ=255) cGCTCTTCACGCCCGACGGGAAAGCCCACAACCAGCCCGGCGAGCAGCAGAAAAAAAGCCAGTCCCACCGCCAGTACGTAAGGCCCGAAc < 1:551200/90‑1 (MQ=255) gacGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCAc > 1:545706/1‑90 (MQ=255) ccggggAAGCCCAAAACCAGCCCGGCGAGCAGCAGCAAAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt < 1:323207/86‑1 (MQ=255) acGTGGAAGGCCACAACCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt < 2:210402/90‑1 (MQ=255) acGTGGAAGCCCACACCCAGCCCGGCGAGCAGCAGCAAAAAACCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt < 1:268666/90‑1 (MQ=255) acGTGGAAGCCCACAACCAGCCCGGCGAGCAGCAGCAAAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt < 1:104388/90‑1 (MQ=255) ccACAACCAGCCCTGCGGGCAGCAGCAAAAAACCCATTCCCCCCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc < 1:802730/90‑1 (MQ=255) ccACAACCAGCCCTGCGAGCAGCAGAAAAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc < 1:397185/90‑1 (MQ=255) ccACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc > 1:381149/1‑90 (MQ=255) ccACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc > 2:62725/1‑90 (MQ=255) ccACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc > 2:813825/1‑90 (MQ=255) cccAGCCCGGCGAGCAGAAGAAAAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGAtt < 1:642568/89‑1 (MQ=255) cccAGACCTGCGAGCAGCAGAAAAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGAtt < 1:558551/89‑1 (MQ=255) aCCAGACCTGCGAGCAGCAGAATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGAtt < 2:343779/90‑1 (MQ=255) ccAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTg < 2:102460/90‑1 (MQ=255) | TCCGTTTTTCTTTTCGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTG > NZ_CP009273/2388834‑2389008 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |