Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,326,718 | 0 | T | A | 42.9% | 30.3 / 12.2 | 28 | Y44F (TAT→TTT) | adiC | arginine/agmatine antiporter |
Reads supporting (aligned to +/- strand): ref base T (12/4); new base A (0/12); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.62e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGG > NZ_CP009273/4326629‑4326802 | aagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat > 2:676075/1‑90 (MQ=255) ccaccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATg > 1:544262/1‑90 (MQ=255) accaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGc > 2:609829/1‑90 (MQ=255) ggACTTGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCAAGTGCACCGATAATCGTCACCAACCATCCAAAAATGGCAATc < 1:40015/90‑1 (MQ=255) ttGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATgg < 2:443957/81‑1 (MQ=255) ttGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATgg > 1:443957/1‑81 (MQ=255) tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGcc > 2:704390/1‑90 (MQ=255) aaCGACTTTTTGGCGTATCCCACGGAGACCCCCAGTCCCCCGATAATCGCCCCCACCCACCCATAAATGGCAATCCCGCCAGTAGAGGcc < 2:178136/90‑1 (MQ=255) ttgggCGTACCCCAGCGAGCCCCCCAGTGCCCCGAAAATCGCCCCCACCCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGc < 2:456800/87‑1 (MQ=255) gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAaca > 1:150712/1‑90 (MQ=255) gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc > 1:325017/1‑90 (MQ=255) gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc > 2:106565/1‑90 (MQ=255) gCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg > 1:240847/1‑90 (MQ=255) ccccaAGTGCACCGAAAACCTCCACCAACCACAAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg < 2:439995/89‑1 (MQ=255) ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa > 2:57885/1‑90 (MQ=255) cccAGTGCCCCGATAACCGCCACCACCCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc < 1:488981/90‑1 (MQ=255) cccAGTGCCCCGAAAACCGTCCCCACCCATCCAAAAATGCCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc < 2:447483/90‑1 (MQ=255) cccAGTGCACCGATAATCGTCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc < 1:57027/90‑1 (MQ=255) cccAGTGCACCGATAATCGTCACCAACCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc < 1:148052/90‑1 (MQ=255) aaGTGCACCGATAATCGTCACCAACCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAccc < 1:422727/89‑1 (MQ=255) cccACCGATAACCGCCCCCACCCACCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCat < 2:693382/88‑1 (MQ=255) aCCGATAATCGTCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat < 2:478397/90‑1 (MQ=255) aaTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTccc > 1:70744/1‑90 (MQ=255) gccaccacCCATCAAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac < 1:212012/88‑1 (MQ=255) tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc > 2:520967/1‑71 (MQ=255) tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc < 1:520967/71‑1 (MQ=255) caccacCCACCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac < 2:240847/90‑1 (MQ=255) cACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac > 1:569606/1‑90 (MQ=255) cACCAACCATCCATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac < 1:676075/90‑1 (MQ=255) ccatcaaaAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg < 2:48745/90‑1 (MQ=255) atccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg > 2:683091/1‑90 (MQ=255) | AAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGG > NZ_CP009273/4326629‑4326802 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |