| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,254,480 | 0 | A | C | 53.1% | 17.4 / 35.7 | 32 | V137G (GTG→GGG) | fruA | PTS fructose transporter subunit IIBC |
| Reads supporting (aligned to +/- strand): ref base A (8/7); new base C (0/17); total (8/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.12e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTC > NZ_CP009273/2254392‑2254563 | tGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTcac > 1:942960/1‑90 (MQ=255) gccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTcacccac < 1:1046852/90‑1 (MQ=255) ccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTcacccacc < 1:504560/90‑1 (MQ=255) aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCa < 2:996873/90‑1 (MQ=255) cAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGttt > 2:1170772/1‑90 (MQ=255) ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCg > 2:811150/1‑90 (MQ=255) ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCg > 1:272140/1‑90 (MQ=255) gtgtTTCCATTCCCCGCGCCAACAGACCCCCGGTTTTCAACTTCCCCCCCCCAGCCCCGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 2:1122931/90‑1 (MQ=255) gggATTGCATTCCCCGCGCCAACAGCACCCCGGGTTTAAACTTTCACCCCCCAGCCACTTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 1:1113302/88‑1 (MQ=255) gggATTGCATTACCCGCCCCAACAGAACCACGGGTTTCAACTTTCCCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 1:1170772/88‑1 (MQ=255) gTGATTGCATTCCCCGCGCCAACAGACCCACGGGTTTCAACTTTCCCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 1:105747/90‑1 (MQ=255) gTGATTGCATTCCCCGCGCCAACAGAACCACGGGTTTCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 1:1241846/90‑1 (MQ=255) gTGATTCCTTTCCCCGCGCCAACAGACCCAGGGGTTTCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 2:259220/90‑1 (MQ=255) gTGATTCCATTCCCCGCGCCAAAAGAACCACGGTTTTCAACTTTCCCCCCCCAGCCCCTTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc < 2:1290925/90‑1 (MQ=255) ttttGCTTTCCCCGCCCAAAAAAACCCAGGGTTTTCAATTTCCCCCCCCCAGCCCCTTTTTTTCGCTTCGTTTTCAATGGCTTCAgccgc < 2:263063/88‑1 (MQ=255) gATTGCATTCCCCGCGCCAACAGACCCACGGGTTTCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgccgc < 1:660341/90‑1 (MQ=255) gATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGt > 1:1209291/1‑72 (MQ=255) gATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGt < 2:1209291/72‑1 (MQ=255) aTTCCCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATaaa < 1:382578/90‑1 (MQ=255) aTTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATaaa < 1:14897/90‑1 (MQ=255) cccGCTCCACAAGACCCAGGGTTTTCATCTTCCCCCCCCCAGCCGCGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtg < 2:778815/90‑1 (MQ=255) cccGCGCCAACAGACCCCCGGGTTTCAACTTTCCCCCCCCAGCCCCGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtg < 1:847434/90‑1 (MQ=255) cccGCCCCAACAAACCCAGGGTTTTAAATTTCCCCCCACCAGCCCCTTTTTTCCGCTTCGTTTTCAATGGCTTCAGCCGCCATAAAGgtg < 2:32916/90‑1 (MQ=255) cgcCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg > 1:265293/1‑90 (MQ=255) cgcCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg > 1:934663/1‑90 (MQ=255) ccccGGTTTTCACCTTCCCCCCACCAGCCTTGTTTTTTCCCTTCGGTTTAAATGGCTCCAGCCGCCATAAAGGTGTGTGCTACGCCAGTc < 2:744139/87‑1 (MQ=255) ccACGGGTTTCACCTTCCCCCCCCCAGCCACTTTTTTTCGCTTCGTTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTc < 2:876491/90‑1 (MQ=255) ccACGGGTTTCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTc < 1:927817/90‑1 (MQ=255) aCGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCgg > 2:395073/1‑90 (MQ=255) gttttCAACTTCCCCCCCCCAGCCACTTTTTTTCGCTTCGGTTTCAAGGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCa < 2:298591/88‑1 (MQ=255) ttttCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCaa < 1:811150/89‑1 (MQ=255) tttCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAg > 2:228121/1‑90 (MQ=255) ttCAACTTTCACCCCCCAGCCACGTTTTTTCCCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGc < 1:1257214/90‑1 (MQ=255) tCACCTTCCCCCCCCCACCCACTTTTTTCCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCa < 2:712710/90‑1 (MQ=255) aCTTCCACCCCCCAGCCACGTTTTTTCGCTTCGTTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTc < 1:1012972/90‑1 (MQ=255) | TGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTC > NZ_CP009273/2254392‑2254563 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |