Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,641,322 C→T 50.7% S319S (AGC→AGT ggt → gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,641,3220CT50.7% ‑4.5 / 62.6 69S319S (AGC→AGTggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (13/21);  new base T (13/22);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CCATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAATACGCCTACGCCGACCGCTCGGAATAT  >  minE/2641259‑2641385
                                                               |                                                               
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:1396951/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:399493/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:1904553/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:1863079/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:455998/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:2129363/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:3261019/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:127379/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:2342589/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:2520212/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:1166506/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcga                                                            >  1:1153412/2‑69 (MQ=255)
tcATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGCAGAAGTACGGCTTTGGCAGTGCcg                                                             >  1:538274/2‑68 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:3169212/69‑1 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:2568139/69‑1 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:213073/69‑1 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:3288293/69‑1 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:485042/69‑1 (MQ=255)
            aaTCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:1301854/69‑1 (MQ=255)
            aaTCCTCAATATTCCGGAAAACCTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAAt                                                <  1:3145921/69‑1 (MQ=255)
                     tatTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGa                                       >  1:2816812/1‑69 (MQ=255)
                     tatTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGa                                       >  1:2230892/1‑69 (MQ=255)
                     tatTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGa                                       >  1:2234848/1‑69 (MQ=255)
                              aaaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGgaga                               <  1:3358199/68‑1 (MQ=255)
                              aaaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGgaga                               <  1:1297889/68‑1 (MQ=255)
                               aaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:2260430/69‑1 (MQ=255)
                               aaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:2829666/69‑1 (MQ=255)
                               aaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:2865501/69‑1 (MQ=255)
                               aaCTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:3284466/69‑1 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCTGAGaaa                             <  1:440372/67‑1 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:343878/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:694863/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:3108234/67‑1 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:1150910/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:577503/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:967298/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:1223491/67‑1 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:1424997/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:1538958/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             >  1:2808215/1‑67 (MQ=255)
                                 cTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:228217/67‑1 (MQ=255)
                                   tCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGgaga                               <  1:2346942/63‑1 (MQ=255)
                                      atatGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAAtacgcc                       <  1:2244549/68‑1 (MQ=255)
                                      atatGAAGAAATAAGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAAtacgcc                       <  1:1937652/68‑1 (MQ=255)
                                           aagaaATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAAtacgcctacgcc                 >  1:3297849/1‑69 (MQ=255)
                                           aagaaATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAAtacgcctacgcc                 >  1:3414064/1‑69 (MQ=255)
                                                    ggCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGaaa                             <  1:2814485/48‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCTAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2880599/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATGCGCCTATGCCGACCGCTCGGAatat  <  1:2468854/68‑1 (MQ=37)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGTAatat  <  1:69728/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:1581933/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:3022625/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:871514/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:737509/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:73118/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:1234178/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:675432/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2560568/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2545132/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2972004/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:1748608/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2304746/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:1957495/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2014555/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:3259687/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2024625/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2888427/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2093237/68‑1 (MQ=255)
                                                           gCAGTGCCGACGCGATGCAAATCATGGCAGAAGCTGAGAAATACGCCTATGCCGACCGCTCGGAatat  <  1:2181940/68‑1 (MQ=255)
                                                               |                                                               
CCATATCGTACAAATCCTCAATATTCTGGAAAACTTCGATATGAAGAAATACGGCTTTGGCAGCGCCGATGCGATGCAAATCATGGCAGAAGCGGAGAAATACGCCTACGCCGACCGCTCGGAATAT  >  minE/2641259‑2641385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: