Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,931 A→C 92.2% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE214,9310AC92.2% 2763.4 / 123.9 1140E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (59/31);  new base C (590/459);  total (650/490)
Fisher's exact test for biased strand distribution p-value = 9.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 1145 total aligned reads displayed.

GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  minE/214866‑214995
                                                                 |                                                                
ctcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                               <  1:1227376/67‑1 (MQ=255)
 tcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                              >  1:2548820/2‑69 (MQ=255)
  cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:1200547/1‑69 (MQ=255)
      tCTGACATAAATGACCTCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:664383/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCTCTACCGCAACCACCgg                                                         <  1:567920/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:447040/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:1924896/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:2262867/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:1148743/69‑1 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        >  1:1449405/1‑69 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                        >  1:2661980/1‑69 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                        >  1:1888701/1‑69 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     <  1:2102034/67‑1 (MQ=255)
             tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:2299672/1‑67 (MQ=255)
             tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:1579999/1‑67 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                               <  1:2233233/54‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        <  1:1519471/61‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        <  1:2165017/61‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     >  1:1817147/1‑64 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:231975/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:3369963/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:981542/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:2814263/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:877416/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:791176/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:1000043/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:3392178/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:403538/69‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:2082538/69‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                    >  1:2974344/1‑65 (MQ=255)
                aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               >  1:746706/1‑69 (MQ=255)
                 tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:2467816/68‑1 (MQ=255)
                 tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:243189/68‑1 (MQ=255)
                   aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:493281/69‑1 (MQ=255)
                     cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:2762377/1‑69 (MQ=255)
                     cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:1672526/1‑69 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               >  1:2497086/1‑60 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               >  1:2788067/1‑60 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:2052850/69‑1 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:2640709/69‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     >  1:2863481/1‑51 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:3336517/56‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               >  1:3129739/1‑57 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:1643061/1‑69 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:3306773/69‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:470141/1‑69 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:970957/1‑69 (MQ=255)
                              cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                  >  1:240553/1‑68 (MQ=255)
                                cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:684311/69‑1 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         >  1:368087/1‑41 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         >  1:945325/1‑41 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:3035297/1‑69 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGCTGCGCACTATGAAg                            >  1:3159600/1‑69 (MQ=255)
                                     aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                            <  1:2937982/67‑1 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:388428/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:1751210/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:3058537/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:1550749/1‑68 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:1262439/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:1297389/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:1720655/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         >  1:640153/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcgc                                     <  1:106117/55‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:2372873/58‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:3194054/69‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:2578175/69‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:1471945/69‑1 (MQ=255)
                                              cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:842673/69‑1 (MQ=255)
                                              cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:1382638/69‑1 (MQ=255)
                                                 aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:2206387/66‑1 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:2838485/65‑1 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:609191/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1079137/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1117648/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:2522657/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1421546/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1183988/1‑68 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                          >  1:1355909/1‑39 (MQ=37)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:524393/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1623520/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1338005/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:3106101/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                >  1:1037444/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:1846364/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:3435284/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:2349330/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:323369/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:3085217/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:1796533/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:1869656/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:2691014/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:2729787/69‑1 (MQ=255)
                                                       ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         >  1:2600390/1‑68 (MQ=255)
                                                       ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:2900743/60‑1 (MQ=255)
                                                          gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                  >  1:917761/3‑56 (MQ=255)
                                                            tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc   <  1:1977844/69‑1 (MQ=255)
                                                             aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:3002395/1‑56 (MQ=255)
                                                             aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:2124180/59‑1 (MQ=255)
                                                             aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  >  1:716318/1‑69 (MQ=255)
                                                             aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  <  1:3274702/69‑1 (MQ=255)
                                                                 |                                                                
GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  minE/214866‑214995

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: