Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,971,856 G→A 20.7% N350N (AAC→AAT speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,971,8560GA20.7% 104.2 / 17.1 72N350N (AAC→AATspeAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (38/19);  new base A (2/13);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 2.86e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

AGCACGGTGTGATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTT  >  minE/1971791‑1971921
                                                                 |                                                                 
aGCACGGTGTGATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGttt                                                                >  1:2876539/1‑69 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTctca                                                        <  1:1465710/69‑1 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:102579/69‑1 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:533952/69‑1 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:2165122/69‑1 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:1431266/69‑1 (MQ=255)
        gtgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:1426949/69‑1 (MQ=255)
         tgATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTcaca                                                        <  1:2892257/68‑1 (MQ=255)
           aTGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGc                                                     >  1:1343666/1‑69 (MQ=255)
           aTGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGc                                                     >  1:731/1‑69 (MQ=255)
           aTGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGc                                                     >  1:489539/1‑69 (MQ=255)
              cgcAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACGCGCATc                                                  >  1:1004015/1‑69 (MQ=255)
              cgcAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATc                                                  >  1:2072249/1‑69 (MQ=255)
              cgcAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATc                                                  >  1:2006646/1‑69 (MQ=255)
                  gTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCa                                              <  1:1427251/69‑1 (MQ=255)
                   tCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCa                                              <  1:2918839/68‑1 (MQ=255)
                     acagCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:2234186/4‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACCCGCATCGCCAAt                                            >  1:2003347/1‑68 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:741371/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:377413/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:2641252/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:3281380/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:96170/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:2829949/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:3239238/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:1515087/1‑69 (MQ=255)
                     aCCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:3224101/1‑69 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATc                                           >  1:1873272/1‑67 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2819737/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2417785/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2627913/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2500588/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2491250/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:798861/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:2165584/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:1452869/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:1760512/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCg                                          >  1:1760810/1‑68 (MQ=255)
                       cGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTGTTCTTCACACGCATCGCCAATc                                           >  1:2533810/1‑67 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGTAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:1122073/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:738248/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:1051191/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:3266813/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:3001552/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:773605/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:2399743/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGAt                                   <  1:2403925/69‑1 (MQ=255)
                             aCCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCAAACGCATCGCCAATCGCCCAGAt                                   <  1:3364/69‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2828910/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2128698/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:1681679/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2947501/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2988097/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:275733/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2720676/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:3168578/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:1075123/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:2542734/68‑1 (MQ=255)
                                agATTCGGTGATTACCGTCGGATGTGGCAGACCATTTTCTTCACACGCATCACCAATCGCCCAgatgat                                <  1:361707/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:190338/1‑68 (MQ=255)
                                            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGt                     >  1:1391267/1‑68 (MQ=255)
                                            tACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTa                    >  1:2424084/1‑69 (MQ=255)
                                                         ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTAtt                  <  1:3362527/58‑1 (MQ=255)
                                                              aCCATTTTCTTCACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGtt  >  1:1743462/1‑69 (MQ=255)
                                                              aCCATTTTCTTCACACGCATCACCAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGtt  >  1:840556/1‑69 (MQ=255)
                                                                 |                                                                 
AGCACGGTGTGATGCGCAGTCACCGCACGACCCGATTCGGTGATTACCGTCGGATGCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTT  >  minE/1971791‑1971921

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: