Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,237,409 A→G 40.0% K318K (AAA→AAG purD → phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,237,4090AG40.0% 43.4 / 74.5 105K318K (AAA→AAGpurDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (39/24);  new base G (13/29);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 2.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 108 total aligned reads displayed.

AACCCAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCGGGTGGA  >  minE/2237340‑2237475
                                                                     |                                                                  
aaCCCAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTTGCCTGGCGGCCTGTGAAAGCaaa                                                                    >  1:3219507/1‑69 (MQ=255)
   ccAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAAc                                                                   >  1:11921/1‑68 (MQ=255)
   ccAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAAc                                                                   >  1:265700/1‑68 (MQ=255)
   ccAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAAc                                                                   >  1:139084/1‑68 (MQ=255)
    cAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTg                                                                 >  1:955298/1‑69 (MQ=255)
    cAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTg                                                                 >  1:559786/1‑69 (MQ=255)
    cAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTg                                                                 >  1:2934985/1‑69 (MQ=255)
    cAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTg                                                                 >  1:725168/1‑69 (MQ=255)
    cAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTg                                                                 >  1:2448942/1‑69 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:535009/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:2220613/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:595949/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:1611180/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:1441341/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:627140/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:1649603/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:1163355/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:494553/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTgg                                                                >  1:860796/1‑66 (MQ=255)
        cGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAGGCAAGCTg                                                                 >  1:1261923/1‑65 (MQ=255)
         gATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACg                                                             >  1:1847876/1‑68 (MQ=255)
             aTGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGaaa                                                        >  1:2725098/1‑69 (MQ=255)
             aTGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGaaa                                                        >  1:2132314/1‑69 (MQ=255)
                     cATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGa                                                >  1:3390681/1‑69 (MQ=255)
                              cGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTgg                                            >  1:532309/1‑64 (MQ=255)
                              cGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTgg                                            >  1:513757/1‑64 (MQ=255)
                              cGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTgg                                            >  1:3174990/1‑64 (MQ=255)
                              cGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTgg                                            >  1:3273545/1‑64 (MQ=255)
                              cGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTgg                                            >  1:2529174/1‑64 (MQ=255)
                                           ctctGCCTGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTtctctc                          <  1:2299520/69‑1 (MQ=255)
                                           ctctGCCTGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTtctctc                          <  1:3316613/69‑1 (MQ=255)
                                           ctctGCCTGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTtctctc                          <  1:1443787/69‑1 (MQ=255)
                                           ctctGCCTGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTtctctc                          <  1:620096/69‑1 (MQ=255)
                                           ctctGCCTGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGAAGAACGCGCTtctctc                          <  1:1958524/69‑1 (MQ=255)
                                                  ttgcggcCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:524317/66‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTTGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1159029/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGATTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1823209/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:3293892/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2740741/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2980341/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2711746/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2989826/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:718787/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:705742/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:307309/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2323808/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:563340/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2279002/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1500688/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:325462/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1628108/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2125037/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:2091371/68‑1 (MQ=255)
                                                  tGGCGGCCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1986884/68‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:956692/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:849635/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:3289885/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:2862189/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1140151/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1234886/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1510867/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:2733066/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1547935/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1785760/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:1933879/69‑1 (MQ=255)
                                                   ggcggcCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtggt                  <  1:2228764/69‑1 (MQ=255)
                                                     cggcCTGTGAAGGCAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:455179/65‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGc                       >  1:2836816/1‑59 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:1070313/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:1494484/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:2895964/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:2718850/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:3165671/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:1595981/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:1486269/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:1472388/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:1380999/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:1264610/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:439337/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:1709605/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:762670/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:1915717/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:3028899/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:2179584/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             <  1:2491556/69‑1 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATgg             >  1:2754343/1‑69 (MQ=255)
                                                        ccTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGCTGATgg             <  1:443223/69‑1 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAAGGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:2368558/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:1601417/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:1651265/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:485021/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:2699861/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:2749609/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:3008525/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:2251838/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCgg       >  1:889520/1‑69 (MQ=255)
                                                              aaaGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCg        <  1:1332570/68‑1 (MQ=255)
                                                                  cAAGCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCgtg                    <  1:1571970/52‑1 (MQ=255)
                                                                   aaaCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCGGGTGGa  <  1:1138875/69‑1 (MQ=255)
                                                                   aaaCTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCGGGTGGa  <  1:1033484/69‑1 (MQ=255)
                                                                     |                                                                  
AACCCAGCCGATTATGCTGCGCATGAAGTCCGATCTGGTTGAGCTCTGCCTGGCGGCCTGTGAAAGCAAACTGGACGAGAAAACGTCCGAGTGGGATGAACGCGCTTCTCTCGGCGTGGTGATGGCTGCGGGTGGA  >  minE/2237340‑2237475

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: