Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 103,919 A→T 52.5% K149N (AAA→AAT ampD → N‑acetyl‑anhydromuranmyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE103,9190AT52.5% 3.5 / 92.4 78K149N (AAA→AATampDN‑acetyl‑anhydromuranmyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base T (17/24);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCA  >  minE/103853‑103976
                                                                  |                                                         
tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                         <  1:2047764/69‑1 (MQ=255)
tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                         <  1:2084648/69‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:1750643/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:3306261/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:3302241/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:1366320/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:1843084/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:3066118/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:799645/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:852018/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:254978/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:1086269/68‑1 (MQ=255)
    cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                                      <  1:2703846/68‑1 (MQ=255)
        aTCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:3025380/1‑69 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                           <  1:1841352/58‑1 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                           <  1:3333844/58‑1 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                           <  1:2691550/58‑1 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                           <  1:1502213/58‑1 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:1667056/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:822743/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:3051682/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:659340/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:2068336/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:2052387/1‑68 (MQ=255)
         tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                                 >  1:2048434/1‑68 (MQ=255)
                cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATtgt                                         <  1:447192/69‑1 (MQ=255)
                cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATtgt                                         <  1:1566049/69‑1 (MQ=255)
                  tGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGa                                       <  1:2026663/69‑1 (MQ=255)
                   gcggcgGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                            >  1:1067478/1‑63 (MQ=255)
                   gcggcgGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                            >  1:2227987/1‑63 (MQ=255)
                           tACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                         <  1:3069161/42‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTCTCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:1537309/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:983483/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:3267144/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:2756030/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:2829391/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:931658/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:396027/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:3428824/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:712750/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:412403/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:432831/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:3211105/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:3215060/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:2656179/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:1218172/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:1724932/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:1789697/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:1991449/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:207720/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:244315/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:254175/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:2663030/68‑1 (MQ=255)
                               cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:2621155/68‑1 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGTCCATTGTGATATTgcgc                               >  1:1163058/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1024783/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:938444/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1120121/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:772481/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1235776/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:424635/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1783944/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:3335886/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1788154/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:1838878/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:3009940/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:2950320/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:2550487/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:3115524/1‑59 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgc                                 >  1:2804756/1‑57 (MQ=255)
                                    aCTGATTGATCGCTATCCAGATATCACTAATAACATGACGGGCCATTGTGATATTgcgc                               >  1:3005960/1‑59 (MQ=255)
                                     cTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGaaa                    >  1:309259/1‑69 (MQ=255)
                                                 tATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCgg         <  1:233810/68‑1 (MQ=255)
                                                 tATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCgg         <  1:86038/68‑1 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGa                           >  1:2748529/1‑46 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg                >  1:1337200/1‑57 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg                >  1:129690/1‑57 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg                <  1:1282494/57‑1 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg                >  1:2593106/1‑57 (MQ=255)
                                                     cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCg                >  1:1924233/1‑57 (MQ=255)
                                                       gATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCa  <  1:1174672/69‑1 (MQ=255)
                                                       gATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCa  <  1:1156916/69‑1 (MQ=255)
                                                       gATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCa  <  1:925660/69‑1 (MQ=255)
                                                       gATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCa  <  1:1936813/69‑1 (MQ=255)
                                                          atCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                           <  1:738949/41‑1 (MQ=255)
                                                                  |                                                         
TGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCA  >  minE/103853‑103976

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: