Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,641 A→G 11.2% I380V (ATT→GTT)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6410AG11.2% 218.0 / 17.5 80I380V (ATT→GTT) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (5/4);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGT  >  minE/189574‑189709
                                                                   |                                                                     
aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:252468/71‑1 (MQ=255)
aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:3511725/71‑1 (MQ=255)
aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:2083622/71‑1 (MQ=255)
 gTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTaa                                                                   <  1:3872704/71‑1 (MQ=255)
  tGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTaaa                                                                  <  1:3674234/71‑1 (MQ=255)
  tGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    >  1:3374996/1‑69 (MQ=255)
      cTTCGCAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                >  1:4006338/1‑69 (MQ=255)
         cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                               <  1:952189/67‑1 (MQ=255)
         cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                           >  1:2787942/1‑70 (MQ=255)
         cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                           <  1:3536201/70‑1 (MQ=255)
          gAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑Gt                                                          <  1:2764191/70‑1 (MQ=255)
           aaaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTAt                                                        <  1:3308721/71‑1 (MQ=255)
           aaaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTAt                                                        <  1:3790587/71‑1 (MQ=255)
            aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                               <  1:273207/64‑1 (MQ=255)
             aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                      >  1:2706938/1‑71 (MQ=255)
              cGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACAt                                                                      <  1:2146877/55‑1 (MQ=255)
              cGCAATCGACCCACAAAATGCTGGCGGCGTTATAGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGAc                                                     >  1:2076044/1‑71 (MQ=255)
                 aaTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAGGGCGA‑GTATGACGa                                                   <  1:2699099/70‑1 (MQ=255)
                 aaTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                  <  1:3639271/71‑1 (MQ=255)
                   tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:1676378/52‑1 (MQ=255)
                   tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                           <  1:272144/60‑1 (MQ=255)
                    cGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                           >  1:3371089/1‑59 (MQ=255)
                    cGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                               >  1:3630981/1‑71 (MQ=255)
                     gACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                               <  1:1569402/55‑1 (MQ=255)
                     gACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                               <  1:2897260/70‑1 (MQ=255)
                     gACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAaga                                                >  1:736254/1‑69 (MQ=255)
                       cccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                >  1:1754145/1‑52 (MQ=255)
                       cccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGcc                                            >  1:766665/1‑71 (MQ=255)
                        tcacaAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTaa                                                                   <  1:1792850/47‑1 (MQ=16)
                         cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTaaa                                                                  <  1:3768577/48‑1 (MQ=255)
                         cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:3719236/46‑1 (MQ=255)
                         cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                    <  1:3724012/46‑1 (MQ=255)
                         cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCtt                                          >  1:2948334/1‑71 (MQ=255)
                         cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCtt                                          <  1:1941925/71‑1 (MQ=255)
                          acaAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAg                                                                 >  1:2122352/1‑48 (MQ=19)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcaca                                                                       >  1:252244/1‑41 (MQ=255)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                >  1:363569/1‑48 (MQ=255)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                >  1:769940/1‑48 (MQ=255)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                >  1:3121588/1‑48 (MQ=255)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                               <  1:2507303/49‑1 (MQ=255)
                           caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTa                                                         >  1:3151360/1‑54 (MQ=255)
                            aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAg                                                                 >  1:3548272/1‑46 (MQ=19)
                            aaaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAaga                                                >  1:3413046/1‑62 (MQ=255)
                             aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTaa                                                                   >  1:1479120/1‑43 (MQ=255)
                             aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                                      <  1:547807/71‑1 (MQ=255)
                             aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                                      >  1:3269386/1‑71 (MQ=255)
                              aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑TTATGACGAAGAGGCCTTTAACg                                     <  1:407097/71‑1 (MQ=255)
                              aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                     <  1:120559/71‑1 (MQ=255)
                                tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                  <  1:2777367/56‑1 (MQ=255)
                                tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                  <  1:2417675/56‑1 (MQ=255)
                                tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                  <  1:3868561/56‑1 (MQ=255)
                                 gCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAg                                  >  1:3504986/1‑71 (MQ=12)
                                     gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCta                              <  1:2406988/71‑2 (MQ=9)
                                     gcggcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGcc                                <  1:1372940/69‑1 (MQ=255)
                                       ggcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCttt                             <  1:651137/70‑1 (MQ=255)
                                        gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                           <  1:3953913/71‑1 (MQ=255)
                                        gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                           <  1:883548/71‑1 (MQ=255)
                                        gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                           <  1:3881690/71‑1 (MQ=255)
                                             aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                               >  1:945348/1‑46 (MQ=255)
                                             aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCt                                           >  1:1069056/1‑50 (MQ=255)
                                             aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                                      >  1:3932462/1‑55 (MQ=255)
                                             aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                       <  1:2064230/70‑1 (MQ=255)
                                                 cAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTa                            <  1:2103034/61‑1 (MQ=255)
                                                  aGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATac                  >  1:181881/1‑70 (MQ=255)
                                                      ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGaa                                   >  1:3954594/1‑49 (MQ=255)
                                                      ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                           <  1:1969773/57‑1 (MQ=255)
                                                      ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATa                   >  1:2795776/1‑65 (MQ=255)
                                                        cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCt                               >  1:3852598/1‑51 (MQ=255)
                                                        cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATagcacca             <  1:2899388/69‑1 (MQ=255)
                                                        cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc           >  1:2828082/1‑71 (MQ=255)
                                                        cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccac               >  1:2562854/1‑67 (MQ=255)
                                                         gCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGc                                 >  1:3038132/1‑48 (MQ=255)
                                                         gCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGc                                 >  1:4004722/1‑48 (MQ=255)
                                                          atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc         >  1:1223342/2‑71 (MQ=9)
                                                           tGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTct        >  1:3372426/1‑70 (MQ=9)
                                                           tGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTaa                                       >  1:1317893/1‑40 (MQ=255)
                                                           tGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccac            >  1:105536/1‑67 (MQ=255)
                                                           tGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccac            >  1:1971569/1‑67 (MQ=255)
                                                            gATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTaa                                       >  1:95492/1‑39 (MQ=255)
                                                              tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                         >  1:2660156/1‑35 (MQ=255)
                                                               ccACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatga                         >  1:1816185/1‑50 (MQ=255)
                                                                 acaTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGt  <  1:2056686/71‑1 (MQ=255)
                                                                  caTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGt  >  1:197777/1‑70 (MQ=255)
                                                                   |                                                                     
AGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGT  >  minE/189574‑189709

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: ‑