Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,154 C→T 15.2% D606D (GAC→GAT thrA → fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,1540CT15.2% 152.8 / 16.5 59D606D (GAC→GATthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base T (5/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTG  >  minE/2091‑2223
                                                               |                                                                     
gATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAAc                                                                 <  1:3592978/70‑1 (MQ=255)
gATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAAc                                                                 <  1:2485792/70‑1 (MQ=255)
gATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAAc                                                                 <  1:3832260/70‑1 (MQ=255)
gATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACg                                                                >  1:1953850/1‑71 (MQ=255)
  tGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGtt                                                              >  1:3787763/1‑71 (MQ=255)
   ggATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTg                                                             >  1:2091643/1‑71 (MQ=255)
    gATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                            >  1:823320/1‑71 (MQ=255)
     aTTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                           >  1:3871754/1‑71 (MQ=255)
      ttACTAACATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggc                                                           >  1:1169609/1‑69 (MQ=255)
        actacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGc                                                         <  1:3518061/70‑1 (MQ=255)
         ctacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGtt                                                              >  1:770098/1‑64 (MQ=255)
             cATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAtt                                                   <  1:3508382/71‑1 (MQ=255)
              aTCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTa                                                  <  1:3560165/71‑1 (MQ=255)
               tCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAc                                                 >  1:3164118/1‑71 (MQ=255)
               tCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAc                                                 >  1:157310/1‑71 (MQ=255)
                 aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAAc                                                                 >  1:3645946/1‑53 (MQ=255)
                 aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTgg                                                            >  1:187416/1‑58 (MQ=255)
                   ttGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGt                                             <  1:1375647/71‑1 (MQ=255)
                   ttGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGt                                             <  1:1786112/71‑1 (MQ=255)
                   ttGCGTTATGCGGCGGAAAAATCGCGGAGTAAATTCCTCTATGATACCAACGtt                                                              <  1:310234/54‑1 (MQ=255)
                     gCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTg                                                             >  1:645579/1‑53 (MQ=255)
                     gCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGtt                                            <  1:370552/70‑1 (MQ=255)
                      cGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTa                                           <  1:1074582/70‑1 (MQ=255)
                      cGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTa                                           <  1:3158151/70‑1 (MQ=255)
                      cGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                           <  1:718565/70‑1 (MQ=255)
                      cGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                           <  1:4107384/70‑1 (MQ=255)
                      cGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                           <  1:2217498/70‑1 (MQ=255)
                       gTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTg                                                       <  1:1430542/57‑1 (MQ=255)
                          aTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGt                                                               >  1:1422376/1‑46 (MQ=255)
                          aTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                         >  1:2782438/1‑68 (MQ=255)
                             cggcggAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                   <  1:3395948/71‑1 (MQ=255)
                                  gAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                           <  1:2741405/42‑1 (MQ=255)
                                    aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTg                                                       >  1:118466/1‑44 (MQ=255)
                                    aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAt                                          >  1:3406229/1‑57 (MQ=255)
                                       aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAt                                                    >  1:2389549/1‑44 (MQ=255)
                                       aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTAcc                                                >  1:1517953/1‑48 (MQ=255)
                                       aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTg                                        <  1:1808120/56‑1 (MQ=255)
                                       aTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATc                         >  1:1937706/1‑71 (MQ=255)
                                         cgcgGCGTAAATTCCTCTATGACACCAACGTTGAGGCTGGATTa                                                  <  1:1292493/44‑1 (MQ=255)
                                          gcggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTg                       <  1:101289/70‑1 (MQ=255)
                                          gcggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTg                       <  1:1404568/70‑1 (MQ=255)
                                          gcggcgTAAATTCCTATATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                     <  1:1323130/56‑1 (MQ=255)
                                            ggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCt                     >  1:1762078/1‑70 (MQ=255)
                                              cgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                            <  1:2565751/61‑1 (MQ=255)
                                              cgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCa                   <  1:3726896/70‑1 (MQ=255)
                                                  aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaaa                           >  1:1944064/1‑58 (MQ=255)
                                                  aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCt                        >  1:1438985/1‑61 (MQ=255)
                                                  aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAAt                 >  1:2350764/1‑68 (MQ=255)
                                                  aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATg                >  1:1719154/1‑69 (MQ=255)
                                                  aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCa              <  1:587608/71‑1 (MQ=255)
                                                    ttCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                            <  1:1888937/55‑1 (MQ=255)
                                                     tCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCaa                  >  1:1533789/1‑64 (MQ=255)
                                                     tCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt           >  1:2874432/1‑71 (MQ=255)
                                                     tCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt           >  1:2295963/1‑71 (MQ=255)
                                                       ctctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCa                   >  1:745267/1‑61 (MQ=255)
                                                          tATGACACCAACGTTGGGGCTGGATTACCGGTTATTgag                                      <  1:1952936/39‑1 (MQ=255)
                                                          tATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                   >  1:1039317/1‑42 (MQ=255)
                                                            tGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg       <  1:624053/68‑1 (MQ=255)
                                                            tGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg       <  1:1259292/68‑1 (MQ=255)
                                                              acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAtt   >  1:3912392/1‑70 (MQ=255)
                                                              acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAtt   <  1:1740662/70‑1 (MQ=255)
                                                              acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAtt   >  1:1062668/1‑70 (MQ=255)
                                                              acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAtt   <  1:2580172/70‑1 (MQ=255)
                                                              acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATtg  <  1:1738052/71‑1 (MQ=255)
                                                               cacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg       <  1:512690/65‑1 (MQ=255)
                                                               |                                                                     
GATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTG  >  minE/2091‑2223

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: