Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,484 T→C 11.8% I327I (ATT→ATC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,4840TC11.8% 199.5 / 15.9 76I327I (ATT→ATCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base C (5/4);  total (42/34)
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.81e-01

CCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCG  >  minE/189414‑189554
                                                                      |                                                                      
ccAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                                                        <  1:3814791/71‑1 (MQ=255)
 cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                                                        >  1:108453/1‑70 (MQ=255)
 cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTc                                                                       <  1:3592842/71‑1 (MQ=255)
 cAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTc                                                                       <  1:1216009/71‑1 (MQ=255)
  aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCa                                                                      >  1:110506/1‑71 (MQ=255)
  aaCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCa                                                                      >  1:2619085/1‑71 (MQ=255)
   aCTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCAc                                                                     <  1:3885216/71‑1 (MQ=255)
      ccACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACtt                                                                   >  1:1136506/1‑70 (MQ=255)
       cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                                  <  1:124671/70‑1 (MQ=255)
       cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                                  <  1:1565380/70‑1 (MQ=255)
       cACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTc                                                                  >  1:1987341/1‑70 (MQ=255)
          cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAtt                                                              >  1:789431/1‑71 (MQ=255)
          cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAtt                                                              >  1:3353835/1‑71 (MQ=255)
          cTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGAt                                                               >  1:647505/1‑70 (MQ=255)
            atgatgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGa                                                                >  1:2435795/1‑67 (MQ=255)
                tgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCt                                                        <  1:2713246/71‑1 (MQ=255)
                  gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                      <  1:168469/71‑1 (MQ=255)
                  gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTg                                                       <  1:2679645/70‑1 (MQ=255)
                   cTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTggg                                                     <  1:3873571/71‑1 (MQ=255)
                   cTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTcgg                                                     <  1:2089225/71‑1 (MQ=255)
                     tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                   >  1:3499844/1‑71 (MQ=255)
                      gCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGc                                                  >  1:87838/1‑71 (MQ=255)
                        tctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGa                                                                >  1:1899268/1‑55 (MQ=255)
                          tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGcctt                                               <  1:3242697/70‑3 (MQ=255)
                          tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                              >  1:2404752/1‑71 (MQ=255)
                          tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                              >  1:3923733/1‑71 (MQ=255)
                          tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                              <  1:3537515/71‑1 (MQ=255)
                          tACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTa                                              <  1:3480662/71‑1 (MQ=255)
                          tACAACACCGACTGGATCAAACAGACTCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                <  1:1280254/69‑1 (MQ=255)
                          tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                <  1:1557950/69‑1 (MQ=255)
                          tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                <  1:1474286/69‑1 (MQ=255)
                          tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                              >  1:540775/1‑71 (MQ=255)
                          tACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                              >  1:476107/1‑71 (MQ=255)
                           acaacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGcctt                                               <  1:673216/69‑3 (MQ=255)
                           acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                              >  1:1123100/1‑70 (MQ=255)
                           acaacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                              >  1:2144761/1‑70 (MQ=255)
                            caacaCCGACTTCATCAACAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                            >  1:3742465/1‑71 (MQ=255)
                               cacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGAtt                                                                        <  1:424376/40‑1 (MQ=255)
                                acagACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCa                                        >  1:669736/4‑71 (MQ=255)
                                acCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCa                                        >  1:3811930/1‑71 (MQ=255)
                                acCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACATAGATTCt                                                            >  1:3890559/1‑51 (MQ=255)
                                 ccGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTa                                              <  1:59962/64‑1 (MQ=255)
                                   gACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAttt                                     <  1:3002159/71‑1 (MQ=255)
                                   gACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAttt                                     <  1:757113/71‑1 (MQ=255)
                                    aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt                                    >  1:3379048/1‑71 (MQ=255)
                                    aCTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt                                    >  1:2482060/1‑71 (MQ=255)
                                     cTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCAtttt                                    <  1:3861866/70‑1 (MQ=255)
                                      tGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                            <  1:1512994/61‑1 (MQ=255)
                                       ggATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                    >  1:781594/1‑52 (MQ=255)
                                        gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAcc                                          >  1:755731/1‑61 (MQ=255)
                                        gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                                  <  1:1265227/69‑1 (MQ=255)
                                        gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                                  <  1:2709764/69‑1 (MQ=255)
                                        gATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCAt                                 <  1:1294762/70‑1 (MQ=255)
                                           cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                            >  1:416130/1‑56 (MQ=255)
                                           cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                            >  1:1569762/1‑56 (MQ=255)
                                             aaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATc                           >  1:382467/1‑71 (MQ=255)
                                             aaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATc                           <  1:2252234/71‑1 (MQ=255)
                                              aCAGACGCTGGATGTCCCGTCGATTGACTTCGATTCTGCCTGGGTGCCGTacac                                           >  1:1430428/1‑54 (MQ=255)
                                              aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCt                          <  1:2792172/71‑1 (MQ=255)
                                              aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTATGCCTGGGTGCCGTACACCCATTTTCAta                                >  1:1193992/1‑64 (MQ=255)
                                               cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCg                              >  1:2710033/1‑66 (MQ=255)
                                                  aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAc                        >  1:483846/1‑69 (MQ=255)
                                                    gCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                      >  1:2573694/1‑37 (MQ=255)
                                                    gCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCa                                  <  1:872736/57‑1 (MQ=255)
                                                        gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTa                 >  1:651933/1‑70 (MQ=255)
                                                        gATGTCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaa               <  1:667737/71‑1 (MQ=255)
                                                         aTGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaa               <  1:3153832/71‑1 (MQ=255)
                                                            tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGt            <  1:2497198/71‑1 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGGCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg          >  1:3355659/1‑71 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                            >  1:683974/1‑53 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCATGGTaa                <  1:2600200/65‑1 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg          <  1:3386711/71‑1 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTgg          >  1:205432/1‑71 (MQ=255)
                                                              ccGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTg           <  1:1167861/70‑1 (MQ=255)
                                                                gTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAggg                   >  1:2070441/1‑60 (MQ=255)
                                                                gTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTa        >  1:3089411/1‑71 (MQ=255)
                                                                 tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                         >  1:1457966/1‑53 (MQ=255)
                                                                 tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaa               >  1:1462377/1‑63 (MQ=255)
                                                                 tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt       >  1:523469/1‑71 (MQ=255)
                                                                      tCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAgcg  >  1:3752459/1‑71 (MQ=255)
                                                                      |                                                                      
CCAACTCCACCTATGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCG  >  minE/189414‑189554

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