Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,665 G→A 10.0% A388T (GCC→ACC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6650GA10.0% 143.3 / 5.6 50A388T (GCC→ACC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (3/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGTT  >  minE/189598‑189730
                                                                    |                                                                 
ccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGc                                                                  >  1:577684/1‑69 (MQ=255)
 cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATCACGAAGAGGCCtt                                                               >  1:2614385/1‑71 (MQ=255)
  acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                                              <  1:3878029/71‑1 (MQ=255)
   caAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTa                                                             <  1:106517/71‑1 (MQ=255)
   caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                                              >  1:109014/1‑70 (MQ=255)
   caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                                              >  1:342552/1‑70 (MQ=255)
   caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                                              >  1:1440890/1‑70 (MQ=255)
    aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTaa                                                            >  1:3327390/1‑71 (MQ=255)
    aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTaa                                                            >  1:3151622/1‑71 (MQ=255)
     aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAAc                                                           <  1:468205/71‑1 (MQ=255)
      aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACg                                                          <  1:1989331/71‑1 (MQ=255)
      aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACg                                                          <  1:2524304/71‑1 (MQ=255)
          cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                                                       <  1:1100066/70‑1 (MQ=255)
          cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                                                       <  1:949683/70‑1 (MQ=255)
          cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                                                       <  1:1421299/70‑1 (MQ=255)
             gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTAACGAAGCCt                                                    <  1:1810438/70‑1 (MQ=9)
                gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTat                                                >  1:1164167/1‑71 (MQ=255)
                   ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat                                             <  1:3050058/71‑1 (MQ=255)
                     aTCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                                            >  1:541136/1‑70 (MQ=255)
                      tCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGc                                           <  1:1846554/70‑1 (MQ=255)
                                cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGc                                                                  >  1:1171744/1‑37 (MQ=255)
                                cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc                                >  1:2429243/1‑71 (MQ=255)
                                cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc                                >  1:1936249/1‑71 (MQ=255)
                                 gCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCt                               >  1:1620737/1‑71 (MQ=255)
                                  atGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                              >  1:2163227/2‑71 (MQ=9)
                                     aTCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCg                             >  1:608139/1‑69 (MQ=255)
                                        cacaTTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAg                        <  1:3307876/71‑1 (MQ=255)
                                         acaTTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCa                         >  1:3013327/1‑69 (MQ=255)
                                          caTTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGtt                      <  1:2203421/71‑1 (MQ=255)
                                          caTTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGtt                      <  1:723364/71‑1 (MQ=255)
                                          caTTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGt                       >  1:14578/1‑70 (MQ=255)
                                            ttAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGc                                                      >  1:3027248/1‑37 (MQ=255)
                                              aaaGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATa                                        <  1:3889274/49‑1 (MQ=255)
                                               aaGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatg                                               <  1:207303/41‑1 (MQ=255)
                                               aaGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatg                                               <  1:204087/41‑1 (MQ=255)
                                               aaGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATccc                 <  1:563082/71‑1 (MQ=255)
                                                aGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATacca                                     <  1:3345295/50‑1 (MQ=255)
                                                  gCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGcc                           <  1:2642849/58‑1 (MQ=255)
                                                   cGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCTCCACCTCGCCCa                         <  1:3931271/59‑1 (MQ=255)
                                                   cGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg             <  1:3083419/71‑1 (MQ=255)
                                                   cGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg             <  1:485159/71‑1 (MQ=255)
                                                    gAGTATGACGATGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCg                             <  1:3155015/54‑1 (MQ=255)
                                                    gAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttgt            <  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)
                                                               aaGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGtt  <  1:990962/71‑1 (MQ=255)
                                                                 gagGCCTTTAAGGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCg     <  1:2411128/66‑1 (MQ=255)
                                                                 gagGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCg     <  1:90137/66‑1 (MQ=255)
                                                                    gCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCgg    >  1:3629703/1‑64 (MQ=255)
                                                                    |                                                                 
CCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGTT  >  minE/189598‑189730

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: