Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,641 A→G 9.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,6410AG 9.2% 146.5 / 5.2 54I380V (ATT→GTT) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (26/23);  new base G (3/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTT  >  minE/189571‑189710
                                                                      |                                                                      
aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcaca                                                                        <  1:2900179/71‑1 (MQ=255)
  aaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACAt                                                                       >  1:3542458/1‑70 (MQ=255)
   aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                     <  1:1927645/71‑1 (MQ=255)
    gTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                     >  1:2282388/1‑70 (MQ=255)
       aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTaaa                                                                   <  1:3768447/69‑1 (MQ=255)
       aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                                                     <  1:238895/67‑1 (MQ=255)
          ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                                                 >  1:2955028/1‑68 (MQ=255)
          ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                                >  1:754662/1‑69 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTTATCCACATTAAAGGc                                                                <  1:3926652/67‑1 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                                <  1:1068284/67‑1 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                                                <  1:375758/67‑1 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                               <  1:1811161/68‑1 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                            >  1:1049595/1‑70 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                            <  1:2076571/70‑1 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                            >  1:3454258/1‑70 (MQ=255)
            cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTAGCTGATCCACATTAAAGGCg                                                               <  1:2447914/68‑1 (MQ=255)
               aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATg                                                        <  1:2130320/71‑1 (MQ=255)
                aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                       >  1:1912756/1‑71 (MQ=255)
                aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                       >  1:2772430/1‑71 (MQ=255)
                 cGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                               <  1:3825564/63‑1 (MQ=255)
                     aTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACg                                                     >  1:3144947/1‑68 (MQ=255)
                      tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                   >  1:30969/1‑69 (MQ=255)
                      tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                   >  1:1795071/1‑69 (MQ=255)
                           ccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGc                                              >  1:577684/1‑69 (MQ=255)
                            cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATCACGAAGAGGCCtt                                           >  1:2614385/1‑71 (MQ=255)
                             acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                          <  1:3878029/71‑1 (MQ=255)
                              caAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTa                                         <  1:106517/71‑1 (MQ=255)
                              caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                          >  1:109014/1‑70 (MQ=255)
                              caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                          >  1:1440890/1‑70 (MQ=255)
                              caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                          >  1:342552/1‑70 (MQ=255)
                               aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTaa                                        >  1:3151622/1‑71 (MQ=255)
                               aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTaa                                        >  1:3327390/1‑71 (MQ=255)
                               aaaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGa                                                    >  1:3602599/1‑59 (MQ=255)
                                aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                                       <  1:468205/71‑1 (MQ=255)
                                 aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                      <  1:2524304/71‑1 (MQ=255)
                                 aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                      <  1:1989331/71‑1 (MQ=255)
                                  aTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTat                                                                    <  1:2838230/41‑2 (MQ=255)
                                   tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                               >  1:489712/1‑45 (MQ=255)
                                   tGCTGGCGGCGTTATCGCAGGCTTCGCTGATACACATTAAAGGCGA‑GTATGACGAag                                                  >  1:874904/1‑57 (MQ=255)
                                     cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg                                   <  1:949683/70‑1 (MQ=255)
                                     cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg                                   <  1:1100066/70‑1 (MQ=255)
                                     cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg                                   <  1:1421299/70‑1 (MQ=255)
                                        gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCt                                <  1:1810438/70‑1 (MQ=9)
                                           gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                            >  1:1164167/1‑71 (MQ=255)
                                              ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                   <  1:2776494/45‑1 (MQ=255)
                                              ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                   <  1:1195601/45‑1 (MQ=255)
                                              ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat                         <  1:3050058/71‑1 (MQ=255)
                                                aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                        >  1:541136/1‑70 (MQ=255)
                                                 tCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGc                       <  1:1846554/70‑1 (MQ=255)
                                                           cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGc                                              >  1:1171744/1‑37 (MQ=255)
                                                           cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc            >  1:2429243/1‑71 (MQ=255)
                                                           cGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc            >  1:1936249/1‑71 (MQ=255)
                                                            gCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCt           >  1:1620737/1‑71 (MQ=255)
                                                             atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc          >  1:2163227/2‑71 (MQ=9)
                                                                aTCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCg         >  1:608139/1‑69 (MQ=255)
                                                                   cacaTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAg    <  1:3307876/71‑1 (MQ=255)
                                                                    acaTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCa     >  1:3013327/1‑69 (MQ=255)
                                                                     caTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGtt  <  1:723364/71‑1 (MQ=255)
                                                                     caTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGtt  <  1:2203421/71‑1 (MQ=255)
                                                                     caTTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGt   >  1:14578/1‑70 (MQ=255)
                                                                      |                                                                      
AAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTT  >  minE/189571‑189710

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