Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 85,795 C→A 36.8% A97A (GCC→GCA murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE85,7950CA36.8% 69.2 / 50.6 73A97A (GCC→GCAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base A (20/7);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGT  >  minE/85735‑85861
                                                            |                                                                  
tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATg                                                             <  1:1475059/68‑1 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:2657213/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:418154/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:1895214/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:2461887/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:679783/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:2465818/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                        >  1:1139901/1‑68 (MQ=255)
     ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgc                                                         >  1:716253/1‑67 (MQ=255)
       aGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGcgcg                                                      <  1:914139/68‑1 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2113430/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:3279938/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2063991/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:1950353/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:1912017/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2980834/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2205583/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:1893919/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:1849214/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2966344/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:3324169/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2851824/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2833015/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:652277/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:668186/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:1503255/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:2744764/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:926644/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc                                                                >  1:270735/1‑51 (MQ=255)
              tCGTTGTTTCCAGCGCGATTTCTGCCCATAACCCGGAAATTGACGCAGCTc                                                                >  1:1114910/1‑51 (MQ=38)
               cGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                             >  1:339875/1‑69 (MQ=255)
               cGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                             >  1:582912/1‑69 (MQ=255)
                gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                            >  1:1047129/1‑69 (MQ=255)
                gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                            >  1:1731210/1‑69 (MQ=255)
                gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                            >  1:1550106/1‑69 (MQ=255)
                  tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                           <  1:829224/68‑1 (MQ=255)
                  tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                           <  1:966003/68‑1 (MQ=255)
                  tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                           <  1:1045018/68‑1 (MQ=255)
                   gttTCCAGCGCTATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:2251874/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:47038/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:97616/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:2635868/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:1252370/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:1828749/56‑1 (MQ=255)
                   gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg                                                      <  1:2065945/56‑1 (MQ=255)
                      tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                      >  1:1341358/1‑69 (MQ=255)
                      tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                      >  1:2470052/1‑69 (MQ=255)
                      tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                      >  1:2193091/1‑69 (MQ=255)
                      tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcg                                       >  1:1099762/1‑68 (MQ=255)
                      tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGCCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt                                      >  1:2812007/1‑69 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGCTCCGTCGTGCCGAACTGGTGGCTGa                  >  1:62575/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:2453500/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:2601393/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:1315032/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:2735061/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:3245296/1‑68 (MQ=255)
                                           aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa                  >  1:2239244/1‑68 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2464190/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2881622/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:728799/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2700594/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2696828/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1358890/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1766462/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1764246/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:517736/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1518461/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2882793/68‑1 (MQ=255)
                                              ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1520824/68‑1 (MQ=255)
                                                   aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:607783/63‑1 (MQ=255)
                                                   aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:1177481/63‑1 (MQ=255)
                                                   aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt               <  1:2883544/63‑1 (MQ=255)
                                                     ttGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtt       <  1:226387/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:1661053/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:182589/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:3353200/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:2887549/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:2593624/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:3068612/69‑1 (MQ=255)
                                                          gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt  <  1:1151842/69‑1 (MQ=255)
                                                            |                                                                  
TGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGT  >  minE/85735‑85861

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: