Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 53,762 A→G 36.2% D342D (GAT→GAC imp ← exported protein required for envelope biosynthesis and integrity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE53,7620AG36.2% 97.5 / 94.0 149D342D (GAT→GACimpexported protein required for envelope biosynthesis and integrity
Reads supporting (aligned to +/- strand):  ref base A (43/52);  new base G (22/32);  total (65/84)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 149 total aligned reads displayed.

CCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACAC  >  minE/53695‑53818
                                                                   |                                                        
ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:2183345/69‑1 (MQ=255)
ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:3430738/69‑1 (MQ=255)
ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:2689377/69‑1 (MQ=255)
ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:3054471/69‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:986253/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:215365/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2077264/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2073292/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:2034162/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2288229/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:1962329/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:2313863/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:248832/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2631282/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:15113/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:754405/68‑1 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:591515/1‑68 (MQ=255)
       gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:137801/1‑68 (MQ=255)
         cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1496126/1‑69 (MQ=255)
         cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:2499399/1‑69 (MQ=255)
         cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1772788/1‑69 (MQ=255)
         cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1672675/1‑69 (MQ=255)
         cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:3125072/1‑69 (MQ=255)
                gAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:1742567/69‑1 (MQ=255)
                gAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:3035443/69‑1 (MQ=255)
                gAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:343134/69‑1 (MQ=255)
                          gttgcgtAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGAc                                <  1:3285998/68‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1017530/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:697937/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:65872/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3144948/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2270104/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2311756/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:335326/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3234362/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3179502/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3203727/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2828964/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2817122/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1125598/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:192655/69‑1 (MQ=255)
                            tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2166252/69‑1 (MQ=255)
                             gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1287654/68‑1 (MQ=255)
                             gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:631914/68‑1 (MQ=255)
                             gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1171011/68‑1 (MQ=255)
                             gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:102086/68‑1 (MQ=255)
                             gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3190794/68‑1 (MQ=255)
                              cgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:975568/67‑1 (MQ=255)
                                        aGTACTGGAACCGTACTTGTTATCGAAATCTTTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCgc                  >  1:620848/1‑67 (MQ=255)
                                        aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:1515753/1‑69 (MQ=255)
                                        aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:1868940/1‑69 (MQ=255)
                                        aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:696079/1‑69 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:38647/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:3126507/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:529168/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1245649/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:293543/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1193483/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:3378098/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1008402/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:2594279/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:776013/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:2523425/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:200293/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGGCGACGt               <  1:2380714/69‑1 (MQ=255)
                                          tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTACTCGACGt               <  1:2571563/69‑1 (MQ=255)
                                            cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:1458305/1‑68 (MQ=255)
                                            cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:193140/1‑68 (MQ=255)
                                            cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:1199632/1‑68 (MQ=255)
                                            cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:899636/1‑68 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTTGGATCGCTGACTTTg                            >  1:834786/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAgg                                         >  1:662616/1‑39 (MQ=38)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:473771/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:387248/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:841725/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1452030/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2022502/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:3195544/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2012637/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1919040/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1569004/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2820069/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1770557/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2765329/1‑52 (MQ=255)
                                              ggAACCGTATTTGTTATCGAAGTCATTCAAGTAGCTAGGATCGCTGACTTTg                            >  1:2935780/1‑52 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:132714/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:896301/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:2525945/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:2545619/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:2698557/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:563198/68‑1 (MQ=255)
                                              ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGAGCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:3141628/68‑1 (MQ=255)
                                                  ccGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGc       <  1:2608109/69‑1 (MQ=255)
                                                  ccGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGc       <  1:2248658/69‑1 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1836879/1‑69 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1606981/1‑69 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1429936/1‑69 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1378874/1‑69 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1316477/1‑69 (MQ=255)
                                                       cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:2069906/1‑69 (MQ=255)
                                                                   |                                                        
CCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACAC  >  minE/53695‑53818

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: