Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,088,157 | G→T | 44.7% | G48C (GGC→TGC) | hisH → | imidazole glycerol phosphate synthase, glutamine amidotransferase subunit with HisF |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,088,157 | 0 | G | T | 44.7% | 11.4 / 50.6 | 45 | G48C (GGC→TGC) | hisH | imidazole glycerol phosphate synthase, glutamine amidotransferase subunit with HisF |
Reads supporting (aligned to +/- strand): ref base G (14/11); new base T (9/11); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01 |
CTGTAATGAACGTGGTGATCCTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTCGAC > CP009273/2088011‑2088297 | cTGTAATGAACGTGGTGATCCTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGc > 2:378106/1‑149 (MQ=255) gTAATGAACGTGGTGATCCTTGATACCGGCTGCGCTAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCAGACGTCGTGTTGCTGGCCGATTAACTGTTTTTACCCGGCGTTTGCAc < 2:355953/149‑1 (MQ=255) tGAACGTGGTGATCCTTGACACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTgcg < 2:256736/149‑1 (MQ=255) aacgccgcgATCCCTGATACCGCCTCCGCCAACCTGAACTCGGTGAAGTCTCCCATTGCGCCTCACGGTGATTAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCAc < 1:249498/135‑1 (MQ=255) gtggtgATCCTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAgc > 2:346651/1‑149 (MQ=255) gtgATCCTTGATACCGGCTGCGCCAGCCTGGACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCACAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAgcggc > 2:328766/1‑149 (MQ=255) tgATCCTTGATACCGGCTGCGACAACCTGCGATCGGAGTGGTATCACCGACCCCCGCCCTGTTAGGAACCCAAAGCCAGTCGTGCCCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCaa < 1:402210/143‑1 (MQ=255) gATCCTTGTTACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGa < 2:190207/149‑1 (MQ=255) gATCCTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGa < 2:519132/149‑1 (MQ=255) cTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGa > 2:579142/1‑149 (MQ=255) cGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGAGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTAc < 2:594572/149‑1 (MQ=255) ccTGCGCCAACCTGAACTCGGTGAAGTCTGCGATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAg > 1:418442/2‑143 (MQ=255) gcCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAgcg > 2:503363/1‑149 (MQ=255) cTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCt > 2:93359/1‑149 (MQ=255) cTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGGACGTGGGCGCGAGCt > 2:222713/1‑149 (MQ=255) tGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATAAGGTACGTGAGCGCGAGCTg > 2:441809/1‑149 (MQ=255) gTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGAtct > 2:515962/1‑149 (MQ=255) tGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAATCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGCCACTGCGCAGCCGGCGCTGGCTCACGTACGTGAGCGCGAGCTGTTTGAgctc > 2:53774/1‑149 (MQ=255) aGTCTGCCATTGCGCGTCACGGTTATGAACCCTAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCtcatc < 2:434651/149‑1 (MQ=255) aGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGAt < 1:394297/143‑1 (MQ=255) cTGCCATTGCGCGTCACGTTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCaaa > 2:94207/1‑149 (MQ=255) tGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGGACGTGAGCGCGAGCTGTTTGATCTCATCAAAg > 2:592027/1‑149 (MQ=255) tGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAg > 2:581539/1‑149 (MQ=255) cATTGCGGGTCCCGGTTATGAACCCAAGGCCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCt < 2:167361/149‑1 (MQ=255) tCACGGTTATGAACCCAAAGAGAACTATGAACGAGCCGGTGCGTTGTTGGCGGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAAc < 2:400499/149‑1 (MQ=255) ttATGAACCCAAAGGCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGACCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGc < 2:405410/149‑1 (MQ=255) ttAGGAACCCAAGGCCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAAc < 1:179805/143‑1 (MQ=255) tGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTgg > 2:170062/1‑149 (MQ=255) aggcAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCt < 2:234758/146‑1 (MQ=255) aGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCGGGGCATCTGCt > 2:227711/1‑149 (MQ=255) gTCAGCCGTGGCCCGGCCGTCGTGTTGCCGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGACCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCtt < 2:318487/149‑1 (MQ=255) gACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCAGGGCATCTGCTTAGGGATGCAActgct > 2:444371/1‑149 (MQ=255) aCGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCaa < 1:534774/143‑1 (MQ=255) gtTGCTGGCCGTTACACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGggcggc < 2:51282/149‑1 (MQ=255) gtTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATAAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAAACGGTGCTGGCCATCTGCTTAGGGATGCAACTGCAgg > 1:76073/1‑143 (MQ=255) cTGGCCGATAAACTGTTTTTACCCGGCGTTGTCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAAACGGTGCTGGGCATCTGCTTAGGGATTCCACTGCTGGggcg > 1:27573/1‑143 (MQ=255) gCCGATAAACTGTTTTTACCCGGCGGTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGa < 2:84119/149‑1 (MQ=255) aTAAACTGGTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAagag < 2:613750/149‑1 (MQ=255) ttttACCCGGCGTTTGCACTGCGTAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGc > 2:351482/1‑149 (MQ=255) ttttACCCGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTCGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGc > 1:136782/1‑143 (MQ=255) cccGGCGTTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACTGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTCGa < 2:88978/149‑1 (MQ=255) accGGCGGTTGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGCACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTCGa < 2:88710/148‑1 (MQ=255) cGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGGAACTGCTGGGGCGGCGCAGCGAAGAGGGCAACGGCg > 1:113174/1‑143 (MQ=255) gCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTc < 1:84233/143‑1 (MQ=255) ttGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTCGAc > 1:261515/1‑143 (MQ=255) | CTGTAATGAACGTGGTGATCCTTGATACCGGCTGCGCCAACCTGAACTCGGTGAAGTCTGCCATTGCGCGTCACGGTTATGAACCCAAAGTCAGCCGTGACCCGGACGTCGTGTTGCTGGCCGATAAACTGTTTTTACCCGGCGTTGGCACTGCGCAAGCGGCGATGGATCAGGTACGTGAGCGCGAGCTGTTTGATCTCATCAAAGCCTGTACCCAACCGGTGCTGGGCATCTGCTTAGGGATGCAACTGCTGGGGCGGCGCAGCGAAGAGAGCAACGGCGTCGAC > CP009273/2088011‑2088297 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |