| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,443,277 | 0 | A | C | 50.0% | 55.3 / 18.4 | 46 | V274G (GTA→GGA) | secY | preprotein translocase membrane subunit |
| Reads supporting (aligned to +/- strand): ref base A (22/1); new base C (0/23); total (22/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC > NC_000913/3443128‑3443425 | cAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑Ac > 1:146148/1‑151 (MQ=255) aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑Ac < 1:489546/150‑1 (MQ=255) aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACAGCCACCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGc > 1:130306/1‑150 (MQ=255) agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCat > 2:544801/1‑151 (MQ=255) taaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTc > 1:437714/2‑149 (MQ=255) ccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg > 2:668673/1‑151 (MQ=255) cagccaaaaaaTACAGCAAAATTGTTCCAAGCCAGTCCCAACCAGCCCCCCCCCCGAACCATGCC‑GCGATGGTCCCCGGAAACAAAAAAATACTGGAAGGAAAATTTGCGGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtgtgt < 1:22302/143‑1 (MQ=255) ggcggaaaagaaaGCAAATTTGTTGCCACCCAGTCCCAACCAG‑CCCCCCCCCCACCCATGGCGGCGGTGGCCGCGGGAAAAAAAAAAATACGGAAAGGGAATATGGCAGAGT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg < 2:390379/138‑1 (MQ=255) ggCTGAAAATACAGCAAATTTGTTGTCAGCCATTCCCAACCGGTCCCCCCCCCGAACCATGCC‑GCGATGGTCGCCGGGAAAAGAATAATCCTGGAGGGAAAGTTGGCCGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc < 1:668749/151‑1 (MQ=255) aaaaACAGCGAAATGGTTGTCAGCCAGCTCCAACCAGCCCCCCCCCGGACCCATGCC‑GGGATGCTCGCGGGAAACAAAATAATACTGAAGGCGAAGTTCGCCGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 2:172187/146‑1 (MQ=255) aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg > 1:370478/1‑151 (MQ=255) tACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc > 1:62842/1‑149 (MQ=255) aGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt > 1:290217/1‑151 (MQ=255) gaattggttgaCAGCCAGTCCCACCCGCACCCCCCCCCCACCCATGCC‑GGGGTGGTCGCCGAGAAGAAAATAATCCTAAAAGCGAAAAGGGGCGGTA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATa < 1:355425/143‑1 (MQ=255) aaaTTTTGGCCAGCCAATCCCAACAAGCCCCCCCCCCAAACAAGGCC‑GCAGTGGCCGCCGAAAAAAAAATAATTTTGGAAGGAAAGTTGGCCGGGT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 1:423684/147‑1 (MQ=255) aaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac > 2:242938/1‑150 (MQ=255) ttgttgTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacg > 1:1141/1‑151 (MQ=255) ttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg > 2:280697/1‑151 (MQ=255) cAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac > 1:577908/1‑149 (MQ=255) cAGCCAGTCCCAACCAGTCCCGCCCCAGAACCATGAC‑GCGGTGGTCGCCGGGAACAGAAAATTATTGAAATAGAAGATTGCGGTTT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc < 1:100772/150‑1 (MQ=255) gCCATTCCCACCCAGCACCCCCCCCAACCAATGCC‑GGGAGGGTCGCCGGAAACAAAAAAATACGGGAAGC‑AAATTTTCCGGGATTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt < 1:657690/150‑1 (MQ=255) acccagccccGCCCCCGACCCAGGGC‑GGGAGGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGTTTGCCGGTT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc < 1:658494/139‑1 (MQ=255) ccAGTCCCCCCCCCAAACCATGCG‑GGGATGGTCGCGGGAAACAAAAAAATACTGAAAGCGAAGATTGCCGGAT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg < 2:321283/151‑1 (MQ=255) aatcccccccccAAACCAGGAG‑GCGAGGGTCGCGGGAAAAAAAAAAATACGGAAAGGGAAGTTGGCGGGTT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg < 1:169762/144‑1 (MQ=255) aGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg > 1:5073/1‑151 (MQ=255) ccccccccAAACCATGAC‑GCGGTGGTCGCCGGAAACAAAAAAATCCGGGAAGCAAATTTGGCCGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa < 1:421774/135‑1 (MQ=255) cccccaaaaCAATGCC‑GCGATGGGGGCGGGGAAAAAAATAATACTGGAAGGAAATTTTGCGGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt < 1:576834/139‑1 (MQ=255) cccccGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt > 2:501262/1‑150 (MQ=255) ccGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc > 1:51876/1‑150 (MQ=255) aaCCATGCC‑GCGGTGGTCGCGGGAAAAAGAATAATACTGGAAGCGAAGTTTGCGGGTT‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc < 1:215505/147‑1 (MQ=255) ccATGCC‑GCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATGGCCGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt < 2:58732/150‑1 (MQ=255) tGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc > 1:155376/1‑151 (MQ=255) ggTCGCCGGGAAAAAAAAAATATTGGAGGCGAAGTTTGCGGGGA‑TC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa < 2:44475/148‑1 (MQ=255) gTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc > 2:553271/1‑151 (MQ=255) cGCCGGGAACAGAATAATACTGGAAGCGAAGATTTCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGCCGTTGACCACGc > 1:585063/1‑151 (MQ=255) ggaaaaaaaaaaataCTGAAAGCAAAGAGTGCGGGGA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca < 1:677638/139‑1 (MQ=255) aaaaaaataatcggggAAGCGAAG‑TTGGCGGGT‑TTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa < 1:465155/137‑1 (MQ=255) aaaaaaataatACGGGAAGCAAAGATTGCGGGTA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca < 2:270482/145‑1 (MQ=255) aCAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa > 1:672704/1‑150 (MQ=255) aCTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa > 2:359633/1‑149 (MQ=255) ggAAGCGAAGTTTGCGGGGT‑TC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa < 1:251705/149‑1 (MQ=255) aaccaaatatggCCGGAA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg < 1:506995/135‑1 (MQ=255) aaGCGAAGATGGCCGGTA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg < 2:644655/150‑1 (MQ=255) gggA‑TT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat < 2:146152/147‑1 (MQ=255) ggA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatacta > 2:500908/1‑150 (MQ=255) tt‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatac > 2:317865/1‑151 (MQ=255) | CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGAC‑GCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGA‑TT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATAC > NC_000913/3443128‑3443425 |
| 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: ‑ |