Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,619 A→C 9.9% L372F (TTA→TTC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6190AC 9.9% 203.2 / 19.9 81L372F (TTA→TTC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (25/48);  new base C (3/5);  total (28/53)
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.03e-01

AGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATG  >  minE/189551‑189688
                                                                    |                                                                      
aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATc                                                                      >  1:205172/1‑71 (MQ=255)
 gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:2460097/71‑1 (MQ=255)
 gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:2736570/71‑1 (MQ=255)
 gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:2194378/71‑1 (MQ=255)
      gAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGc                                                                <  1:76774/71‑1 (MQ=255)
       aGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCt                                                               <  1:3730107/71‑1 (MQ=255)
        gCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTTGCGGCGTTATCGCAGGCtt                                                              <  1:724395/71‑1 (MQ=255)
         cGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTc                                                             <  1:1918075/71‑1 (MQ=255)
             ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCt                                                          >  1:2914073/1‑70 (MQ=255)
             ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                         <  1:880447/71‑1 (MQ=255)
                cGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGTCGTTATCGCAGGCTTCGCTGATc                                                      <  1:73316/71‑1 (MQ=255)
                  ggAAAAGTGATCTTCGAAACGCAATCGACTCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca                                                    <  1:1438213/71‑1 (MQ=255)
                  ggAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca                                                    <  1:530234/71‑1 (MQ=255)
                    aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                                                                   >  1:1606991/1‑54 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTAt                                                                       <  1:1609684/47‑1 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:252056/49‑1 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               <  1:252468/71‑1 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               <  1:2083622/71‑1 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               <  1:3511725/71‑1 (MQ=255)
                        gTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTaa                                              <  1:3872704/71‑1 (MQ=255)
                         tGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTc                                                             <  1:194028/55‑1 (MQ=255)
                         tGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTaaa                                             <  1:3674234/71‑1 (MQ=255)
                         tGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               >  1:3374996/1‑69 (MQ=255)
                             cTTCGCAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                           >  1:4006338/1‑69 (MQ=255)
                               tCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                         <  1:4125853/53‑1 (MQ=255)
                               tCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                         <  1:480638/53‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCtt                                                              >  1:644382/1‑47 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGAt                                                       >  1:2475113/1‑54 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca                                                    >  1:3113847/1‑57 (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)
                                  aaaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCtt                                                              <  1:1859091/45‑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:3121588/1‑48 (MQ=255)
                                                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                           >  1:769940/1‑48 (MQ=255)
                                                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                           >  1:363569/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:3269386/1‑71 (MQ=255)
                                                    aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                 <  1:547807/71‑1 (MQ=255)
                                                     aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑TTATGACGAAGAGGCCTTTAACg                <  1:407097/71‑1 (MQ=255)
                                                     aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                <  1:120559/71‑1 (MQ=255)
                                                       tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             <  1:3868561/56‑1 (MQ=255)
                                                       tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             <  1:2417675/56‑1 (MQ=255)
                                                       tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             <  1:2777367/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)
                                                                    |                                                                      
AGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATG  >  minE/189551‑189688

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