Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,053 T→A 53.0% F517Y (TTC→TAC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0530TA53.0% 9.6 / 98.0 85F517Y (TTC→TAC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base A (20/25);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTA  >  minE/189994‑190118
                                                           |                                                                 
gATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcgacga                                                          >  1:1037892/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:2039235/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:2081719/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:683068/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:2357575/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:3121509/1‑69 (MQ=255)
  tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAAc                                                        >  1:1459216/1‑69 (MQ=255)
    gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGt                                                      >  1:2998876/1‑69 (MQ=255)
    gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGt                                                      >  1:1875943/1‑69 (MQ=255)
                    tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                       >  1:319451/1‑68 (MQ=255)
                            taggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTATTAGAGAAAACCGGc                              <  1:950277/68‑1 (MQ=255)
                              ggggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc                             >  1:1894187/1‑68 (MQ=255)
                              ggggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGc                              <  1:197886/67‑1 (MQ=255)
                                  aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtat                        <  1:51934/69‑1 (MQ=255)
                                   tCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:618687/69‑1 (MQ=255)
                                   tCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:196227/69‑1 (MQ=255)
                                    cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                    >  1:27723/1‑55 (MQ=255)
                                    cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                    >  1:154331/1‑55 (MQ=255)
                                    cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                    >  1:514741/1‑55 (MQ=255)
                                    cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                      >  1:1257167/1‑69 (MQ=255)
                                       ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGTGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:894249/65‑1 (MQ=255)
                                       ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:2644101/65‑1 (MQ=255)
                                       ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:463346/65‑1 (MQ=255)
                                       ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:106897/65‑1 (MQ=255)
                                           gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctt               <  1:1834328/69‑2 (MQ=255)
                                           gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:2675369/69‑1 (MQ=255)
                                           gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:1821718/69‑1 (MQ=255)
                                           gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:149214/69‑1 (MQ=255)
                                           gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:621761/69‑1 (MQ=255)
                                            cgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:228627/68‑1 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                        >  1:2641366/1‑36 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:2144436/61‑1 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc         >  1:431183/1‑67 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc         >  1:1491608/1‑67 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc         >  1:872027/1‑67 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc         >  1:1747043/1‑67 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc         >  1:2442649/1‑67 (MQ=255)
                                                   aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtct          >  1:3346748/1‑66 (MQ=255)
                                                    gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCttt      <  1:2359602/69‑1 (MQ=255)
                                                       aaaTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt   <  1:786128/69‑1 (MQ=255)
                                                        aaTACCTCGACGCACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3303229/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATTGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:206573/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCgg                               >  1:221566/1‑40 (MQ=16)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3077589/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2926966/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:953241/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:592151/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:870975/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3224580/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3314474/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3352934/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:616781/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2403666/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1888849/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1904420/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1846561/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2072512/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1607158/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2151660/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1468180/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1419785/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2663649/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1233416/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2668939/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2680412/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2854087/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2842915/69‑3 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1307497/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:114925/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:838202/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1278889/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1528336/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1791176/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:566402/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2017844/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:209702/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:397673/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:211024/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3317817/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2225287/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3243803/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2485818/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3196280/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3183370/1‑67 (MQ=255)
                                                        aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3028215/1‑67 (MQ=255)
                                                           |                                                                 
GATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTA  >  minE/189994‑190118

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: