Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 87,821 G→A 36.7% V283V (GTG→GTA ddlB → D‑Alanine:D‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE87,8210GA36.7% 67.4 / 66.3 106V283V (GTG→GTAddlBD‑Alanine:D‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (42/25);  new base A (2/37);  total (44/62)
Fisher's exact test for biased strand distribution p-value = 1.00e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 106 total aligned reads displayed.

CAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGAA  >  minE/87755‑87873
                                                                  |                                                    
cAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGcc                                                    <  1:2556163/69‑1 (MQ=255)
cAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGcc                                                    <  1:2453952/69‑1 (MQ=255)
cAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGc                                                     >  1:1553716/1‑68 (MQ=255)
cAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTAACCGGGTATGACCAGCCACAGCCTGGTGcc                                                    >  1:372876/1‑69 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2196436/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1644362/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1795218/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1810736/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1929358/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1929925/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1972048/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:201733/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2037704/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2142380/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:103136/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2248565/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2341452/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:2447082/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:3196969/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:3306706/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:388803/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:496534/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1483269/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1414063/69‑1 (MQ=255)
   cGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1457685/69‑1 (MQ=255)
    gATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGAt                                                 <  1:1635044/68‑1 (MQ=255)
        gACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATggcgg                                            >  1:2686526/1‑69 (MQ=255)
        gACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATggcgg                                            >  1:1114102/1‑69 (MQ=255)
               ttATCTTCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:636215/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:2856473/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:207868/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:394753/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:1345730/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:245538/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:2493408/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:1459094/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:325545/1‑57 (MQ=255)
               ttATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGAt                                                 >  1:164873/1‑57 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:3260578/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1067100/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1099735/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:3121079/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:534404/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:85474/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1605129/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1426082/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1466677/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:837593/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:831214/66‑1 (MQ=255)
                        ggAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1911491/66‑1 (MQ=255)
                         gAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:2166453/65‑1 (MQ=255)
                         gAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:318924/65‑1 (MQ=255)
                         gAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:1350501/65‑1 (MQ=255)
                                 tACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:1739899/1‑68 (MQ=255)
                                 tACCTCACCGGGTATGACCAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTtc                     >  1:1930627/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               >  1:1228189/1‑56 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               >  1:3270267/1‑56 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGt                              >  1:441641/1‑57 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:2101581/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:765245/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:1094879/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:265170/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:2213447/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:1396566/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTtct                    >  1:2920693/1‑67 (MQ=255)
                                  aCCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACCTCAGGCAGGTATGAGCTtct                    >  1:1079634/1‑67 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGa                          >  1:2794610/1‑60 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:326596/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:3336998/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:1520563/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:1368489/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:1543735/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:2690628/1‑69 (MQ=255)
                                   ccTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGc                 >  1:1128905/1‑69 (MQ=255)
                                      cACCGGGTATGACCAGCCACAGCCTGGGCCGATGGCGGCACGTcaggcagg                               <  1:572835/51‑1 (MQ=255)
                                          gggTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgt          <  1:2433792/69‑1 (MQ=255)
                                          gggTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgt          <  1:2379880/69‑1 (MQ=255)
                                          gggTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgt          <  1:528902/69‑1 (MQ=255)
                                                  ccAGCCACAGCCTGGTGCCGATGGCGGCACGTcaggcagg                               <  1:2818456/40‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:874022/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:580533/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:1415244/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:1377202/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:1893510/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:2989046/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:1011655/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:395865/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:2311815/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:234481/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:3280183/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:3269908/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:3202804/68‑1 (MQ=255)
                                                   cAGCCACAGTCTGGTACCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGaa  <  1:2708617/68‑1 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:2182439/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:2405973/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:2424660/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:3197151/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:1878406/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:184971/1‑49 (MQ=255)
                                                                gTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGgtag        >  1:3195655/1‑49 (MQ=255)
                                                                  |                                                    
CAGCGATGGACAGTTTTATCTGCTGGAAGCCAATACCTCACCGGGTATGACCAGCCACAGCCTGGTGCCGATGGCGGCACGTCAGGCAGGTATGAGCTTCTCGCAGTTGGTAGTACGAA  >  minE/87755‑87873

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: