Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,580 C→T 9.4% I359I (ATC→ATT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,5800CT 9.4% 149.1 / 5.7 53I359I (ATC→ATTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base T (3/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAA  >  minE/189511‑189646
                                                                     |                                                                  
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATc                                                                    <  1:3546463/70‑1 (MQ=255)
 acCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCtt                                                                  <  1:507438/71‑1 (MQ=255)
    cATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGa                                                               <  1:779238/71‑1 (MQ=255)
     aTTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGaa                                                              <  1:3336562/71‑1 (MQ=255)
     aTTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGa                                                               >  1:367910/1‑70 (MQ=255)
     aTTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGa                                                               >  1:3424102/1‑70 (MQ=255)
     aTTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGa                                                               >  1:1892658/1‑70 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACg                                                           >  1:2187315/1‑71 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTCGTATGAGCGGCGAGCGTGATGCGGGAAAAGTGATCTTCGaa                                                              >  1:2855256/1‑68 (MQ=255)
            tCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAAc                                                            <  1:3640658/66‑1 (MQ=255)
              cGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATc                                                      >  1:1163969/1‑70 (MQ=255)
               gATTTACGAAGGTAAATGCGGTATGAGCGGTGGCCGTGTAGAAGGGAAAGTGATTTACGaaa                                                             <  1:573051/62‑1 (MQ=255)
               gATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGa                                                    >  1:1372304/1‑71 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGa                                                    >  1:1228178/1‑69 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                  >  1:895463/1‑71 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                  >  1:1263497/1‑71 (MQ=255)
                     ccAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                                         <  1:1081263/60‑1 (MQ=255)
                     ccAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                                         <  1:1621085/60‑1 (MQ=255)
                         ggTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAAt                                          >  1:1494879/1‑71 (MQ=255)
                          gTAAATGCGGTATGAGCGGTGGCCGTGTAGAAGGGAAAGTGATTTACGAAACCCa                                                         <  1:3047840/55‑1 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATg                                         >  1:2345305/1‑71 (MQ=255)
                           tAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGa                                                               <  1:2213161/48‑1 (MQ=255)
                                gTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTggcg                                   <  1:1708418/71‑1 (MQ=255)
                                      tGAGCGGTGGCCGTGTAGAAGGGAAAGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGttc                             >  1:2339175/1‑70 (MQ=255)
                                      tGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAAt                                                       >  1:3239481/1‑45 (MQ=255)
                                      tGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGtt                              <  1:2899027/70‑1 (MQ=255)
                                      tGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGtt                              <  1:1459859/70‑1 (MQ=255)
                                        aGCGGCGAGCTTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATc                           >  1:31746/1‑71 (MQ=255)
                                         gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCgt                                                    <  1:2919381/45‑2 (MQ=255)
                                         gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                          <  1:1353903/71‑1 (MQ=255)
                                           ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                        <  1:3767439/71‑1 (MQ=255)
                                           ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                        <  1:168812/71‑1 (MQ=255)
                                           ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                        <  1:3518502/71‑1 (MQ=255)
                                            gcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGCCGGCGTTATCGCa                        <  1:2304605/70‑1 (MQ=255)
                                              gAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaa                                             <  1:2937092/47‑1 (MQ=255)
                                                 cGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTc                  >  1:602715/1‑71 (MQ=255)
                                                    gtTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCtt                   <  1:3125983/67‑1 (MQ=255)
                                                     ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCTACAAAATGCTGGCGGCGTTATCGCAGGCt                    <  1:598359/65‑1 (MQ=255)
                                                     ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCt               >  1:2457588/1‑70 (MQ=255)
                                                      tGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg              <  1:555725/70‑1 (MQ=255)
                                                         gggAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGc                                        <  1:702812/41‑1 (MQ=255)
                                                         gggAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATc           >  1:875236/1‑70 (MQ=255)
                                                          gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca         >  1:3431554/4‑71 (MQ=255)
                                                          gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca         >  1:1732390/4‑71 (MQ=255)
                                                           gAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcac        <  1:1232237/71‑1 (MQ=255)
                                                            aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcaca       <  1:2900179/71‑1 (MQ=255)
                                                             aaaGTGATCTTCGAAACGCAATCGACCCACAAAAt                                          >  1:3676241/1‑35 (MQ=255)
                                                             aaaGTGATCTTCGAAACGCAATCGACCCACAAAAt                                          >  1:215213/1‑35 (MQ=255)
                                                             aaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                          <  1:2430194/51‑1 (MQ=255)
                                                              aaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACAt      >  1:3542458/1‑70 (MQ=255)
                                                               aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa    <  1:1927645/71‑1 (MQ=255)
                                                                gTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa    >  1:2282388/1‑70 (MQ=255)
                                                                   aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTaaa  <  1:3768447/69‑1 (MQ=255)
                                                                   aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa    <  1:238895/67‑1 (MQ=255)
                                                                     |                                                                  
CACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAA  >  minE/189511‑189646

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: