Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,987,051 T→C 30.9% P224P (CCA→CCG ansB ← periplasmic L‑asparaginase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,987,0510TC30.9% 73.1 / 40.7 84P224P (CCA→CCGansBperiplasmic L‑asparaginase II
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (11/15);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTAACC  >  minE/1986998‑1987119
                                                     |                                                                    
tAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGt                                                        >  1:1530912/1‑68 (MQ=255)
tAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTa                                                       >  1:1006591/1‑69 (MQ=255)
    tataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCt                                                   <  1:2839745/69‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:1879865/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2096322/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2142474/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:1744491/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:3180790/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:362454/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:1407043/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:1183329/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2581750/62‑1 (MQ=255)
     ataAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:969495/62‑1 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:2488577/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:276603/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:742070/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:3325630/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:692814/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:1672317/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:3377579/1‑69 (MQ=255)
          cAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTg                                             >  1:1694766/1‑69 (MQ=255)
           aaTGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGc                                            <  1:3052328/69‑1 (MQ=255)
              gCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCgg                                         >  1:735802/1‑69 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCTCTgg                                                         <  1:1645576/52‑1 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2874759/52‑1 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2185392/52‑1 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2123450/52‑1 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:2432778/52‑1 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         >  1:2727391/1‑52 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         >  1:1112687/1‑52 (MQ=255)
               ccGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTgg                                                         <  1:1611767/52‑1 (MQ=255)
                      tCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCggg                                        >  1:1828350/1‑62 (MQ=255)
                      tCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCggg                                        >  1:3182031/1‑62 (MQ=255)
                      tCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCggg                                        >  1:2660590/1‑62 (MQ=255)
                      tCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCggg                                        >  1:230536/1‑62 (MQ=255)
                      tCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCgg                                         >  1:2146014/1‑61 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:1567993/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:2585911/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:1255278/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:105158/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:2030374/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:387053/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:2692422/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:1769901/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:3321722/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTa                             >  1:3186810/1‑68 (MQ=255)
                           aGTTCATTCAGCTTAGAGACATCGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGTACGTTGGTAg                            >  1:2073426/1‑68 (MQ=255)
                                   caGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCtt                    >  1:509154/1‑69 (MQ=255)
                                   caGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCtt                    >  1:2894493/1‑69 (MQ=255)
                                   caGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCtt                    >  1:2301292/1‑69 (MQ=255)
                                              cATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGaa         >  1:926271/1‑69 (MQ=255)
                                              cATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGaa         >  1:2570172/1‑69 (MQ=255)
                                              cATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGaa         >  1:999020/1‑69 (MQ=255)
                                              cATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGaa         >  1:2675514/1‑69 (MQ=255)
                                              cATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGaa         >  1:2172885/1‑69 (MQ=255)
                                                 cGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGt      >  1:640065/1‑69 (MQ=255)
                                                 cGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGt      >  1:2329585/1‑69 (MQ=255)
                                                 cGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGt      >  1:2615751/1‑69 (MQ=255)
                                                 cGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGt      >  1:257021/1‑69 (MQ=255)
                                                 cGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGt      >  1:1313076/1‑69 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1685137/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1043050/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:114768/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:868286/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1191071/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1425128/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1446185/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:150972/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:2652967/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:208358/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:2171164/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:2235658/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:2595921/69‑1 (MQ=255)
                                                 cGAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACCTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:1043819/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:1780537/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:3292428/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:1455442/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:724096/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:1296992/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:3101124/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:809599/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:2927329/69‑1 (MQ=255)
                                                  gAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTa     <  1:2637750/69‑1 (MQ=255)
                                                  gAACGGCGTGTCGCTGGTGTGCTTACGTGCCGGGGTACGTTGGTAGTCAATCTTACCGTTGTGAATGt      <  1:317566/68‑1 (MQ=255)
                                                     tGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTAAcc  <  1:1877028/69‑1 (MQ=255)
                                                     |                                                                    
TAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCAGCTTAGAGACATCGAATGGCGTGTCGCTGGTATGCTTACGTGCCGGGGTACGCTGGTAGTCAATCTTACCGTTGTGAATGTAACC  >  minE/1986998‑1987119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: