Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,123 C→T 37.2% G540G (GGC→GGT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1230CT37.2% 36.3 / 40.3 43G540G (GGC→GGTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base T (5/11);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCA  >  minE/190055‑190190
                                                                    |                                                                   
cTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTAt                                                                   <  1:2701062/71‑1 (MQ=255)
cTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTAt                                                                   <  1:2073456/71‑1 (MQ=255)
  cgacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                                                                  >  1:2661880/1‑70 (MQ=255)
   gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCg                                                                 <  1:3073030/70‑1 (MQ=255)
   gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGa                                                                >  1:1524719/1‑71 (MQ=255)
   gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg                                                                 <  1:498848/70‑1 (MQ=255)
    acgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                               <  1:742748/71‑1 (MQ=255)
     cgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                               >  1:2295930/1‑70 (MQ=255)
        aCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                           <  1:3567069/71‑1 (MQ=255)
           tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                                                        <  1:1405728/71‑1 (MQ=255)
            ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                               <  1:1101192/60‑1 (MQ=12)
             gcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaaa                                                      <  1:2785490/69‑1 (MQ=255)
              caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                                                     <  1:3691571/70‑1 (MQ=255)
              caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                                                     <  1:3696164/70‑1 (MQ=255)
               aTCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaa                                                       <  1:92617/68‑1 (MQ=255)
                tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                   <  1:1021879/71‑1 (MQ=255)
                tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                   >  1:2204207/1‑71 (MQ=255)
                tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                   <  1:3732692/71‑1 (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:1819542/1‑47 (MQ=255)
                                        ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCac                                                  >  1:955875/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)
                                                                    |                                                                   
CTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCA  >  minE/190055‑190190

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: