Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,075 A→T 58.5% V524V (GTA→GTT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0750AT58.5% 36.0 / 74.9 77V524V (GTA→GTTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (13/19);  new base T (20/25);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTA  >  minE/190014‑190118
                                                             |                                           
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:514741/1‑55 (MQ=255)
                cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                    >  1:27723/1‑55 (MQ=255)
                cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                    >  1:154331/1‑55 (MQ=255)
                cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                      >  1:1257167/1‑69 (MQ=255)
                   ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGTGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:894249/65‑1 (MQ=255)
                   ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:106897/65‑1 (MQ=255)
                   ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:2644101/65‑1 (MQ=255)
                   ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                       <  1:463346/65‑1 (MQ=255)
                       gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctt               <  1:1834328/69‑2 (MQ=255)
                       gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:621761/69‑1 (MQ=255)
                       gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:2675369/69‑1 (MQ=255)
                       gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:1821718/69‑1 (MQ=255)
                       gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:149214/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:2442649/1‑67 (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)
                               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:1419785/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2926966/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:616781/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3077589/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:592151/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1468180/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1233416/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3224580/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2680412/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3314474/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:953241/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:3352934/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2854087/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2842915/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1607158/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2668939/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2663649/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1846561/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1888849/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2403666/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:1904420/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:870975/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2151660/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca  <  1:2072512/69‑3 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:114925/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:566402/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:838202/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:397673/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3317817/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:3243803/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)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2485818/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2225287/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:211024/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:209702/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:2017844/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1791176/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1528336/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1307497/1‑67 (MQ=255)
                                    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg    >  1:1278889/1‑67 (MQ=255)
                                                         tCGTAGTAGAGAAAACCGGCCCTTATAACctgctg               <  1:972656/35‑1 (MQ=255)
                                                             |                                           
TGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTA  >  minE/190014‑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: