Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,736,011 T→A 50.9% A385A (GCT→GCA fusA → protein chain elongation factor EF‑G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,736,0110TA50.9% ‑4.7 / 63.9 53A385A (GCT→GCAfusAprotein chain elongation factor EF‑G
Reads supporting (aligned to +/- strand):  ref base T (15/11);  new base A (14/13);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CAGATGCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTGG  >  minE/2735949‑2736075
                                                              |                                                                
cAGATGCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGt                                                            >  1:1009222/1‑69 (MQ=255)
cAGATGCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGt                                                            >  1:2574099/1‑69 (MQ=255)
cAGATGCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGt                                                            >  1:888432/1‑69 (MQ=255)
     gCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGa                                                        <  1:1321748/68‑1 (MQ=255)
           tAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGt                                                 <  1:1898897/69‑1 (MQ=255)
            aaCAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTa                                                >  1:3312875/1‑69 (MQ=255)
                  cGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACt                                          <  1:1583513/69‑1 (MQ=255)
                  cGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACt                                          <  1:1259291/69‑1 (MQ=255)
                  cGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACt                                          <  1:3187745/69‑1 (MQ=255)
                  cGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACt                                          <  1:1412047/69‑1 (MQ=255)
                     gAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGt                                       <  1:3248031/69‑1 (MQ=255)
                     gAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGt                                       <  1:2573010/69‑1 (MQ=255)
                     gAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGt                                       <  1:2007716/69‑1 (MQ=255)
                     gAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGt                                       <  1:2759570/69‑1 (MQ=255)
                     gAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGGCTGAAAGACGTAACCACTGGt                                       <  1:3155562/69‑1 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:2028124/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:2548755/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1969725/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:2155741/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:193716/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:243705/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:2445469/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1030640/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1074228/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1076641/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1487777/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:555407/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:1728818/1‑68 (MQ=255)
                      aaGAGATCAAAGAAGTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGt                                       >  1:776713/1‑68 (MQ=255)
                                agaagTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAAc                                              >  1:840819/1‑51 (MQ=255)
                                agaagTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAAc                                              >  1:2662036/1‑51 (MQ=255)
                                agaagTGCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAAc                                              >  1:440949/1‑51 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:2917511/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:942725/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:576203/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:426018/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:340028/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:333841/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:3235527/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:2493129/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:1240973/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:2171980/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:2494438/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:1913059/69‑1 (MQ=255)
                                 gaagTTCGCGCGGGCGACATCGCAGCTGCAATCGGTCTGAAAGACGTAACCACTGGTGACACCCtgtgt                           <  1:2114198/69‑1 (MQ=255)
                                    gTTCGCGCTGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGAc                        >  1:3369025/1‑69 (MQ=255)
                                    gTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGAc                        >  1:2896035/1‑69 (MQ=255)
                                    gTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGAc                        >  1:1305174/1‑69 (MQ=255)
                                                   aTCGCTGCTGCTATCGGTCTGAAAGACGTACCCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATc         >  1:2372305/1‑69 (MQ=255)
                                                   aTCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATc         >  1:301641/1‑69 (MQ=255)
                                                   aTCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATc         >  1:1193926/1‑69 (MQ=255)
                                                    tCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCa        >  1:106995/1‑69 (MQ=255)
                                                    tCGCTGCTGCTATCGGTCTGAAAGACGTAAACACTGGTGACACCCTGTGTGACCCGGATGCGCCGATc         >  1:1573040/1‑68 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGAGGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:1255360/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:354465/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:890817/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:1521353/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:2902614/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:3369394/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:2857480/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGAt          <  1:3292399/65‑1 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:376904/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:140969/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:525050/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:165380/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:1760966/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:2094315/1‑69 (MQ=255)
                                                      gctgctgctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCAtt      >  1:1004957/1‑69 (MQ=255)
                                                            gctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTgg  <  1:34886/67‑1 (MQ=255)
                                                            gctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTgg  <  1:59521/67‑1 (MQ=255)
                                                            gctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTgg  <  1:184366/67‑1 (MQ=255)
                                                            gctATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTgg  <  1:1973457/67‑1 (MQ=255)
                                                              |                                                                
CAGATGCACGCTAACAAACGTGAAGAGATCAAAGAAGTTCGCGCGGGCGACATCGCTGCTGCTATCGGTCTGAAAGACGTAACCACTGGTGACACCCTGTGTGACCCGGATGCGCCGATCATTCTGG  >  minE/2735949‑2736075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: