| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 276,021 | 0 | A | G | 40.7% | 31.4 / 64.0 | 54 | pseudogene (53/249 nt) | BW25113_RS01360 | transposase |
| Reads supporting (aligned to +/- strand): ref base A (18/14); new base G (13/9); total (31/23) | |||||||||||
| 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.65e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTC > NZ_CP009273/275931‑276106 | ggACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATcaacaa < 1:419600/91‑1 (MQ=255) aCCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAAcc < 2:102466/91‑1 (MQ=255) aCCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAAcc < 2:590000/91‑1 (MQ=255) gTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGt > 2:802241/1‑91 (MQ=255) gTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGt > 2:317498/1‑91 (MQ=255) gTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGt > 2:230043/1‑91 (MQ=255) gTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGt > 1:420192/1‑91 (MQ=255) tACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGGGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGtt > 1:3571/1‑91 (MQ=255) ttATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGc > 2:461477/1‑91 (MQ=255) aTCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCtt > 2:36927/1‑91 (MQ=255) tCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTc < 2:252339/91‑1 (MQ=255) tCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTc < 1:605639/91‑1 (MQ=255) cAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAg > 2:611040/1‑91 (MQ=255) cAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAg < 1:461477/91‑1 (MQ=255) cAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAg < 1:230043/91‑1 (MQ=255) aGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGt > 2:419991/1‑91 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGtt > 1:274194/1‑91 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGtt > 2:64321/1‑91 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGtt > 1:242378/1‑91 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGa < 1:225717/80‑1 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGa > 2:225717/1‑80 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGtt > 2:398887/1‑91 (MQ=255) ggCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGtt > 1:137575/1‑91 (MQ=255) aaaGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTaaa > 2:662377/1‑91 (MQ=255) aaaGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTaaa > 1:809557/1‑91 (MQ=255) aaaGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTaaa > 1:609877/1‑91 (MQ=255) cTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCtt > 1:692135/1‑91 (MQ=255) cTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCtt > 2:547611/1‑91 (MQ=255) cTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCtt > 1:526629/1‑91 (MQ=255) tgtgGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTc < 2:811435/91‑1 (MQ=255) tgtgGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTc > 2:739300/1‑91 (MQ=255) tgtgGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTc < 1:250136/91‑1 (MQ=255) tggtgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATa < 2:764161/91‑1 (MQ=255) tggtgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATa < 1:419991/91‑1 (MQ=255) atatCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCt > 1:627337/1‑91 (MQ=255) tatCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCtt < 1:26106/91‑1 (MQ=255) cGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTAcc < 1:406170/91‑1 (MQ=255) cGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTAcc < 2:137575/91‑1 (MQ=255) agaCGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGc < 2:800418/91‑1 (MQ=255) cGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCa > 1:203458/1‑91 (MQ=255) gtgtTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCAt > 1:368464/1‑91 (MQ=255) tgtTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATg > 2:778956/1‑91 (MQ=255) tgtTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATg > 1:578295/1‑91 (MQ=255) tGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGc < 1:784728/91‑1 (MQ=255) cGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCAccgc < 1:547611/91‑3 (MQ=255) cGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACctc < 2:535575/91‑1 (MQ=255) gCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACctct > 1:798036/1‑91 (MQ=255) aCAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAg < 1:36927/91‑1 (MQ=255) cAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGc < 2:756944/91‑1 (MQ=255) cAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGc < 1:133814/91‑1 (MQ=255) tAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCAttt > 2:398850/1‑91 (MQ=255) tAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCAttt > 2:280474/1‑91 (MQ=255) gagaTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAAt < 2:789294/91‑1 (MQ=255) gagaTCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAAt < 2:573301/91‑1 (MQ=255) gaTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGa > 1:593979/1‑91 (MQ=255) caacaaCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTc < 1:354399/91‑1 (MQ=255) | GGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTC > NZ_CP009273/275931‑276106 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |