Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 977,525 G→A 18.3% S540S (AGC→AGT dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE977,5250GA18.3% 242.9 / 34.3 151S540S (AGC→AGTdcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base G (84/39);  new base A (21/7);  total (105/46)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

GCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAT  >  minE/977461‑977583
                                                                |                                                          
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:1984596/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:2042102/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:207932/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:3170036/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:2137584/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:526430/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:1130025/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:1589048/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:1554782/1‑69 (MQ=255)
gCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg                                                        >  1:2654461/1‑69 (MQ=255)
 ccAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGc                                                       >  1:33531/1‑69 (MQ=255)
 ccAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGc                                                       >  1:1710027/1‑69 (MQ=255)
 ccAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGc                                                       >  1:453015/1‑69 (MQ=255)
 ccAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGc                                                       >  1:1195012/1‑69 (MQ=255)
 ccAGCGAATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGc                                                       >  1:443020/1‑69 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2855824/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:134540/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2798905/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:794457/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:1530517/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2459495/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2983953/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2196852/69‑1 (MQ=255)
   aGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt                                                     <  1:2922640/69‑1 (MQ=255)
         tatCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               <  1:59086/69‑1 (MQ=255)
         tatCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               <  1:1019444/69‑1 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          <  1:2295023/56‑1 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:3092846/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          <  1:3305263/56‑1 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:2515075/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:335872/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:1026717/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:2516472/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:1806399/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          >  1:176029/1‑56 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCt                                                          <  1:1598574/56‑1 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAtt                                             >  1:28339/1‑69 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAtt                                             >  1:2729944/1‑69 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAtt                                             >  1:1359952/1‑69 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAt                                              >  1:3047058/1‑68 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAt                                              >  1:410544/1‑68 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAt                                              >  1:1090692/1‑68 (MQ=255)
           tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAt                                              >  1:1642023/1‑68 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGc                                                >  1:1173871/1‑57 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               >  1:2319653/1‑58 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               >  1:2015705/1‑58 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               >  1:2049957/1‑58 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCa                                               >  1:2102106/1‑58 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGc                                    >  1:1078988/1‑69 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGc                                    >  1:348859/1‑69 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGc                                    >  1:2914876/1‑69 (MQ=255)
                    gcggcgCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCCGGCATTACGCATTTTCTGTTGc                                    >  1:1577119/1‑69 (MQ=255)
                                  ttCGCTCATCTCTTACCCTTTGTTGAACAGGCTGGCATTACTCa                                               <  1:476386/44‑1 (MQ=38)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:3165033/69‑1 (MQ=255)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:1267039/69‑1 (MQ=255)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:632177/69‑1 (MQ=255)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:3341267/69‑1 (MQ=255)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:376365/69‑1 (MQ=255)
                                  ttCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:1152535/69‑1 (MQ=255)
                                    cGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAgg                         >  1:1786753/1‑64 (MQ=255)
                                    cGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAgg                         >  1:1379855/1‑64 (MQ=255)
                                    cGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAgg                         >  1:1507889/1‑64 (MQ=255)
                                    cGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAgg                         >  1:3237462/1‑64 (MQ=255)
                                    cGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCAt                      <  1:2686172/67‑1 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTGAGGCATtgc                   >  1:1060269/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:3201526/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:3223760/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:1101073/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:596261/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:1255391/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:2750454/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:1336297/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:1541722/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:2342984/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:2340616/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:524421/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtgc                   >  1:733131/1‑69 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtg                    >  1:793010/1‑68 (MQ=255)
                                     gCTCATCTCATACCCTTTATTGAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATtg                    >  1:1760342/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCATTTCGTCAGGCATTGCTGCCCCGctct        >  1:2476270/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:753034/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:2407031/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:1070059/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:3124243/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:1740445/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:2577732/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctct        >  1:2321485/1‑68 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctc         >  1:824416/1‑67 (MQ=255)
                                                 cccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGctat        >  1:1101133/1‑66 (MQ=255)
                                                  ccTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTGGTCAGGCATTGCTGCCCCGCTCTg       >  1:1961856/1‑68 (MQ=255)
                                                   cTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGAt     <  1:2152458/69‑1 (MQ=255)
                                                      tgttgAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAt  >  1:2585276/1‑69 (MQ=255)
                                                      tgttgAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAt  >  1:65144/1‑69 (MQ=255)
                                                      tgttgAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAt  >  1:1544772/1‑69 (MQ=255)
                                                      tgttgAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAt  >  1:2861823/1‑69 (MQ=255)
                                                      tgttgAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAt  >  1:3201456/1‑69 (MQ=255)
                                                      tattgAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATTGCGGCCCCGCTCTg       <  1:3313959/62‑1 (MQ=255)
                                                      tattgAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATTGCGGCCCCGCTCTg       <  1:1935502/62‑1 (MQ=255)
                                                      tattgAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATTGCGGCCCCGCTCTg       <  1:1608971/62‑1 (MQ=255)
                                                      tattgAACAGACTGGCATTACGCATTTTCTGTTGAAGTTCGTCAGGCATTGCGGCCCCGCTCTg       <  1:3297709/62‑1 (MQ=255)
                                                                |                                                          
GCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTTTCTGTTGCAGTTCGTCAGGCATTGCTGCCCCGCTCTGATAAT  >  minE/977461‑977583

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: