Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 87,506 G→A 10.9% L178L (TTG→TTA ddlB → D‑Alanine:D‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE87,5060GA10.9% 202.4 / 12.6 73L178L (TTG→TTAddlBD‑Alanine:D‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base A (6/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAAAGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTC  >  minE/87437‑87571
                                                                     |                                                                 
aaaaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattga                                                                  >  1:1910941/1‑71 (MQ=255)
aaaaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattg                                                                   <  1:3926107/70‑1 (MQ=255)
  aaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattgatt                                                                >  1:561298/1‑71 (MQ=255)
   agtagtAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattgatt                                                                >  1:2736590/1‑70 (MQ=255)
     tagtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaa                                                             <  1:2202706/71‑1 (MQ=255)
     tagtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaa                                                             <  1:3522590/71‑1 (MQ=255)
       gtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaaa                                                            >  1:1901573/1‑70 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt                                                         <  1:2796325/70‑1 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt                                                         <  1:3783991/70‑1 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt                                                         <  1:3653775/70‑1 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg                                                        >  1:1022684/1‑71 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg                                                        >  1:49429/1‑71 (MQ=255)
          gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg                                                        >  1:3663242/1‑71 (MQ=255)
             gAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGc                                                      >  1:496292/1‑70 (MQ=255)
              aaaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCt                                                     >  1:3328293/1‑70 (MQ=255)
              aaaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa                                                    <  1:523555/71‑1 (MQ=255)
               aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg                                                        >  1:1423266/1‑66 (MQ=255)
               aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTaa                                                   <  1:688994/71‑1 (MQ=255)
               aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa                                                    <  1:2203262/70‑1 (MQ=255)
               aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa                                                    <  1:399325/70‑1 (MQ=255)
                aaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAg                                                  >  1:2776234/1‑71 (MQ=255)
                    ctctACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTggg                                              <  1:714186/71‑1 (MQ=255)
                       tACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGc                                             <  1:1853581/69‑1 (MQ=255)
                       tACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg                                           >  1:952001/1‑71 (MQ=255)
                        aCAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaaaaa                                                          <  1:1668945/55‑1 (MQ=255)
                        aCAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg                                           >  1:2252793/1‑70 (MQ=255)
                         cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg                                           <  1:850058/69‑1 (MQ=255)
                         cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg                                           >  1:874656/1‑69 (MQ=255)
                         cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa                                         <  1:1086445/71‑1 (MQ=255)
                            gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      >  1:2212623/1‑71 (MQ=255)
                            gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      >  1:2697714/1‑71 (MQ=255)
                            gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      <  1:2061174/71‑1 (MQ=255)
                            gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      >  1:873590/1‑71 (MQ=255)
                              tGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGcc                                            >  1:1283155/1‑63 (MQ=255)
                               gCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCAc                                   >  1:3363820/1‑71 (MQ=255)
                                cATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCAc                                   <  1:1503440/70‑1 (MQ=255)
                                  ttAAGATTGGCATTTCAGCACGATGAAGAAGTAttgattga                                                              <  1:3033176/41‑1 (MQ=255)
                                   tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGcc                                            >  1:717445/1‑58 (MQ=255)
                                   tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt                               >  1:3404997/1‑71 (MQ=255)
                                   tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt                               >  1:2015352/1‑71 (MQ=255)
                                   tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt                               >  1:416369/1‑71 (MQ=255)
                                   tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGt                                <  1:724142/70‑1 (MQ=255)
                                    aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACgg                                 >  1:1952423/1‑68 (MQ=255)
                                    aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt                               >  1:755904/1‑70 (MQ=255)
                                    aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTg                              >  1:1783941/1‑71 (MQ=255)
                                    aaGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaa                                                             >  1:1389254/1‑40 (MQ=255)
                                    aaGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaa                                                             >  1:1066798/1‑40 (MQ=255)
                                      gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa                                                    >  1:169781/1‑47 (MQ=255)
                                      gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg                            <  1:1394584/71‑1 (MQ=255)
                                      gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg                            <  1:880899/71‑1 (MQ=255)
                                      gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg                            <  1:2446568/71‑1 (MQ=255)
                                       aTTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa                           <  1:3248154/71‑1 (MQ=255)
                                       aTTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa                           <  1:760291/71‑1 (MQ=255)
                                         tGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTgg                                               >  1:70090/1‑49 (MQ=255)
                                          ggCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATAc                        >  1:1439692/1‑71 (MQ=255)
                                             aTTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCg                     >  1:3496096/1‑71 (MQ=255)
                                              tttCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAATTGGGCCGGAGTTCACGGTTGCGATACTCg                     >  1:3614355/1‑70 (MQ=255)
                                               ttCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt                   >  1:1017626/1‑71 (MQ=255)
                                               ttCAGCACGATGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt                   >  1:2909913/1‑71 (MQ=255)
                                               ttCAGCACGATGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt                   >  1:1226640/1‑71 (MQ=255)
                                                tCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt                   <  1:2991558/70‑1 (MQ=255)
                                                 cAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa                                         >  1:2146666/1‑47 (MQ=255)
                                                 cAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa                                         >  1:3112522/1‑47 (MQ=255)
                                                  aGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTga                 <  1:3895381/70‑1 (MQ=255)
                                                   gCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCgg                    <  1:1948444/66‑1 (MQ=255)
                                                    cACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      <  1:2065807/47‑1 (MQ=255)
                                                    cACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTgaag               >  1:926755/1‑70 (MQ=255)
                                                     aCGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt                                      >  1:2352595/1‑46 (MQ=255)
                                                       gATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAAt           >  1:2343965/1‑71 (MQ=255)
                                                         tGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAAttt         >  1:3328773/1‑71 (MQ=255)
                                                            agaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTgaagaa             <  1:1334574/64‑1 (MQ=255)
                                                             gaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTAcc     >  1:2772154/1‑71 (MQ=255)
                                                              aagTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCg    <  1:1257842/71‑1 (MQ=255)
                                                               agTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa                           >  1:513025/1‑47 (MQ=255)
                                                                gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGt   >  1:3549608/1‑70 (MQ=255)
                                                                gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc  <  1:1529224/71‑1 (MQ=255)
                                                                gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc  <  1:466497/71‑1 (MQ=255)
                                                                gTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc  <  1:952263/71‑1 (MQ=255)
                                                                     gattgaAAAATGGCTAAGTGGGCCGGAGTTCACgg                                 >  1:1096686/1‑35 (MQ=255)
                                                                     |                                                                 
AAAAGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTC  >  minE/87437‑87571

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: