Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,334,438 A→G 24.2% S131P (TCG→CCG)  nudD ← GDP‑mannose mannosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,334,4380AG24.2% 84.4 / 24.0 74S131P (TCG→CCG) nudDGDP‑mannose mannosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base A (35/21);  new base G (11/7);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

AAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAC  >  minE/1334385‑1334498
                                                     |                                                            
aaaaTAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTa                                               <  1:502303/69‑1 (MQ=255)
     aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGt                                          <  1:1001904/69‑1 (MQ=255)
     aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGt                                          <  1:1988604/69‑1 (MQ=255)
     aGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGt                                          <  1:724032/69‑1 (MQ=255)
                        aaCATTATCACTGGCGCGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3288847/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGc                         >  1:983981/1‑67 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:552811/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:641162/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3373542/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3349906/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:1038860/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3195427/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3159793/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:3135166/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:2972497/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:2317814/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:2080432/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:1739378/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:1551795/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:1437673/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCa                        >  1:1085484/1‑68 (MQ=255)
                        aaCATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTAATGCTGCTCATCCGGCa                        >  1:1128183/1‑68 (MQ=255)
                          cATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGt                      <  1:2523121/68‑1 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACa                  >  1:2820847/1‑67 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:1319184/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:2972164/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:1270234/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:3018909/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:481572/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:3118537/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:1151478/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:2511104/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:2417774/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:112998/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:896679/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:3158138/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:945119/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:3258425/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGAAACAGc                >  1:3378093/1‑69 (MQ=255)
                               tCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                >  1:225798/1‑68 (MQ=255)
                                cACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:498710/68‑1 (MQ=255)
                                cACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2896472/68‑1 (MQ=255)
                                cACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2505350/68‑1 (MQ=255)
                                cACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2102558/68‑1 (MQ=255)
                                cACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:95148/68‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:3317063/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:1167951/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:1653177/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:920669/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2869623/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2333093/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2659914/65‑1 (MQ=255)
                                   tGGCGAGCAGCGCGTCCGACGTCAGCCAACGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:65861/65‑1 (MQ=255)
                                    ggCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:2225593/64‑1 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTcatctt         >  1:2964729/1‑68 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:531452/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:1073927/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:2131622/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:2224605/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:802088/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:1582876/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:2789668/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:411628/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:400577/1‑69 (MQ=255)
                                       gAGCAACGCGTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCatcttc        >  1:292608/1‑69 (MQ=255)
                                        agcagcGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTcc       <  1:289626/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:1359736/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:2670075/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:806554/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:873079/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:3126667/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:1736024/69‑1 (MQ=255)
                                             cgcgTCCGGCGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCATCTTCCGCTAc  <  1:1366034/69‑1 (MQ=255)
                                               cgTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGc                <  1:1650504/53‑1 (MQ=255)
                                                     |                                                            
AAAATAGGCGCGGCTGTTAGCATGAACATTATCACTGGCGAGCAGCGCGTCCGACGTCAGCCAGCGGTAATCGTCATGCTGCTCATCCGGCAGTAACAGCTCTTCTTCCGATAC  >  minE/1334385‑1334498

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