Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,971,889 A→G 73.7% I339I (ATT→ATC speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,971,8890AG73.7% 78.9 / 23.4 68I339I (ATT→ATCspeAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (13/5);  new base G (21/29);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 5.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATA  >  minE/1971823‑1971957
                                                                  |                                                                    
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:275733/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2720676/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:1075123/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2828910/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2947501/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2542734/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2988097/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2128698/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:3168578/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:361707/68‑1 (MQ=255)
agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:1681679/68‑1 (MQ=255)
 gATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:982124/68‑1 (MQ=255)
 gATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                                                    <  1:2637126/68‑1 (MQ=255)
    tCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAAtgttg                                                                >  1:1436253/1‑69 (MQ=255)
    tCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAAtgtt                                                                 >  1:3333425/1‑68 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCg                                                          >  1:1304559/1‑67 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGt                                                         >  1:3109256/1‑68 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGt                                                         >  1:257588/1‑68 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGt                                                         >  1:1391267/1‑68 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGt                                                         >  1:190338/1‑68 (MQ=255)
            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTa                                                        >  1:2424084/1‑69 (MQ=255)
                         ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTAtt                                                      <  1:3362527/58‑1 (MQ=255)
                              aCCATTTTCTTCACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGtt                                      >  1:840556/1‑69 (MQ=255)
                              aCCATTTTCTTCACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGtt                                      >  1:1743462/1‑69 (MQ=255)
                                   tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCg                                 <  1:653378/69‑1 (MQ=255)
                                   tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCg                                 <  1:338206/69‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCGTTGAGGCCGTAGTTTAc                                   <  1:3272108/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:871705/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2758550/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:799789/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2838465/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2873702/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:688253/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:301720/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:3052200/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:3093049/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:592336/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:3194919/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:3380881/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2753399/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1125042/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1350291/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:14841/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1619740/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1656996/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1744018/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1957408/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1960460/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2702946/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2616220/65‑1 (MQ=255)
                                     tcttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2577680/65‑1 (MQ=255)
                                      cttcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2971062/64‑1 (MQ=255)
                                        tcACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:319876/62‑1 (MQ=255)
                                         cacaCGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGATCAGt                           >  1:2756158/1‑69 (MQ=255)
                                         cacaCGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGt                           >  1:255302/1‑69 (MQ=255)
                                         cacaCGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGt                           >  1:1919032/1‑69 (MQ=255)
                                                aTCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:1401356/54‑1 (MQ=255)
                                                   accAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTAc                                   <  1:2936259/50‑1 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:3113155/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:902332/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:1196652/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:1229318/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:3010784/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:782268/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:68579/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:2784022/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:529345/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:473729/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:389537/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:380015/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:370606/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:2639000/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:3423026/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:1695219/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:3331415/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:1923788/1‑69 (MQ=255)
                                                    ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                >  1:230099/1‑69 (MQ=255)
                                                      aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACg              <  1:1523618/69‑1 (MQ=255)
                                                      aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACg              <  1:2124176/69‑1 (MQ=255)
                                                         cGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGt           >  1:2364673/1‑69 (MQ=255)
                                                         cGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGt           >  1:3208619/1‑69 (MQ=255)
                                                                  aaTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa  <  1:2483269/69‑1 (MQ=255)
                                                                  aaTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa  <  1:2487157/69‑1 (MQ=255)
                                                                  |                                                                    
CGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATA  >  minE/1971823‑1971957

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: