Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 87,241 C→T 49.3% A90V (GCA→GTA)  ddlB → D‑Alanine:D‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE87,2410CT49.3% 2.5 / 70.2 79A90V (GCA→GTA) ddlBD‑Alanine:D‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (12/28);  new base T (25/14);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 3.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-03

TGAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGTGCCGGTTTACCGGTCGCG  >  minE/87173‑87308
                                                                    |                                                                   
tGAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGc                                                                     <  1:177858/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGGTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1430592/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGGTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1435043/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1177004/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:2777175/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:2896194/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:3358109/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1565270/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:680739/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1441580/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1412773/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:711635/69‑1 (MQ=255)
 gAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:1289172/69‑1 (MQ=255)
          aCGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCt                                                            <  1:2660932/68‑1 (MQ=255)
          aCGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCt                                                            <  1:2661635/68‑1 (MQ=255)
          aCGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCt                                                            <  1:909188/68‑1 (MQ=255)
          aCGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCt                                                            <  1:564881/68‑1 (MQ=255)
          aCGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCt                                                            <  1:770273/68‑1 (MQ=255)
                    agATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATa                                                 <  1:1549176/68‑1 (MQ=255)
                     gATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATa                                                 <  1:2914712/68‑1 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCttt                                                          >  1:1211423/1‑57 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATaaa                                               >  1:1752219/1‑68 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:3427018/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:362277/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:50780/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:585048/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:671072/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2876063/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2830071/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2783119/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:983592/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2739871/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:737484/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:68132/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2464804/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2354780/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1477273/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1554133/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1562011/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1758501/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1885709/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:1917600/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2023137/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2177574/1‑69 (MQ=255)
                       tGCTCGAGCTGATGGGCTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAAc                                              >  1:2145397/1‑69 (MQ=255)
                         cTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTa                                            <  1:2295534/69‑1 (MQ=255)
                         cTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTa                                            <  1:218619/69‑1 (MQ=255)
                            gAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCa                                                                    <  1:2574131/42‑1 (MQ=255)
                                 gATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTc                                                         >  1:1148637/1‑48 (MQ=255)
                                 gATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTc                                                         >  1:212941/1‑48 (MQ=255)
                                 gATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTc                                                         >  1:2298665/1‑48 (MQ=255)
                                 gATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTc                                                         >  1:2493064/1‑48 (MQ=255)
                                 gATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTc                                                         >  1:2779171/1‑48 (MQ=255)
                                     ggCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACt                                 >  1:9178/1‑68 (MQ=255)
                                     ggCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACt                                 >  1:257707/1‑68 (MQ=255)
                                     ggCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACt                                 >  1:2905800/1‑68 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:1609926/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:890946/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:872006/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:1394410/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:1463684/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:2276055/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:611836/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:2400540/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:460410/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:2587483/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:2222270/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:2676939/69‑1 (MQ=255)
                                       cTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTc                               <  1:1691438/68‑1 (MQ=255)
                                       cTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTc                               <  1:144027/68‑1 (MQ=255)
                                       cTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTc                               <  1:2809981/68‑1 (MQ=255)
                                        ttGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTc                               <  1:1673509/67‑1 (MQ=255)
                                             ttATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGc                         <  1:1551411/68‑1 (MQ=255)
                                             ttATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGc                         <  1:1271221/68‑1 (MQ=255)
                                               ataCCGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:1629282/61‑1 (MQ=255)
                                                 aCCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGt                    >  1:3428459/1‑69 (MQ=255)
                                                   cGGCAGCGGAGTGATGGTATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCt                              <  1:1349629/57‑1 (MQ=255)
                                                                   gCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGTGCCGGTTTACCGGTcgcg  >  1:2237331/1‑69 (MQ=255)
                                                                   gCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGTGCCGGTTTACCGGTcgcg  >  1:311967/1‑69 (MQ=255)
                                                                   gCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGTGCCGGTTTACCGGTcgcg  >  1:2865439/1‑69 (MQ=255)
                                                                    |                                                                   
TGAAGATGGTACGCTGCAGGGGATGCTCGAGCTGATGGGCTTGCCTTATACCGGAAGCGGAGTGATGGCATCTGCGCTTTCAATGGATAAACTACGCAGCAAACTTCTATGGCAAGGTGCCGGTTTACCGGTCGCG  >  minE/87173‑87308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: