Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,582 T→A 9.6% F360Y (TTC→TAC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,5820TA 9.6% 137.8 / 7.0 52F360Y (TTC→TAC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base A (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGC  >  minE/189512‑189649
                                                                      |                                                                   
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:3518502/71‑1 (MQ=255)
                                          ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                           <  1:168812/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:215213/1‑35 (MQ=255)
                                                            aaaGTGATCTTCGAAACGCAATCGACCCACAAAAt                                             >  1:3676241/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)
                                                                     ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                 <  1:2168222/54‑1 (MQ=255)
                                                                     ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg   >  1:2955028/1‑68 (MQ=255)
                                                                     ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc  >  1:754662/1‑69 (MQ=255)
                                                                      |                                                                   
ACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGC  >  minE/189512‑189649

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: