Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,132 A→G 25.6% K543K (AAA→AAG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1320AG25.6% 59.3 / 22.1 39K543K (AAA→AAGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (14/15);  new base G (4/6);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGAT  >  minE/190063‑190200
                                                                     |                                                                    
aCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                                     <  1:3567069/71‑1 (MQ=255)
   tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                                                                  <  1:1405728/71‑1 (MQ=255)
     gcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaaa                                                                <  1:2785490/69‑1 (MQ=255)
      caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                                                               <  1:3696164/70‑1 (MQ=255)
      caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                                                               <  1:3691571/70‑1 (MQ=255)
       aTCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaa                                                                 <  1:92617/68‑1 (MQ=255)
        tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                             <  1:1021879/71‑1 (MQ=255)
        tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                             <  1:3732692/71‑1 (MQ=255)
        tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                             >  1:2204207/1‑71 (MQ=255)
             gttgagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                                                 >  1:3158362/4‑62 (MQ=255)
             gtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATgg                                                         <  1:720342/70‑1 (MQ=255)
             gtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATgg                                                         <  1:1373617/70‑1 (MQ=255)
                gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAAt                                                           >  1:2195073/1‑65 (MQ=255)
                   aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                                              <  1:1483645/59‑1 (MQ=255)
                   aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATgg                                                         <  1:3110597/64‑1 (MQ=255)
                       ccGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGt                                              <  1:2901634/71‑1 (MQ=255)
                        cGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                                             >  1:3067743/1‑71 (MQ=255)
                            ccTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGtt                                         >  1:2788236/1‑71 (MQ=255)
                            ccTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGtt                                         >  1:2851227/1‑71 (MQ=255)
                                ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatta                                                    <  1:2108835/56‑1 (MQ=255)
                                ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                                     >  1:1323502/1‑71 (MQ=255)
                                ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGACTATTGCGTGGGTTGACg                                     >  1:588582/1‑71 (MQ=255)
                                ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCac                                                            >  1:955875/1‑47 (MQ=255)
                                ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCac                                                            >  1:1819542/1‑47 (MQ=255)
                                   aCCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAAt                                                           >  1:3657983/1‑46 (MQ=255)
                                     ctgctgTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACgg                                    <  1:1493677/67‑1 (MQ=255)
                                           tttCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaaa                            <  1:2156852/69‑1 (MQ=255)
                                             tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGct                        >  1:2558698/1‑71 (MQ=255)
                                                 tttAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAAGGCTCtt                     >  1:1746693/1‑70 (MQ=255)
                                                   tAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTAc                   >  1:220243/1‑70 (MQ=255)
                                                   tAGTATTGGCATAGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTAc                   >  1:1208391/1‑70 (MQ=255)
                                                    aGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTc                               <  1:1201368/57‑1 (MQ=255)
                                                        ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtct              >  1:2883999/1‑70 (MQ=255)
                                                          ggCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCa            <  1:3641069/70‑1 (MQ=255)
                                                             aTCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCt        >  1:2271686/1‑71 (MQ=255)
                                                              tCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGct                        <  1:3290706/54‑1 (MQ=255)
                                                                gATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaa                             >  1:3387435/1‑47 (MQ=255)
                                                                gATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCa                              <  1:3887518/46‑1 (MQ=255)
                                                                gATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCg     <  1:3364433/71‑1 (MQ=255)
                                                                 aTAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCgg    <  1:2420863/71‑1 (MQ=255)
                                                                   aaaaCCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGAt  >  1:3441137/1‑71 (MQ=255)
                                                                    aaaCCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCgg    >  1:315715/1‑68 (MQ=255)
                                                                    aaaCCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGAt  >  1:1165228/1‑70 (MQ=255)
                                                                     |                                                                    
ACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGCGGAT  >  minE/190063‑190200

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: