Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,649 C→T 14.5% G382G (GGC→GGT)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6490CT14.5% 143.0 / 10.4 55G382G (GGC→GGT) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (3/5);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTG  >  minE/189583‑189719
                                                                  |                                                                       
cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                                        <  1:4468663/68‑1 (MQ=255)
cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                                     >  1:2828260/1‑70 (MQ=255)
cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                                     >  1:2389013/1‑70 (MQ=255)
cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATAGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                                     <  1:3524215/70‑1 (MQ=255)
 gAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTa                                                                   >  1:4136647/1‑71 (MQ=255)
   aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATg                                                                 <  1:3577523/71‑1 (MQ=255)
    aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                                >  1:4396090/1‑71 (MQ=255)
      gCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGAc                                                               >  1:1814148/1‑70 (MQ=255)
         aTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACg                                                              >  1:3528236/1‑68 (MQ=255)
            gACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAgg                                                        >  1:969445/1‑71 (MQ=255)
                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                                                     <  1:4881224/53‑1 (MQ=255)
                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            >  1:4640438/1‑62 (MQ=255)
                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                                   <  1:5283223/71‑1 (MQ=255)
                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                                   <  1:1393737/71‑1 (MQ=255)
                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCtt                                                    >  1:165013/1‑70 (MQ=255)
                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACg                                                              <  1:5583990/59‑1 (MQ=255)
                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTa                                                  <  1:5347006/71‑1 (MQ=255)
                  caAAATGCTGGCGGCGGTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGAc                                                               <  1:3759399/58‑1 (MQ=255)
                    aactGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAAc                                                <  1:3267072/68‑1 (MQ=255)
                    aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACg                                                              <  1:4168114/57‑1 (MQ=255)
                    aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            >  1:4537899/1‑59 (MQ=255)
                     actGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAAc                                                >  1:2664385/3‑70 (MQ=255)
                      aTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGa                                              <  1:1993556/71‑1 (MQ=255)
                          tGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGcc                                          >  1:5052938/1‑71 (MQ=9)
                            gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCta                                        <  1:3187422/71‑2 (MQ=9)
                            gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCta                                        <  1:1625746/71‑2 (MQ=9)
                             cggcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCttt                                       <  1:3467769/71‑1 (MQ=255)
                                 gTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatga                                   <  1:1687361/71‑1 (MQ=255)
                                    aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            >  1:374156/1‑43 (MQ=255)
                                    aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                                 >  1:5344576/1‑70 (MQ=255)
                                    aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                                 >  1:200255/1‑70 (MQ=255)
                                     tCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            <  1:676933/42‑1 (MQ=255)
                                     tCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            <  1:2421135/42‑1 (MQ=255)
                                     tCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                            <  1:4873787/42‑1 (MQ=255)
                                       gCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAaga                                                          >  1:3115107/1‑42 (MQ=255)
                                       gCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAaga                                                          >  1:952053/1‑42 (MQ=255)
                                       gCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                                     >  1:2311735/1‑63 (MQ=255)
                                         aGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                                         >  1:810037/1‑41 (MQ=255)
                                          ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                               >  1:4697611/1‑50 (MQ=255)
                                          ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGa                                              >  1:3355102/1‑51 (MQ=255)
                                          ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATacca                          >  1:5404280/1‑71 (MQ=255)
                                             ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                                         >  1:260674/1‑37 (MQ=255)
                                             ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATac                            >  1:1271987/1‑66 (MQ=255)
                                                gCTGATCCACAATAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc                     <  1:4619047/70‑1 (MQ=255)
                                                 cTGATCCACATTAAAGGCGA‑TTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat                                  <  1:224533/56‑1 (MQ=255)
                                                 atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                   >  1:3772016/2‑71 (MQ=9)
                                                 atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                   <  1:4482478/70‑1 (MQ=9)
                                                 atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                   <  1:4488805/70‑1 (MQ=9)
                                                     tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGccc               <  1:5266652/71‑1 (MQ=255)
                                                     tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGccc               <  1:3725837/71‑1 (MQ=255)
                                                     tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGcc                <  1:1177215/70‑1 (MQ=255)
                                                      ccACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatga                                   >  1:4814425/1‑50 (MQ=255)
                                                           ttAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGc                                           <  1:5376660/37‑1 (MQ=255)
                                                             aaaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATc        >  1:5594564/1‑70 (MQ=255)
                                                                 gCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAtt   >  1:3416724/1‑71 (MQ=255)
                                                                  cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg  <  1:3962792/71‑1 (MQ=255)
                                                                  cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAtt   <  1:3617773/70‑1 (MQ=255)
                                                                  |                                                                       
CGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTG  >  minE/189583‑189719

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