| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 516,893 | 0 | G | C | 30.3% | 53.3 / 16.1 | 34 | G104A (GGA→GCA) | fetB | iron export ABC transporter permease; peroxide resistance protein |
| Reads supporting (aligned to +/- strand): ref base G (16/7); new base C (0/10); total (16/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.97e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCT > NC_000913/516830‑516957 | aaCGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAAtt > 2:819172/1‑68 (MQ=255) gCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTAcc > 2:728980/1‑68 (MQ=255) ttAAATATATTGCTAATGCTTTTACCTCACCGTTTATTGCCATTACGGCCGGGGCGGCAATTACCCTg < 1:699694/67‑1 (MQ=255) aaaTATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGc > 2:955448/1‑68 (MQ=255) tatatTGCTAACGATTTTATCTCATCTTTATTTGCCATAAGGGTCGGGGCGGTAATTACCCTGGCGGt < 2:982197/68‑1 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCTTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGt < 2:622492/68‑1 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGt > 2:714218/1‑68 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGt < 2:933499/68‑1 (MQ=255) aatatTGCTAAAGCTTTAAACTCATAGTTTATAGCCATTACGGTCGGTGCGGGAATTACACTGGCGGt > 2:375871/2‑68 (MQ=255) atatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTg < 2:974645/68‑1 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc > 2:562140/1‑68 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc > 2:595430/1‑68 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc > 1:352263/1‑68 (MQ=255) tGCTAAAGCTTTTATCTCATCGTTAATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGa < 2:914685/68‑1 (MQ=255) tGCAAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGa < 2:916161/68‑1 (MQ=255) gCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGAt > 2:433776/1‑68 (MQ=255) atgCTTTTATCTCATCGTTTATTGCCATTACGGGCGGGGCGTGAATTACCCTGGCGGTGCTGATtctc < 1:520502/66‑1 (MQ=255) aatcTTTTATCTCATCGTTATTTGCCTTTACGGGCGGGGCGGCAATTACCCTGGCTGTGCTGATtctc < 2:170704/65‑1 (MQ=255) aaGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATtctc < 2:969748/68‑1 (MQ=255) aaGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGATtctc < 2:931554/68‑1 (MQ=255) cTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCa > 2:267373/1‑68 (MQ=255) tctcATCTTTTATTGCCATAACGGGCGGGGCGGCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg < 2:556118/68‑1 (MQ=38) tctcATCGTTTATTGCCATTTCGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg > 2:989401/1‑68 (MQ=255) tctcATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg < 2:462637/68‑1 (MQ=255) cGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGaa < 1:595425/68‑1 (MQ=255) cGTTTATTGCCATTACGGTCGGGGCGGCATTTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGaa < 2:753236/68‑1 (MQ=255) gCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATccc > 1:605301/1‑68 (MQ=255) aTTCCGGCCGGGGCGGCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGAt < 1:433625/68‑1 (MQ=255) tACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGAGTCTCTCAGGGTCGATTGTATTTATCCCGATGc > 1:63795/1‑68 (MQ=255) ggTCGGGGCGGGAAATACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAAATTATCCCGCTGCAgg > 2:771378/1‑68 (MQ=255) gTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGt < 2:568607/68‑1 (MQ=255) tCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTg > 1:377200/1‑68 (MQ=255) cGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGGATTTATCCCGATGCAGGTGATCCCt > 1:221146/1‑68 (MQ=255) cGGCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCt < 2:379697/68‑1 (MQ=255) | AACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCT > NC_000913/516830‑516957 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |