Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 393 T→G 30.0% R19R (CGT→CGG thrA → fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE3930TG30.0% 80.1 / 51.9 80R19R (CGT→CGGthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base G (17/7);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 4.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCC  >  minE/331‑455
                                                              |                                                              
acaacCATGCGAGTGTTTAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          <  1:2464267/69‑1 (MQ=255)
acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          >  1:3197801/1‑69 (MQ=255)
acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          <  1:2820599/69‑1 (MQ=255)
acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          <  1:2231170/69‑1 (MQ=255)
acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          >  1:1718633/1‑69 (MQ=255)
acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          <  1:897886/69‑1 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:2055195/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:2816110/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:2685420/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:255094/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:3398961/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:2088381/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:1662511/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:1043343/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                                       >  1:1189324/1‑69 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                                        >  1:693059/1‑68 (MQ=255)
   acCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                                        >  1:3247282/1‑68 (MQ=255)
    ccATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGata                                                      <  1:1692370/69‑1 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:3147334/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:726509/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1049558/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:683499/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:3367740/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1103086/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1120464/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1185001/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:2207995/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:873864/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:2088957/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:3407117/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1837327/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1732993/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1319638/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:172681/1‑69 (MQ=255)
       tGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTc                                                   >  1:1602234/1‑69 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:966389/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:1313981/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:310282/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:3328688/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:2927767/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:1263356/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:2806566/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:3351718/1‑68 (MQ=255)
              gtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                                             >  1:1182336/1‑68 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:498559/69‑1 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:854887/69‑1 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:2917990/69‑1 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:363729/69‑1 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:1394436/69‑1 (MQ=255)
               ttGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                                           <  1:1291060/69‑1 (MQ=255)
                    gTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                                      <  1:1099375/69‑1 (MQ=255)
                    gTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                                      <  1:1014484/69‑1 (MQ=255)
                    gTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                                      <  1:2322819/69‑1 (MQ=255)
                    gTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                                      <  1:2980747/69‑1 (MQ=255)
                        ggcggTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                                          <  1:3417115/45‑1 (MQ=255)
                              aCATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                              >  1:2267480/1‑67 (MQ=255)
                              aCATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg                              >  1:2270019/1‑67 (MQ=255)
                                 tCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag                         <  1:832816/69‑1 (MQ=255)
                                 tCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag                         <  1:1064736/69‑1 (MQ=255)
                                      ggCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                    <  1:1778234/69‑1 (MQ=255)
                                      ggCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                    <  1:1235706/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTTTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:2796997/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:3071549/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:3048433/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:2950463/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:1443762/69‑1 (MQ=255)
                                                gAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:677171/69‑1 (MQ=255)
                                                 aaCGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:991818/68‑1 (MQ=255)
                                                  aCGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctc        >  1:1297782/1‑69 (MQ=255)
                                                   cGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc          <  1:876310/66‑1 (MQ=255)
                                                       ttCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:58587/69‑1 (MQ=255)
                                                       ttCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:2927332/69‑1 (MQ=255)
                                                       ttCTGCGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:733242/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCTGGTGGCCACCGTCCTCTCTGccc   <  1:2729012/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:788058/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:2805918/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:2176890/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:1525567/69‑1 (MQ=255)
                                                       ttCTGCGGGTTGCAGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccc   <  1:3236306/69‑1 (MQ=255)
                                                         cTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGcccc  <  1:1224052/68‑1 (MQ=255)
                                                              |                                                              
ACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCC  >  minE/331‑455

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: