Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,078,085 A→G 39.3% F144F (TTT→TTC)  uxaA ← altronate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,078,0850AG39.3% 22.4 / 36.5 56F144F (TTT→TTC) uxaAaltronate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (11/23);  new base G (5/17);  total (16/40)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAG  >  minE/2078028‑2078147
                                                         |                                                              
ggTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTg                                                     >  1:2567235/1‑69 (MQ=255)
  tgtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                                    <  1:2947421/68‑1 (MQ=255)
   gtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                                    <  1:709345/67‑1 (MQ=255)
              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                                       >  1:2084722/1‑69 (MQ=255)
              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                                       >  1:3329225/1‑69 (MQ=255)
              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                                       >  1:3159527/1‑69 (MQ=255)
              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                                       >  1:1852938/1‑69 (MQ=255)
              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                                       >  1:1685572/1‑69 (MQ=255)
               gAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                                    <  1:1841839/55‑1 (MQ=255)
                 acacGCGGTCTGTACCTTCAGCGTTGTTGGTCTCTTTCTGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:3044091/69‑1 (MQ=255)
                 acacGCCGTCGGTGCCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:1798262/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:1422785/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:2737380/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:867023/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:248017/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:3058572/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:728283/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:1039965/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:2067543/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:1862837/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:2997281/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:818916/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:1727676/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:2901926/69‑1 (MQ=255)
                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTCGATCTGCCGCGCGATGc                                    <  1:180504/69‑1 (MQ=255)
                  cacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                                    <  1:604527/68‑1 (MQ=255)
                  cacGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGcc                                   >  1:3428281/1‑69 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTg                                               >  1:955038/1‑54 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCg                                  >  1:998354/1‑67 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt                                >  1:1506129/1‑69 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt                                >  1:1663247/1‑69 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt                                 >  1:846910/1‑68 (MQ=255)
                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt                                 >  1:1911962/1‑68 (MQ=255)
                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa                              <  1:1505307/69‑1 (MQ=255)
                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa                              <  1:2281934/69‑1 (MQ=255)
                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa                              <  1:3360240/69‑1 (MQ=255)
                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa                              <  1:835379/69‑1 (MQ=255)
                           aGTACGTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:278673/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:407657/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:274742/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:3164099/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:2696541/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:2612671/54‑1 (MQ=255)
                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:1114633/54‑1 (MQ=255)
                           aGTACCTTCGGCGTT‑TTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                                         <  1:759826/53‑1 (MQ=255)
                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc                      >  1:3149105/1‑69 (MQ=255)
                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc                      >  1:939995/1‑69 (MQ=255)
                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc                      >  1:1162840/1‑69 (MQ=255)
                                       gttgttGGTCTCTTTCAGAAAACGGTGCTGGATCTGCCGCGCGATGCCGTTGa                              <  1:610352/53‑1 (MQ=255)
                                        ttgttgGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGc             <  1:225874/69‑1 (MQ=255)
                                            tgGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAgg          <  1:1778056/68‑1 (MQ=255)
                                            tgGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAgg          <  1:1641539/68‑1 (MQ=255)
                                                     tCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAg  <  1:3188625/67‑1 (MQ=255)
                                                     tCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAg  <  1:449716/67‑1 (MQ=255)
                                                     tCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAg  <  1:1717310/67‑1 (MQ=255)
                                                     tCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAg  <  1:1425472/67‑1 (MQ=255)
                                                         |                                                              
GGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCCACGGTTGGCAGGATCCACAG  >  minE/2078028‑2078147

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: ‑