Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,350,843 T→C 18.6% intergenic (‑22/+70) dcd ← / ← udk 2'‑deoxycytidine 5'‑triphosphate deaminase/uridine/cytidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,350,8430TC18.6% 182.5 / 35.3 86intergenic (‑22/+70)dcd/udk2'‑deoxycytidine 5'‑triphosphate deaminase/uridine/cytidine kinase
Reads supporting (aligned to +/- strand):  ref base T (48/22);  new base C (7/9);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 8.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

AACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAA  >  minE/1350773‑1350910
                                                                      |                                                                   
aaCGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCAt                                                                     >  1:1993530/1‑71 (MQ=255)
aaCGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCAt                                                                     >  1:2876854/1‑71 (MQ=255)
 aCGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCAt                                                                     >  1:2313582/1‑70 (MQ=255)
  cGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCAtta                                                                   >  1:3631240/1‑69 (MQ=255)
    gCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACt                                                                 <  1:768425/71‑1 (MQ=255)
     ccTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTc                                                                >  1:34156/1‑71 (MQ=255)
      cTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTcc                                                               >  1:2433422/1‑71 (MQ=255)
        tcatcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCt                                                             >  1:2862699/1‑71 (MQ=255)
        tcatcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCt                                                             >  1:1602167/1‑71 (MQ=255)
        tcataaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCt                                                             >  1:2014897/1‑71 (MQ=255)
         catcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCt                                                             >  1:2772304/1‑70 (MQ=255)
          atcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCtt                                                            >  1:1764420/1‑70 (MQ=255)
          atcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCtt                                                            <  1:1577843/70‑1 (MQ=255)
          atcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTa                                                           <  1:1208749/71‑1 (MQ=255)
          atcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTa                                                           <  1:1065560/71‑1 (MQ=255)
           tcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          <  1:3698015/71‑1 (MQ=255)
           tcaAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          <  1:3946557/71‑1 (MQ=255)
            caAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          >  1:357220/1‑70 (MQ=255)
                 cAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCaa                                                    <  1:2548687/71‑1 (MQ=255)
                    gCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTa                                                 <  1:2554530/71‑1 (MQ=255)
                      ttCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGcc                                                      <  1:1099389/64‑1 (MQ=255)
                       tCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGt                                              >  1:1326534/1‑71 (MQ=255)
                        cAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGt                                              >  1:1821159/1‑70 (MQ=255)
                        cAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGt                                              >  1:1328152/1‑70 (MQ=255)
                        cAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGt                                              >  1:838185/1‑70 (MQ=255)
                         aaTATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTc                                             <  1:176054/70‑1 (MQ=255)
                         aaTATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCa                                            <  1:1469027/71‑1 (MQ=255)
                         aaTATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACGACTCCCTTAACGCCAATTACGTCa                                            <  1:3767254/71‑1 (MQ=255)
                         aaTATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTc                                             <  1:16986/70‑1 (MQ=255)
                          atatCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCaa                                           <  1:629737/71‑1 (MQ=255)
                           tatCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCaa                                           <  1:389817/70‑1 (MQ=255)
                             tctcGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTccc                                                              >  1:1831198/1‑49 (MQ=255)
                              ctcGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAggg                                        >  1:1131014/1‑70 (MQ=255)
                               tcGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          >  1:559899/1‑51 (MQ=255)
                                 ggTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCAt                                     >  1:1804456/1‑70 (MQ=255)
                                 ggTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATa                                    >  1:198943/1‑71 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAgg                                         <  1:1845187/65‑1 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAgg                                         <  1:1350703/65‑1 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATAc                                   <  1:1187550/71‑1 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACg                                                        >  1:3683967/1‑50 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGc                                                       >  1:2141933/1‑51 (MQ=255)
                                  gTCACACAGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGc                                                       >  1:2422955/1‑51 (MQ=255)
                                   tCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCa                                            <  1:3924138/61‑1 (MQ=255)
                                    cacacaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGc                                       >  1:3862087/1‑65 (MQ=255)
                                     acacaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAAtt                                                  <  1:3904639/53‑1 (MQ=255)
                                     acacaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTa                                 >  1:449216/1‑70 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          >  1:2818577/1‑43 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTaa                                                          >  1:3599084/1‑43 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGt                              >  1:3213197/1‑71 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGt                              >  1:1236359/1‑71 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGt                              >  1:3406327/1‑71 (MQ=255)
                                       acaGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGt                              >  1:2040442/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTccc                                                              >  1:729150/1‑37 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCtt                                                            >  1:424587/1‑39 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            >  1:2056818/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            >  1:306375/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            >  1:311026/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            >  1:3790755/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            >  1:3274417/1‑71 (MQ=255)
                                         aGACGCATGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTa                                                 >  1:2377825/1‑50 (MQ=255)
                                          gACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTAt                           <  1:2797224/71‑1 (MQ=255)
                                           aCGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTAt                           >  1:2115895/1‑70 (MQ=255)
                                           aCGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTAt                           >  1:2946111/1‑70 (MQ=255)
                                            cGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTa                            <  1:2407175/68‑1 (MQ=255)
                                             gCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAg                                          >  1:2868526/1‑53 (MQ=255)
                                               aTGGCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACt                       <  1:2461623/70‑1 (MQ=255)
                                                tGGTATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCAc                        <  1:3431698/68‑1 (MQ=255)
                                                 ggCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGaa                    >  1:1512103/1‑71 (MQ=255)
                                                  gCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCAc                        >  1:2614547/1‑66 (MQ=255)
                                                  gCATTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTg                      >  1:3849321/1‑68 (MQ=255)
                                                    aTTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAg                                          <  1:1873494/46‑1 (MQ=255)
                                                    aTTTTCTCCTTTCGCGCACCACTCCCTTAACGCCAATTACGTCAAg                                          <  1:1457832/46‑1 (MQ=255)
                                                      tttCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTa               >  1:685638/1‑71 (MQ=255)
                                                      tttCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTa               >  1:884207/1‑71 (MQ=255)
                                                      tttCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTa               >  1:1312930/1‑71 (MQ=255)
                                                         ctcCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACgg                 <  1:3822375/66‑1 (MQ=255)
                                                          tcCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTacac            >  1:150559/1‑70 (MQ=255)
                                                             tttCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTaca             <  1:3362000/66‑1 (MQ=255)
                                                             tttCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTaca             <  1:2031850/66‑1 (MQ=255)
                                                               tCGCGCATCACTCCCTTAACGCCAATTACGTCAAg                                          >  1:3528252/1‑35 (MQ=255)
                                                                cgcgCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTa      <  1:3207891/70‑1 (MQ=255)
                                                                cgcgCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTAt     >  1:951921/1‑71 (MQ=255)
                                                                cgcgCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTAt     >  1:2681486/1‑71 (MQ=255)
                                                                cgcgCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTAt     >  1:2592553/1‑71 (MQ=255)
                                                                cgcgCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTAt     >  1:1696218/1‑71 (MQ=255)
                                                                   gcATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTAt     >  1:400007/1‑68 (MQ=255)
                                                                   gcATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATc    >  1:3138585/1‑69 (MQ=255)
                                                                   gcATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATc    >  1:2043158/1‑69 (MQ=255)
                                                                   gcACCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCaa  <  1:414716/71‑1 (MQ=255)
                                                                    cATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGCACGGt                >  1:1403502/1‑56 (MQ=255)
                                                                      |                                                                   
AACGGCCTTCATCAAGCCAGGCTTCAATATCTCGGTCACACAGACGCATGGCATTTTCTCCTTTCGCGCATCACTCCCTTAACGCCAATTACGTCAAGGGCATACTAGGTTATCACTGAACGGTACACAATTTATCAA  >  minE/1350773‑1350910

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: