Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | ACB143 | 5,386,507:1 | +C | 100% | intergenic (‑516/+1135) | CJ019_05070 ← / ← nptII | hypothetical protein/neomycin phosphotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 5,386,507 | 1 | . | C | 94.8% | 57.9 / ‑3.5 | 20 | intergenic (‑516/+1135) | CJ019_05070/nptII | hypothetical protein/neomycin phosphotransferase |
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (7/11); minor base A (0/1); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CG‑T‑CA‑‑GGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG > ACB143/5386410‑5386604 | ccTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGc < 5:214843/136‑1 (MQ=255) ttGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGaa < 8:6834/135‑1 (MQ=255) ttGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGaa < 2:246551/135‑1 (MQ=255) ttGATAACAAAAATTGATCGAATTATTTCATTAAAGACCCTAAATTGACTTCTCGGCGGGCTTGTCGACGACGGCGGTCTCCGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGaa < 5:191975/135‑1 (MQ=255) aaaTTCATCGAATTCTTTAATTTAACACCCTAAAATCACATCTCGGCCGGCTTGGCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc < 6:167200/136‑1 (MQ=255) accGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGTGGTG‑T‑CC‑‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg < 7:123231/133‑1 (MQ=255) aTTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc < 8:169747/136‑1 (MQ=255) atatGACATCTCGGGCTGCTTGTCGACGACGTCGGTCTCAGT‑‑GGTT‑TAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg < 8:230599/136‑1 (MQ=255) atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg < 3:157320/136‑1 (MQ=255) atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg < 3:291710/136‑1 (MQ=255) atatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg < 6:271704/136‑1 (MQ=255) atatGACATCTCGGCCGGCTTGCCTACGACGGCGGTCTCAGT‑‑GG‑T‑GAACGGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAgg < 6:269843/136‑1 (MQ=255) tatGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAggg > 6:23530/1‑136 (MQ=255) tatGACATCTCGGCCGGCTTGTCCACGACCGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCGAGCAACGCATGAGAAAGCCCCCTGAAGCACACCGTCCGGCGGCTGTTTTAAAGCCCTGCGGGTGCCAgcg > 6:88787/1‑134 (MQ=255) ctcGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt > 2:29701/1‑136 (MQ=255) ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc > 2:144580/1‑95 (MQ=255) ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTcc < 1:144580/95‑1 (MQ=255) ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt > 6:223733/1‑136 (MQ=255) ctcGGCCGGCTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCt > 1:244680/1‑136 (MQ=255) ggccggcTTGTCGACGACGGCGGTCTCAGT‑‑GG‑TGTAC‑GGTACAAACCCCGACGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTgg > 1:47700/1‑136 (MQ=255) | CCTTTTTTATTGATAACAAAAATTGATCGAATTCTTTCATTTAAGACCCTAATATGACATCTCGGCCGGCTTGTCGACGACGGCGGTCTCCGT‑‑CG‑T‑CA‑‑GGATC‑‑ATCCGGGCGCTAGCAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGCGGGTGCCAGGGCGTGCCCTTGG > ACB143/5386410‑5386604 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |