Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 103,899 T→C 51.3% C143R (TGC→CGC)  ampD → N‑acetyl‑anhydromuranmyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE103,8990TC51.3% 1.9 / 67.7 78C143R (TGC→CGC) ampDN‑acetyl‑anhydromuranmyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base T (12/26);  new base C (17/23);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAA  >  minE/103833‑103958
                                                                  |                                                           
gATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGc                                                           >  1:2684806/1‑69 (MQ=255)
      cTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCg                                                     <  1:1073327/69‑1 (MQ=255)
      cTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCg                                                     <  1:3031815/69‑1 (MQ=255)
                ccGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTa                                                         <  1:973610/55‑1 (MQ=255)
                ccGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTa                                                         <  1:2346060/55‑1 (MQ=255)
                ccGATGCGAAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTa                                                         <  1:562738/55‑1 (MQ=255)
                    tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                       <  1:2084648/69‑1 (MQ=255)
                    tGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaaaa                                       <  1:2047764/69‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:799645/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:254978/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:2703846/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:1086269/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:852018/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:1843084/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:3066118/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:3302241/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:1750643/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:1366320/68‑1 (MQ=255)
                        cAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACAt                                    <  1:3306261/68‑1 (MQ=255)
                            aTCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:3025380/1‑69 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                         <  1:3333844/58‑1 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                         <  1:1841352/58‑1 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                         <  1:2691550/58‑1 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTaaa                                         <  1:1502213/58‑1 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:2052387/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:822743/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:659340/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:2068336/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:2048434/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:1667056/1‑68 (MQ=255)
                             tCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACgg                               >  1:3051682/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: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:412403/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:3211105/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:3215060/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:2621155/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:1789697/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:3428824/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:396027/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:3267144/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:432831/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:1724932/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:712750/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:1218172/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:931658/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:983483/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:2656179/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:1991449/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:207720/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:2663030/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:2829391/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:254175/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:2756030/68‑1 (MQ=255)
                                                   cgGGCACTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTGCGCCGGa         <  1:244315/68‑1 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGTCCATTGTGATATTgcgc             >  1:1163058/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:3115524/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:2550487/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:2950320/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:3335886/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:1783944/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:1788154/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:3009940/1‑59 (MQ=255)
                                                        aCTGATTGATCGCTATCCAGATATCGCTAATAACATGACGGGCCATTGTGATATTgcgc             >  1:1838878/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)
                                                                  |                                                           
GATACGCTGGCGTATACCGATGCGCAGTATCAACAGCTTGCGGCGGTTACGCGGGCACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAA  >  minE/103833‑103958

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