Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,656 G→A 9.0% D385N (GAC→AAC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6560GA 9.0% 162.1 / 5.4 56D385N (GAC→AAC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base A (3/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

AACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTT  >  minE/189586‑189725
                                                                       |                                                                     
aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATg                                                                       <  1:2130320/71‑1 (MQ=255)
 aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                                      >  1:2772430/1‑71 (MQ=255)
 aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                                                      >  1:1912756/1‑71 (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:3327390/1‑71 (MQ=255)
                aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTaa                                                       >  1:3151622/1‑71 (MQ=255)
                aaaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGa                                                                   >  1:3602599/1‑59 (MQ=255)
                 aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                                                      <  1:468205/71‑1 (MQ=255)
                  aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                                     <  1:1989331/71‑1 (MQ=255)
                  aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                                     <  1:2524304/71‑1 (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:1195601/45‑1 (MQ=255)
                               ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                                                                  <  1:2776494/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)
                                                        ttAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGc                                                 >  1:3027248/1‑37 (MQ=255)
                                                          aaaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATa                                   <  1:3889274/49‑1 (MQ=255)
                                                           aaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatg                                          <  1:207303/41‑1 (MQ=255)
                                                           aaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatg                                          <  1:204087/41‑1 (MQ=255)
                                                           aaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATccc            <  1:563082/71‑1 (MQ=255)
                                                            aGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATacca                                <  1:3345295/50‑1 (MQ=255)
                                                              gCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGcc                      <  1:2642849/58‑1 (MQ=255)
                                                               cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCTCCACCTCGCCCa                    <  1:3931271/59‑1 (MQ=255)
                                                               cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg        <  1:3083419/71‑1 (MQ=255)
                                                               cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg        <  1:485159/71‑1 (MQ=255)
                                                                gA‑GTATGACGATGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCg                        <  1:3155015/54‑1 (MQ=255)
                                                                gA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttgt       <  1:2170908/71‑1 (MQ=255)
                                                                   gTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttgttg     >  1:425362/1‑71 (MQ=255)
                                                                      tGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCtt  >  1:1524118/1‑71 (MQ=255)
                                                                      tGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCtt  >  1:822975/1‑71 (MQ=255)
                                                                      tGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCt   >  1:818377/1‑70 (MQ=255)
                                                                       |                                                                     
AACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTT  >  minE/189586‑189725

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