| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,984,441 | 0 | A | C | 47.3% | 79.6 / 20.7 | 55 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (25/4); new base C (0/26); total (25/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984292‑2984591 | gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCacac < 1:855699/151‑1 (MQ=255) atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACAccc < 1:616138/151‑1 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTg < 1:814709/151‑1 (MQ=255) gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTgg < 1:63229/150‑1 (MQ=255) gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGc > 1:739977/1‑151 (MQ=255) gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCccg > 2:417599/1‑151 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCccgccg > 2:56113/1‑151 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCccgccg > 2:251891/1‑151 (MQ=255) ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGAt > 1:259326/1‑150 (MQ=255) aaaaaGGGGGGAGCCAAAAAAATGTATGTCGTCTTGCCAGCCATATTTATTCAACAAATTAATGACGGCAGATATTTATGGGA‑AAGGAGGAAATGCAGAGCTATTCGTCACTTTAAAGGGTCAAT‑GCCACACCCTGCCCCCCGCCGATa < 2:1010804/149‑1 (MQ=255) aaggAGGGACCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATCCAAAAAATAAAGGCGGGAAGATATTTATGGGA‑AAGGAGAAAAGGCAGACCTTTTCTTCCCTTTAAATGTTAAGT‑GCCCCCCCCTGCCCCCCGCCGATACAAAg < 2:880017/148‑1 (MQ=255) aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACaa > 2:995225/1‑149 (MQ=255) gggAGCCAAAAAAATGTATGCCGACTGGCCAGCAATAATTATTCAACAGATTAAGGACGGCAGATTTTTACCGGT‑AAGGAGAAAATGAAGACCTTTTCGTCCCGTTAATTGGACAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGtt < 2:1045179/150‑1 (MQ=255) cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc > 1:166517/1‑151 (MQ=255) aaaaaaaTGTATGTCGATTTGCCAGCAATCTTTATTCAACAGATTAATGACGGAAAATATTTACCGGA‑AAGGAGGAAATGCAGAGCTTTTCGTCACTTTCAAGGGTCAAT‑CCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTcc < 1:13487/148‑1 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAAT‑GCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGt < 2:1129812/148‑1 (MQ=255) aaaaTGTTTGTCGACTGGCCAGCAATAATTATTCAAAAGATAAAGGCCGGCAAATATTAATGGGA‑AAGGAGAAAAGGAAGAGCTTTTCGTACCGTTCAATGGTAAAG‑CCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt < 2:854055/151‑1 (MQ=255) aaaaTGTATTTCGCCTTGCCACCCATATTTATTCAACATATTATTGACGGCAGATTTTTATCGGA‑AGGGAGAAAATGCAGGCCTTTTCATCACTTTCAAGGGACAAC‑CCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTcc < 2:1096212/148‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 2:1066514/1‑151 (MQ=255) aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg > 2:45373/1‑151 (MQ=255) tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc > 2:622994/1‑151 (MQ=255) tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAAATTTTTATCGGT‑AGGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAAT‑CCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg < 2:946646/150‑1 (MQ=255) aaCAAATTAATGCCGGCAGATTTTTATCGGT‑AAGGAGAAAATGCAGAGCTATTCGTCCCTTTCATTGGTCAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 2:911246/149‑1 (MQ=255) aCAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag > 2:785828/1‑149 (MQ=255) gATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCACTTCATGAACCAGAGaa > 2:863724/1‑150 (MQ=255) aTTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc > 1:676832/1‑151 (MQ=255) ttAATGCCGGCAAATATTTA‑CGGTAAAGGAGAAAATGCAGTGCTTTTCGTCCGTTTCATTGCCCATT‑CCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag < 2:254220/146‑1 (MQ=255) aaaggcggaaaatttttAGCGGA‑AAGGAGAAAAAGCAGACTTTTTCCTCCCGTTAAATGGTAAAT‑CCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc < 2:1027613/134‑1 (MQ=255) aaTGCCGGCGGTTTTTTATCGGT‑AAGGAGGAAATGCAGACCTTTTCGTCCCGTTAAATGGTAAAC‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg < 2:656389/151‑1 (MQ=255) ggCAGATTTTTATGGGT‑AAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAAGGGTCAAT‑GCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa < 1:707156/147‑1 (MQ=255) ggCAGATATTTATGGGA‑AGGGAGAAAATGCAGAGCTATTCGTCCCTTTC‑ATGGGCAATCCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt < 2:216799/148‑1 (MQ=255) ggCAGATATTTATCGGT‑AAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAAC‑CCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc < 1:146507/150‑1 (MQ=255) gCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGGAACCAGGATTCgg > 1:475866/1‑151 (MQ=255) tttATCGGT‑AAGGAGAAAAGGAAGAGCTTTTCGTCCCGTTCAATGGTCAAC‑CCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 2:186059/149‑1 (MQ=255) tttAGCGGT‑AAGGAGAAAATGCAGGGCTATTCGTCACTTTCAATGGTCAAT‑GCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg < 1:806129/150‑1 (MQ=255) tttAGCGGT‑AAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAAGGATCAAG‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 1:651673/149‑1 (MQ=255) tttAGCGGA‑AAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAAT‑CCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 1:576340/149‑1 (MQ=255) tCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc > 1:420054/1‑151 (MQ=255) aggggaagggagaaAATGAGGAGCTTTTCGTCCCGTTCAAGGGTAAAT‑GCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg < 1:1092453/143‑1 (MQ=255) tgggaaaGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAAC‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa < 2:259275/144‑1 (MQ=255) ggT‑AAGGAGGAAAGGCAGAGCTATTCTTCACGTTCAAGGGTCAAT‑GCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc < 2:274256/149‑1 (MQ=255) aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 1:427385/1‑150 (MQ=255) aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg > 2:22881/1‑151 (MQ=255) aaaagcgggatttttCGTCTCTTTAAACGGTAAAT‑CCCCCCCCCGGCCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt < 1:374232/137‑1 (MQ=255) aTGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAGCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc > 1:25837/1‑151 (MQ=255) gCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:390947/1‑151 (MQ=255) tATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 1:267919/1‑150 (MQ=255) aTTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt > 1:80500/1‑151 (MQ=255) ttCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGTGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc > 2:1048140/1‑151 (MQ=255) ttCGCCACTTTCAATGGTCAAT‑CCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa < 1:427705/149‑1 (MQ=255) tCGTCACGTTCAATGGTCAGC‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc < 2:226714/150‑1 (MQ=255) taacgttgatggatCACT‑GCCCCCCCCTGCCCGCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg < 1:393482/137‑1 (MQ=255) cACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 2:968511/1‑151 (MQ=255) ttCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAACGATGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGAGACCGAGTGCGATCGCGccacc > 2:450089/1‑150 (MQ=255) tCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt > 1:930278/1‑151 (MQ=255) ccACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc < 2:1117695/150‑1 (MQ=255) acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc > 1:797741/1‑151 (MQ=255) | GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984292‑2984591 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |