Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,534,296 G→A 10.4% A147T (GCA→ACA)  yfcN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,534,2960GA10.4% 195.7 / 11.4 67A147T (GCA→ACA) yfcNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/36);  new base A (3/4);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TTTGCGCCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTG  >  minE/1534231‑1534364
                                                                 |                                                                    
tttGCGCCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGACacaa                                                                 >  1:1193087/1‑71 (MQ=255)
 ttGCGCCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCacaa                                                                 >  1:3323498/1‑70 (MQ=255)
  tGCGCCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCacaac                                                                <  1:1196773/70‑1 (MQ=255)
   gcgcCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCacaaca                                                               <  1:1221960/70‑1 (MQ=255)
   gcgcCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCacaaca                                                               <  1:2493083/70‑1 (MQ=255)
         gCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCAt                                                        >  1:73622/1‑71 (MQ=255)
            tGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCAtgt                                                      <  1:1698570/70‑1 (MQ=255)
                gCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGtggctggc                                                                     <  1:1386806/51‑1 (MQ=255)
                gCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATgg                                                 >  1:1744499/1‑71 (MQ=255)
                 cATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGc                                                          <  1:3500981/61‑1 (MQ=255)
                 cATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGc                                                          <  1:2331876/61‑1 (MQ=255)
                 cATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATgg                                                 <  1:2781242/70‑1 (MQ=255)
                  aTGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGa                                                    >  1:799037/1‑66 (MQ=255)
                   tGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATgg                                                 <  1:3026759/68‑1 (MQ=255)
                    gggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCtt                                             >  1:1447977/1‑71 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATcc                                                            <  1:2463289/55‑1 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATc                                                             >  1:2011672/1‑54 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGa                                                    >  1:3129456/1‑63 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCttt                                            >  1:3607884/1‑71 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCttt                                            >  1:133178/1‑71 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCttt                                            >  1:1813580/1‑71 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCAAACATCCGCATGTGATGGCCttt                                            >  1:2580322/1‑70 (MQ=255)
                     ggCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGACACAACATCCGCATGTGATGGCCttt                                            >  1:3450016/1‑71 (MQ=255)
                      gCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCg                                                           <  1:2285177/55‑1 (MQ=255)
                      gCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCttt                                            >  1:3172975/1‑70 (MQ=255)
                       cACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTtca                                          <  1:3171487/71‑1 (MQ=255)
                              aGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCTCATGTGATGGCCTTTCATCAg                                      <  1:2577082/68‑1 (MQ=255)
                              aGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAg                                      <  1:3649562/68‑1 (MQ=255)
                               gCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCAtgt                                                      <  1:721647/51‑1 (MQ=255)
                                cATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCt                                              <  1:248920/58‑1 (MQ=255)
                                  tatTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGc                                                          <  1:1741204/44‑1 (MQ=255)
                                   atTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaa                               <  1:2545427/70‑1 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACAt                                                              <  1:1970247/38‑1 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACAt                                                              <  1:973662/38‑1 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGa                                                    >  1:932970/1‑48 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGa                                                    >  1:973468/1‑48 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCtt                                             <  1:2373616/55‑1 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaa                               >  1:3342539/1‑69 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGGCACAACATACGCATGTGATGGCCTTTCATCAGGCACCaa                               >  1:656150/1‑69 (MQ=255)
                                    ttttGAAGCAACAAACACCGCTGTGGCTGACACAACATCCGCATGTGAt                                                   >  1:641767/1‑49 (MQ=255)
                                      ttGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaaaa                             >  1:4143992/1‑69 (MQ=255)
                                        gAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAagag                          <  1:3055737/70‑1 (MQ=255)
                                        gAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAagag                          <  1:515477/70‑1 (MQ=255)
                                        gAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAagag                          >  1:1611070/1‑70 (MQ=255)
                                        gAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAagag                          <  1:240111/70‑1 (MQ=255)
                                           gCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATg                      <  1:1838431/71‑1 (MQ=255)
                                            caacaaACACCGCTGTGGCTGGCACAACATCCGGATGTGATGGCTTTtcat                                         <  1:532607/51‑1 (MQ=255)
                                             aacaaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCtt                                             <  1:1653070/46‑1 (MQ=255)
                                              acaaacaCCGCTGTGGCTGACACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCg                   <  1:2396392/71‑1 (MQ=255)
                                               caaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATgg                                                 <  1:2459454/40‑1 (MQ=255)
                                               caaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTtca                                          >  1:1533525/1‑47 (MQ=255)
                                                aaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAgg                                     <  1:3939236/51‑1 (MQ=255)
                                                aaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:1366032/71‑1 (MQ=255)
                                                aaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:1734249/71‑1 (MQ=255)
                                                aaacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:1137641/71‑1 (MQ=255)
                                                aaacaCCGCTGTGGCTGACACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:3655426/71‑1 (MQ=255)
                                                 aacaCCGCTGTGGCTGGCACAACATCCGCATGTGa                                                    >  1:3162234/1‑35 (MQ=255)
                                                 aacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:2550365/70‑1 (MQ=255)
                                                 aacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 >  1:3359799/1‑70 (MQ=255)
                                                 aacaCCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:2555241/70‑1 (MQ=255)
                                                 aacaCCGCTGTGGCGGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGt                 <  1:2656430/70‑1 (MQ=255)
                                                   cacCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGt                         >  1:130118/1‑60 (MQ=255)
                                                   cacCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTAt                       <  1:3548498/62‑1 (MQ=255)
                                                   cacCGCTGTGGCTGACACAACATCCGCATGTGATGGCCTTtca                                          <  1:2438490/43‑1 (MQ=255)
                                                     ccGCTGTGGCTGGCACAACATCCGCATGTGATGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATgcg           <  1:1917853/71‑1 (MQ=255)
                                                       gCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATgcg           >  1:2972632/1‑70 (MQ=255)
                                                         tgtgGCTGACACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTAt                       <  1:2581895/56‑1 (MQ=255)
                                                         tgtgCCTGGCACAACATCCCCATGTGATGGCCTTTCATCAGGc                                    >  1:2557681/1‑43 (MQ=255)
                                                           tggctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaaaa                             <  1:2577951/48‑1 (MQ=255)
                                                           tggctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaaaa                             <  1:2048178/48‑1 (MQ=255)
                                                           tggctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCgtt      >  1:1691801/1‑71 (MQ=255)
                                                            ggctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCaa                               <  1:101428/45‑1 (MQ=255)
                                                             gctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTa                        >  1:219338/1‑51 (MQ=255)
                                                              ctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCgttgtt   <  1:3553834/71‑1 (MQ=255)
                                                              ctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCgttgtt   <  1:3225999/71‑1 (MQ=255)
                                                              ctggcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCgttgtt   <  1:255174/71‑1 (MQ=255)
                                                                 gcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTa                        >  1:98955/1‑47 (MQ=255)
                                                                 gcACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGttgttg  <  1:328163/69‑1 (MQ=255)
                                                                 |                                                                    
TTTGCGCCTGCGTGATGCATGGGCACGGGAAGCATATTTTGAAGCAACAAACACCGCTGTGGCTGGCACAACATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTG  >  minE/1534231‑1534364

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: