Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 103,907 G→A 50.0% P145P (CCG→CCA ampD → N‑acetyl‑anhydromuranmyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE103,9070GA50.0% ‑2.3 / 72.9 80P145P (CCG→CCAampDN‑acetyl‑anhydromuranmyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base A (17/23);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

CTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGG  >  minE/103839‑103969
                                                                    |                                                              
cTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCg                                                                <  1:1073327/69‑1 (MQ=255)
cTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCg                                                                <  1:3031815/69‑1 (MQ=255)
              tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                  <  1:2084648/69‑1 (MQ=255)
              tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                  <  1:2047764/69‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:2703846/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:3302241/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:3306261/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:1843084/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:254978/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:3066118/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:1750643/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:1366320/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:1086269/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:852018/68‑1 (MQ=255)
                  cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                               <  1:799645/68‑1 (MQ=255)
                      aTCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:3025380/1‑69 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                    <  1:1502213/58‑1 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                    <  1:2691550/58‑1 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                    <  1:1841352/58‑1 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                                    <  1:3333844/58‑1 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:659340/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:3051682/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:822743/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:1667056/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:2048434/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:2052387/1‑68 (MQ=255)
                       tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                                          >  1:2068336/1‑68 (MQ=255)
                              cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATtgt                                  <  1:1566049/69‑1 (MQ=255)
                              cTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATtgt                                  <  1:447192/69‑1 (MQ=255)
                                tGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGa                                <  1:2026663/69‑1 (MQ=255)
                                 gcggcgGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                     >  1:2227987/1‑63 (MQ=255)
                                 gcggcgGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCAt                                     >  1:1067478/1‑63 (MQ=255)
                                         tACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                                  <  1:3069161/42‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTCTCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:1537309/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:3267144/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:3428824/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:3215060/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:3211105/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:396027/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:412403/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:432831/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:712750/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:2829391/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:931658/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:983483/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:2756030/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:254175/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:1724932/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:1218172/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:2621155/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:2656179/68‑1 (MQ=255)
                                             cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa                    <  1:2663030/68‑1 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGTCCATTGTGATATTgcgc                        >  1:1163058/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1235776/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1838878/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:772481/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1120121/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:424635/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:938444/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1783944/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1788154/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:1024783/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:3335886/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:3115524/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:2550487/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:3009940/1‑59 (MQ=255)
                                                  aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc                        >  1:2950320/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:1924233/1‑57 (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)
                                                                    |                                                              
CTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGG  >  minE/103839‑103969

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: