Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 85,885 G→A 41.3% T127T (ACG→ACA murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE85,8850GA41.3% 30.2 / 51.6 75T127T (ACG→ACAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (30/14);  new base A (22/9);  total (52/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGA  >  minE/85825‑85944
                                                            |                                                           
tcgtGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaa                                                     <  1:3073241/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:3066349/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:3212417/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:1598602/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:774810/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:1537775/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:1403721/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:1407708/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACACGGCAAAACGACa                                              <  1:3329471/69‑1 (MQ=255)
       gAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACACAAGGCAAAACGACa                                              <  1:2012412/69‑1 (MQ=255)
         aaTGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAac                                            <  1:2958649/69‑1 (MQ=255)
          aTGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                                           <  1:725070/69‑1 (MQ=255)
          aTGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                                           <  1:2059822/69‑1 (MQ=255)
          aTGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                                           <  1:1215175/69‑1 (MQ=255)
            gctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                                         >  1:800026/1‑69 (MQ=255)
            gctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                                         >  1:2123769/1‑69 (MQ=255)
            gctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                                         >  1:3211866/1‑69 (MQ=255)
            gctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                                         >  1:300477/1‑69 (MQ=255)
             ctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc                                        >  1:2937936/1‑69 (MQ=255)
             ctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc                                        >  1:19620/1‑69 (MQ=255)
             ctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc                                        >  1:2001241/1‑69 (MQ=255)
             ctggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccac                                         >  1:1960029/1‑68 (MQ=255)
               ggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc                                      >  1:2575560/1‑69 (MQ=255)
               ggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc                                      >  1:2008693/1‑69 (MQ=255)
               ggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc                                      >  1:2452970/1‑69 (MQ=255)
               ggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc                                      >  1:592620/1‑69 (MQ=255)
               ggctgAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgc                                      >  1:814556/1‑69 (MQ=255)
                      ttAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt                               >  1:1928510/1‑69 (MQ=255)
                      ttAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt                               >  1:1424436/1‑69 (MQ=255)
                           gCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAac                                            <  1:280229/51‑1 (MQ=255)
                           gCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAg                          <  1:738603/69‑1 (MQ=255)
                           gCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAg                          <  1:2886001/69‑1 (MQ=255)
                           gCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAg                          <  1:985012/69‑1 (MQ=255)
                           gCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAg                          <  1:989245/69‑1 (MQ=255)
                                           aTCGCCATTGCCGGAACGCACGGCAAAACGACAac                                            <  1:2433118/35‑1 (MQ=255)
                                              gCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCgggg       >  1:1014283/1‑69 (MQ=255)
                                              gCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCgggg       >  1:1790283/1‑69 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt                               <  1:564268/39‑1 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt                               <  1:2166519/39‑1 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt                               <  1:1220883/39‑1 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCGTCTACGCAGAAGCGGGGCTCGa  >  1:799994/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:679734/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:2804260/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1428378/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1222044/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1920492/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1286073/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:918600/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1765575/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:1633829/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:2885142/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:2369836/1‑68 (MQ=255)
                                                    gCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGa  >  1:524908/1‑68 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCa                        >  1:3276458/1‑41 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:785147/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:758368/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:1059410/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:1276388/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:524247/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:812216/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:906706/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:646090/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:1450874/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:591617/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:575014/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:1451565/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:1602947/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:3049122/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2006317/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2031961/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2899708/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2137410/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2779422/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2560408/1‑42 (MQ=255)
                                                         aaCACACGGCAAAACGACAACCACCGCGATGGTTTCCAGCAt                       >  1:2393718/1‑42 (MQ=255)
                                                            |                                                           
TCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACGCAGAAGCGGGGCTCGA  >  minE/85825‑85944

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: