Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 341,944 | 0 | T | C | 56.1% | 2.7 / 83.6 | 57 | G35G (GGA→GGG) | yahN | LysE family transporter |
Reads supporting (aligned to +/- strand): ref base T (11/14); new base C (15/17); total (26/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAG > NZ_CP009273/341808‑342078 | gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGcccc > 1:47998/1‑139 (MQ=255) gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGcccc > 1:358770/1‑139 (MQ=255) gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTcc > 1:191535/1‑139 (MQ=255) cGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCcccc > 1:573675/1‑139 (MQ=255) gTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCg > 2:187567/1‑139 (MQ=255) gCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTaa > 1:466111/1‑138 (MQ=255) gTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTaaaa > 2:135672/1‑138 (MQ=255) gTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTaaaa > 2:483918/1‑138 (MQ=255) cAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAg < 2:533365/138‑1 (MQ=255) acaaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGaa < 1:568958/139‑1 (MQ=255) acaaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGaa < 1:483918/139‑1 (MQ=255) acaaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCAGGATTAAAAAAAGTAATCACGaa < 2:20180/139‑1 (MQ=255) caaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGCTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAAc < 2:235741/139‑1 (MQ=255) caacccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCt < 2:466111/139‑1 (MQ=255) cccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACg < 2:358770/139‑1 (MQ=255) cGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt < 1:537360/139‑1 (MQ=255) cGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt < 2:86398/139‑1 (MQ=255) cGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt < 2:47998/139‑1 (MQ=255) gAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAg < 2:117436/115‑1 (MQ=255) gAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAg > 1:117436/1‑115 (MQ=255) aTCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAg > 2:340895/1‑104 (MQ=255) aTCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAg < 1:340895/104‑1 (MQ=255) aTCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt < 2:545647/126‑1 (MQ=255) aTCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt > 1:545647/1‑126 (MQ=255) aTCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCa > 1:441994/1‑139 (MQ=255) ccAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCaa < 1:154570/139‑1 (MQ=255) ccAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCaa < 1:187567/139‑1 (MQ=255) cGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGaa > 1:111771/1‑100 (MQ=255) cGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGaa < 2:111771/100‑1 (MQ=255) aGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGaa > 1:374887/1‑95 (MQ=255) aGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGaa < 2:374887/95‑1 (MQ=255) aGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGa > 1:66039/1‑139 (MQ=255) gCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGAt > 2:293498/1‑139 (MQ=255) ccGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGAttt < 1:360220/139‑1 (MQ=255) ggTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGaa > 1:31410/1‑90 (MQ=255) ggTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGaa < 2:31410/90‑1 (MQ=255) gTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCAt < 1:305626/139‑1 (MQ=255) gTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCAt < 2:66039/139‑1 (MQ=255) gTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt > 2:526737/1‑101 (MQ=255) gTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGt < 1:526737/101‑1 (MQ=255) cTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGa < 1:74579/113‑1 (MQ=255) cTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGa > 2:74579/1‑113 (MQ=255) cTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaataa < 1:204785/139‑1 (MQ=255) cTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaat > 2:204785/1‑137 (MQ=255) cTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaataa < 2:441994/139‑1 (MQ=255) cGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaa < 1:414035/56‑1 (MQ=255) cGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaa > 2:414035/1‑56 (MQ=255) cGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGaact > 2:319530/1‑138 (MQ=255) gTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTAc > 1:122804/1‑139 (MQ=255) gTTTGTACTACCACAAAGAGATTGGCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTAc > 1:159198/1‑139 (MQ=255) tGTACTACCACACAGAGATTGGCCCCCGGATTAAGACAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTAccc > 2:207078/1‑137 (MQ=255) tactacCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCaaa > 1:46504/1‑138 (MQ=255) ccACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGaactaa > 1:411380/1‑112 (MQ=255) ccACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGaactaa < 2:411380/112‑1 (MQ=255) cacaAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATc < 2:191535/139‑1 (MQ=255) agaTTGGCCCCCGGATTAAAAAAAGTAATCACGaa > 2:290359/1‑35 (MQ=255) agaTTGGCCCCCGGATTAAAAAAAGTAATCACGaa < 1:290359/35‑1 (MQ=255) aTTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGAtt < 2:405059/79‑1 (MQ=255) aTTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGAtt > 1:405059/1‑79 (MQ=255) ggCCCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTACATCATATTTCTACCAAAAATAATCGGGTGCgagag < 1:135672/138‑1 (MQ=255) | GCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAG > NZ_CP009273/341808‑342078 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |