Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,534,073 G→C 33.3% V72V (GTG→GTC yfcN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,534,0730GC33.3% 84.0 / 64.4 84V72V (GTG→GTCyfcNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (21/7);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGAT  >  minE/1534007‑1534130
                                                                  |                                                         
acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         <  1:1514270/68‑1 (MQ=255)
acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         <  1:2584306/68‑1 (MQ=255)
acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         <  1:3106075/68‑1 (MQ=255)
 cAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGaaa                                                        <  1:1333011/69‑1 (MQ=255)
 cAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGaaa                                                        <  1:2199187/69‑1 (MQ=255)
  aGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         <  1:337820/67‑1 (MQ=255)
    gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                                     <  1:1551542/69‑1 (MQ=255)
    gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                                     <  1:1593718/69‑1 (MQ=255)
    gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                                     <  1:2186767/69‑1 (MQ=255)
    gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                                     <  1:2190577/69‑1 (MQ=255)
     cTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                                     <  1:1869473/68‑1 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2199969/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2218377/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2224488/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:3034801/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:3095200/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2965009/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2438598/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1729705/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2484031/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:2661651/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1072854/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:383993/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:128879/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:973003/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:86578/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1348984/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1365489/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1412495/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:583138/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1485954/1‑49 (MQ=255)
                    attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                                         >  1:1149314/1‑49 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                     >  1:1433596/1‑69 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                     >  1:3201637/1‑69 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                     >  1:2737347/1‑69 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                     >  1:162614/1‑69 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                     >  1:1188298/1‑69 (MQ=255)
                    attTCTCCGATGAGTTTCAGCCGGTATTAAATACCGAAGGTCCGGTg                                                           >  1:2533740/1‑47 (MQ=255)
                         tcCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCa                                 <  1:3212765/68‑1 (MQ=255)
                                 gTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCg                                            <  1:1695137/49‑1 (MQ=255)
                                   ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg                       <  1:1362953/68‑1 (MQ=255)
                                   ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg                       <  1:2707855/68‑1 (MQ=255)
                                   ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg                       <  1:3063179/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2565434/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2696780/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:3074387/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1346116/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:863651/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1429718/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1619204/68‑1 (MQ=255)
                                       gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1615601/68‑1 (MQ=255)
                                        ccGTTTTTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa                  <  1:1154672/68‑1 (MQ=255)
                                        ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa                  <  1:2504684/68‑1 (MQ=255)
                                        ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa                  <  1:1779427/68‑1 (MQ=255)
                                        ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa                  <  1:2607978/68‑1 (MQ=255)
                                        ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa                  <  1:1510733/68‑1 (MQ=255)
                                        ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1113906/67‑1 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:3137782/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:1200591/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:774439/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:709485/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:1438039/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:520871/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:1788819/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:1859499/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:2282599/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:259949/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:281739/1‑69 (MQ=255)
                                             attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg            >  1:2940162/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2431136/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1337647/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2444715/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:525383/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2833084/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:1954245/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga                   <  1:2316064/54‑1 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2913065/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:1538934/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2704072/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2048704/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2761621/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2147620/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:1381999/1‑69 (MQ=255)
                                                     ccGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCg    >  1:2282250/1‑69 (MQ=255)
                                                       gAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGat  <  1:760823/69‑1 (MQ=255)
                                                                  |                                                         
GCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGAT  >  minE/1534007‑1534130

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: