Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,665 A→C 33.6% V443V (GTA→GTC thrA → fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,6650AC33.6% 104.7 / 108.7 134V443V (GTA→GTCthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base A (42/47);  new base C (11/34);  total (53/81)
Fisher's exact test for biased strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 139 total aligned reads displayed.

AATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTA  >  minE/1606‑1732
                                                           |                                                                   
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:1834485/69‑1 (MQ=255)
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:2851375/69‑1 (MQ=255)
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:2773508/69‑1 (MQ=255)
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:949031/69‑1 (MQ=255)
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:2617138/69‑1 (MQ=255)
aaTATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:20619/69‑1 (MQ=255)
 atatCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:1307525/68‑1 (MQ=255)
    tCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgatgat                                                         >  1:1235010/1‑68 (MQ=255)
    tCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgatgat                                                         >  1:2122549/1‑68 (MQ=255)
    tCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgatgat                                                         >  1:1935487/1‑68 (MQ=255)
    tCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:2847486/65‑1 (MQ=255)
      aaCATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCg                                                      >  1:1496799/1‑69 (MQ=255)
            gTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACt                                                >  1:1047345/1‑69 (MQ=255)
            gTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACt                                                >  1:383642/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGt                                           >  1:567089/1‑68 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGTg                                          >  1:2703972/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGTg                                          >  1:2990284/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGTg                                          >  1:572579/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGTg                                          >  1:855227/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGATGATGCGACCACTGGCGTg                                          >  1:1166596/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAAAGATGATGCGACCACTGGCGTg                                          >  1:2229966/1‑69 (MQ=255)
                  aTTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTg                                          <  1:3049758/69‑1 (MQ=255)
                    tGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgcg                                        >  1:2433766/1‑69 (MQ=255)
                    tGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgc                                         >  1:1807028/1‑68 (MQ=255)
                     gCTGAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:1586014/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:3296027/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:1255095/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:2660480/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:2180813/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:1988892/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:2722066/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:3364516/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:103353/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:65092/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            >  1:2854018/1‑48 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACga                                                             >  1:1814442/1‑47 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgcgc                                       >  1:1455621/1‑69 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgcgc                                       >  1:1940917/1‑69 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTgcgc                                       >  1:281943/1‑69 (MQ=255)
                     gCTCAGGGATCTTCTGAACGCTCAATCTCAGTCGTGGAAAATAACgat                                                            >  1:2420158/1‑48 (MQ=255)
                           ggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTAc                                  >  1:2008719/1‑68 (MQ=255)
                           ggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTAc                                  >  1:1823622/1‑68 (MQ=255)
                           ggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACt                                 >  1:449772/1‑69 (MQ=255)
                             aTCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACgat                                                            <  1:1224980/40‑1 (MQ=255)
                               cttctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtc                                >  1:2455374/1‑66 (MQ=255)
                               cttctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtc                                >  1:2264882/1‑66 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:162776/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:3376067/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:633818/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:546164/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:1908995/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:459952/65‑1 (MQ=255)
                                    gAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            <  1:925908/65‑1 (MQ=255)
                                          tCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTc                  <  1:3393605/69‑1 (MQ=255)
                                          tCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTc                  <  1:115593/69‑1 (MQ=255)
                                          tCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTc                  <  1:3355718/69‑1 (MQ=255)
                                           cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACtcatca                            >  1:3137445/1‑58 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:377914/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:3275058/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:921385/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:3194633/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:3278714/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:824169/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:512137/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:69195/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:952489/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2783189/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1680813/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1765801/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1517803/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1965131/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2054166/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1324065/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1256542/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2224376/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2366380/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2347995/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1171658/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2587970/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:2799519/67‑1 (MQ=255)
                                                      gTGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:268125/67‑1 (MQ=255)
                                                      gTGGTCAATAAAGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:156797/67‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:725030/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:3147618/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:1780050/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:555975/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:1783972/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:727680/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:3423095/68‑1 (MQ=255)
                                                      gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa       <  1:923257/68‑1 (MQ=255)
                                                       tGGTCAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATc        <  1:1532401/66‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:1078727/69‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:1834394/69‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:1948063/69‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:2205636/69‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:2547604/69‑1 (MQ=255)
                                                         gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGtt   <  1:2634018/69‑1 (MQ=255)
                                                           aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTa  <  1:139926/68‑1 (MQ=255)
                                                           aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTa  <  1:207223/68‑1 (MQ=255)
                                                           aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTa  <  1:3378461/68‑1 (MQ=255)
                                                           |                                                                   
AATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTA  >  minE/1606‑1732

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: