| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,020,300 | 0 | A | C | 51.3% | 47.5 / 57.7 | 38 | E391D (GAA→GAC) | guaD | guanine deaminase |
| Reads supporting (aligned to +/- strand): ref base A (16/2); new base C (0/19); total (16/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATG > NZ_CP009273/3020214‑3020388 | aCGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAccc > 2:1336421/1‑90 (MQ=255) gCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCAc > 1:165645/1‑90 (MQ=255) gggAGCGAAATCTCGGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCAc < 2:452874/89‑1 (MQ=255) cggAGCAAAATCTCTGGGCTTTGACGTTTTGTTTGGCAATTTTTTCCCTGGAAAAGGGGTTGTTTTCGTGGTGATGGACCCCACCGCCAc < 1:33732/90‑1 (MQ=255) gAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc > 1:203969/1‑90 (MQ=255) gAAATCTCGGGGCTTTGACGTTTTGATTGGCACTTTTTTCCCGGGAAAAGGGGTTGTTTTGGGGGTGAGGGACCCCCCCGCCACTCCGCt < 1:785783/90‑1 (MQ=255) aTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAca > 2:840402/1‑90 (MQ=255) aTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAca > 2:43248/1‑90 (MQ=255) ctctGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTAcagc < 2:78412/90‑1 (MQ=255) tctGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagca > 2:340767/1‑90 (MQ=255) gggCCTTGACGTTTTGTTTGGCAATTTTTTCCCTGGCAAAGGGGTTGTTTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCt < 1:652202/90‑1 (MQ=255) gggCCTTGACGATTTGATTGGCAATTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCt < 2:794989/90‑1 (MQ=255) gggCCTTGACGATTTGATTGGCAACTTTTTCCCGGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCt < 1:585356/90‑1 (MQ=255) cTTGACGTTTTGATTGGCAATTTTTTCCCGGGAAAAGAGGCTGTTTTCGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTgcgc < 1:340767/90‑1 (MQ=255) ttGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt < 2:165645/90‑1 (MQ=255) tGATTGGCAATTTTTTCCCGGGCAAAGGGGCTGTTTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct < 2:1034931/90‑1 (MQ=255) tttCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt < 1:365851/90‑1 (MQ=255) tttCCCGGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt < 2:203969/90‑1 (MQ=255) tGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaa > 1:876622/1‑90 (MQ=255) ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa > 1:1312700/1‑90 (MQ=255) ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa > 2:1399924/1‑90 (MQ=255) ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACaaa > 2:1285939/1‑90 (MQ=255) ggCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAACAGCTGCGCTATGCCAACTCTGTGTCTTTAGTCGACaaa > 1:1057667/1‑90 (MQ=255) gaaaaGGGGTTGTTTTGGTGGTGAGGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt < 2:1312700/88‑1 (MQ=255) ggCTGTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt < 1:840402/90‑1 (MQ=255) ggCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt < 2:94555/90‑1 (MQ=255) gCTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg < 2:63491/90‑1 (MQ=255) cTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtga > 2:939397/1‑90 (MQ=255) tGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgat > 2:230014/1‑90 (MQ=255) ttCGTGGTACGGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg < 2:1415085/90‑1 (MQ=255) tggtggtgATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga < 1:529240/88‑1 (MQ=255) tCGTGGTGAGGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga < 1:654353/90‑1 (MQ=255) cGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc > 2:385404/1‑90 (MQ=255) ggtgATGGATCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt < 2:219845/90‑1 (MQ=255) ggtgATGCCCCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt < 2:1578346/90‑1 (MQ=255) ggtgATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt < 1:92616/90‑1 (MQ=255) tgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg > 2:498254/1‑90 (MQ=255) gAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGAt > 1:619352/1‑90 (MQ=255) aaCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATg > 2:164980/1‑90 (MQ=255) | ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATG > NZ_CP009273/3020214‑3020388 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |