Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,493 T→C 10.4% D330D (GAT→GAC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,4930TC10.4% 173.0 / 9.7 67D330D (GAT→GACldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (5/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

ATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCG  >  minE/189426‑189561
                                                                   |                                                                     
atgatgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:4954634/1‑68 (MQ=255)
 tgatgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTg                                                                   >  1:5642568/1‑71 (MQ=255)
   atgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGcc                                                                 >  1:1373105/1‑71 (MQ=255)
    tgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCt                                                                <  1:4418760/71‑1 (MQ=255)
      gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTg                                                               <  1:614753/70‑1 (MQ=255)
        ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            >  1:2674363/1‑71 (MQ=255)
        ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            >  1:2428289/1‑71 (MQ=255)
         tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:69485/1‑71 (MQ=255)
         tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:636499/1‑71 (MQ=255)
         tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:130997/1‑71 (MQ=255)
          gCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGc                                                          >  1:2076418/1‑71 (MQ=255)
           ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                         <  1:5310417/71‑1 (MQ=255)
           ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                         >  1:2328671/1‑71 (MQ=255)
             ctACAACACCGACTGGATCAAACAGCCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                        >  1:2405386/1‑70 (MQ=255)
                caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                    >  1:19933/1‑71 (MQ=255)
                caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                    >  1:2258878/1‑71 (MQ=255)
                caacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                    >  1:485429/1‑71 (MQ=255)
                 aacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                   >  1:2438796/1‑71 (MQ=255)
                  acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:3466722/61‑1 (MQ=255)
                  acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:5215632/61‑1 (MQ=255)
                  acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:302618/61‑1 (MQ=255)
                  acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:4547776/61‑1 (MQ=255)
                   cacCGACTGGACCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAccc                                                 <  1:5215005/71‑1 (MQ=255)
                    acCGTCTGGATCAAAGAGACGCTGGATGTCCCGTCGATTGACTTCGATTGTGCCTGGGTGCCGTACACCCa                                                <  1:5465383/71‑1 (MQ=255)
                     ccGACTGGAACAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGc                                                                  >  1:770251/1‑52 (MQ=255)
                      cGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTAGATTCTGCCTGGGTGCCGTACACCCAtt                                             >  1:7307/1‑71 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:2722826/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:591502/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:1169679/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:4881135/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:1960778/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:853719/1‑44 (MQ=255)
                        aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                                       >  1:3646242/1‑44 (MQ=255)
                               cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGa                                    >  1:3387964/1‑71 (MQ=255)
                                aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:801927/1‑45 (MQ=255)
                                aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:4588438/1‑45 (MQ=255)
                                aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:2494008/1‑45 (MQ=255)
                                aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:2431789/1‑45 (MQ=255)
                                  aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                   >  1:2049515/1‑54 (MQ=255)
                                  aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                    >  1:4164848/1‑53 (MQ=255)
                                   cAGACGCTGTATGTCCCGTCGATTCACTTCTATTCTg                                                                   <  1:5132272/37‑1 (MQ=38)
                                   cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                <  1:3199739/71‑1 (MQ=255)
                                   cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                >  1:1274597/1‑71 (MQ=255)
                                      aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:4213022/1‑65 (MQ=255)
                                       cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:1872445/70‑1 (MQ=255)
                                       cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                            >  1:317984/1‑71 (MQ=255)
                                            gATGTCCCGTCGATTCACTTTGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:5687877/71‑1 (MQ=255)
                                            gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                       <  1:3289155/71‑1 (MQ=255)
                                            gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                       <  1:5057101/71‑1 (MQ=255)
                                            gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:4871286/71‑1 (MQ=255)
                                              tGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaaa                     <  1:862260/71‑1 (MQ=255)
                                                tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAg                    <  1:1131259/70‑1 (MQ=255)
                                                   cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:684431/58‑1 (MQ=255)
                                                   cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:4289751/58‑1 (MQ=255)
                                                   cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:4283153/58‑1 (MQ=255)
                                                     tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:1856686/55‑1 (MQ=255)
                                                     tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:4846197/55‑1 (MQ=255)
                                                     tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:359900/55‑1 (MQ=255)
                                                       gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTa               >  1:4615643/1‑68 (MQ=255)
                                                       gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt              >  1:2351111/1‑69 (MQ=255)
                                                       gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt              >  1:792084/1‑69 (MQ=255)
                                                       gATTCACTTCGATTCTGCCTAGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGa            <  1:4226539/71‑1 (MQ=255)
                                                          tCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAgcg         >  1:2957342/1‑71 (MQ=255)
                                                           cACTTTGACTCCGCGTTGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcg         <  1:2040140/70‑1 (MQ=255)
                                                           cACTTTGACTCCGCGTGGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcgg        <  1:609115/71‑1 (MQ=255)
                                                            aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:5018399/55‑1 (MQ=255)
                                                            aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:1966890/55‑1 (MQ=255)
                                                             cTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAc                               <  1:3561015/46‑1 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:3369274/1‑39 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:2358286/1‑39 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:2805639/45‑1 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                           >  1:4457879/1‑47 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                           >  1:5009936/1‑47 (MQ=255)
                                                                cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                            >  1:4993203/1‑46 (MQ=255)
                                                                 gATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCg  <  1:5686929/71‑1 (MQ=255)
                                                                   |                                                                     
ATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCG  >  minE/189426‑189561

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: