Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,516 T→A 14.6% R60R (CGA→CGT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5160TA14.6% 249.7 / 32.2 96R60R (CGA→CGTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (28/54);  new base A (7/7);  total (35/61)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTT  >  minE/2919450‑2919584
                                                                  |                                                                    
gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCAt                                                                      <  1:318720/67‑1 (MQ=255)
gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATc                                                                     <  1:356214/68‑1 (MQ=255)
gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCggg                                                                  <  1:3869213/71‑1 (MQ=255)
gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCggg                                                                  <  1:1261065/71‑1 (MQ=255)
  tGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTa                                                                >  1:4116215/1‑71 (MQ=255)
   gCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                                               >  1:581768/1‑71 (MQ=255)
    ccGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                                               <  1:1212922/70‑1 (MQ=255)
     cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAATAGCATCGGGTAcga                                                             <  1:1028387/71‑1 (MQ=255)
     cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                                             <  1:3570112/71‑1 (MQ=255)
     cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                                             <  1:504669/71‑1 (MQ=255)
     cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                                             <  1:683439/71‑1 (MQ=255)
         tGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                                          >  1:724116/1‑70 (MQ=255)
          gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                                          <  1:2970995/69‑1 (MQ=255)
            gCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCggg                                                                  <  1:1637446/59‑1 (MQ=255)
              cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                                        >  1:3624535/1‑67 (MQ=255)
               gagaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtt                                                   >  1:1292363/1‑71 (MQ=255)
               gagaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtt                                                   >  1:3634090/1‑71 (MQ=255)
                agaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCttt                                                  <  1:728890/71‑1 (MQ=255)
                agaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCttt                                                  <  1:3711207/71‑1 (MQ=255)
                 gaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtttt                                                 >  1:2520862/1‑71 (MQ=255)
                 aaTAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCtttt                                                 >  1:3865966/2‑71 (MQ=255)
                  aTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                                          <  1:964656/61‑1 (MQ=255)
                   tAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                                             <  1:1678567/57‑1 (MQ=255)
                    aaGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                                        >  1:2968530/1‑61 (MQ=255)
                      gTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:676591/71‑1 (MQ=255)
                      gTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:2659332/71‑1 (MQ=255)
                      gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                                        >  1:1808338/1‑59 (MQ=255)
                      gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:3770557/71‑1 (MQ=255)
                      gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:2809018/71‑1 (MQ=255)
                      gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:2317538/71‑1 (MQ=255)
                       tAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGGCTTTTTCGAAg                                           <  1:1855672/71‑1 (MQ=255)
                        aaCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCAt                                                                      >  1:3342155/1‑43 (MQ=255)
                        aaCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCg                                              >  1:2955093/1‑67 (MQ=255)
                        aaCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAg                                           <  1:100459/70‑1 (MQ=255)
                         aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:403243/68‑1 (MQ=255)
                         aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:3656974/68‑1 (MQ=255)
                          cgcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATg                                        <  1:1830995/71‑1 (MQ=255)
                          cgcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATg                                        <  1:2924942/71‑1 (MQ=255)
                            cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt                                         >  1:2283473/1‑68 (MQ=255)
                            cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGa                                       >  1:2698401/1‑70 (MQ=255)
                            cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAg                                      >  1:2666823/1‑71 (MQ=255)
                            cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTAGAAGATg                                        >  1:688405/1‑69 (MQ=255)
                               gCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                                        >  1:555993/1‑50 (MQ=255)
                                cgcccTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGcgcg                                  <  1:1884764/69‑1 (MQ=255)
                                   ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                                               <  1:1242539/39‑1 (MQ=255)
                                   ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                                               <  1:30519/39‑1 (MQ=255)
                                   ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtttt                                                 <  1:2121779/53‑1 (MQ=255)
                                   ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATg                                        <  1:2486895/62‑1 (MQ=255)
                                   ccTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATg                               <  1:1777323/71‑1 (MQ=255)
                                    cTGGTCATATGCGGCAACTTCGAAAGAGCATCggg                                                                  <  1:3336708/35‑1 (MQ=255)
                                    cTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGa                                          >  1:2282761/1‑59 (MQ=255)
                                    cTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATg                                        <  1:329397/61‑1 (MQ=255)
                                     tGGTCATATGCGGCAACTTCGAAAGAGCATCGCGTACGAGTCGAGTCt                                                    >  1:2178340/1‑48 (MQ=255)
                                     tGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGtt                             >  1:1681307/1‑71 (MQ=255)
                                      ggTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAcgagt                                                           >  1:3059741/1‑40 (MQ=255)
                                      ggTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGttt                            >  1:3268068/1‑71 (MQ=255)
                                       gTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGtt                             <  1:2059077/69‑1 (MQ=255)
                                        tCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGttttt                          <  1:2276910/71‑1 (MQ=255)
                                          atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            <  1:1994963/51‑1 (MQ=255)
                                          atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGa                                 >  1:143603/1‑62 (MQ=255)
                                          atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                        <  1:3429740/71‑1 (MQ=255)
                                          atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                        >  1:2623331/1‑71 (MQ=255)
                                           tatGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                                          <  1:956578/36‑1 (MQ=255)
                                           tatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGttttt                          <  1:3025664/68‑1 (MQ=255)
                                             tGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                                        >  1:3419397/1‑36 (MQ=255)
                                             tGCGGCAACTTCGAAAGAGCATCGGGTACgagtcgagtc                                                     >  1:3520141/1‑39 (MQ=255)
                                             tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            >  1:1185471/1‑48 (MQ=255)
                                              gCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAccgg                     >  1:3511448/1‑68 (MQ=255)
                                               cGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGtt                             <  1:4107318/61‑1 (MQ=255)
                                                ggCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGAAGTTTTTACCAGt                    <  1:3453768/69‑1 (MQ=255)
                                                 gCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGtttt                           <  1:892166/61‑1 (MQ=255)
                                                   aaCTTCGAAAGAGCATCGGGTACGAGTCGAGTCttt                                                  <  1:3286013/36‑1 (MQ=255)
                                                   aaCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTc                                               >  1:1690998/1‑39 (MQ=255)
                                                   aaCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                                            >  1:815364/1‑42 (MQ=255)
                                                     cTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg             <  1:457845/71‑1 (MQ=255)
                                                     cTTCGAAAGAGCAACGGGTACGATTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttc             <  1:2919292/71‑2 (MQ=255)
                                                      ttCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                 <  1:2729500/66‑1 (MQ=255)
                                                      ttCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg             >  1:3770343/1‑70 (MQ=255)
                                                       tCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                       >  1:1375800/1‑59 (MQ=255)
                                                        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg           >  1:1068263/1‑67 (MQ=255)
                                                        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg           >  1:392637/1‑67 (MQ=255)
                                                        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttct            <  1:1149638/69‑3 (MQ=255)
                                                         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAg                                           <  1:2062205/37‑1 (MQ=255)
                                                         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAg                                           <  1:130122/37‑1 (MQ=255)
                                                         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGctt         <  1:952633/71‑1 (MQ=255)
                                                         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct          <  1:1519431/70‑1 (MQ=255)
                                                           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct       <  1:467653/71‑1 (MQ=255)
                                                           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct       <  1:3823534/71‑1 (MQ=255)
                                                           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttc        <  1:3578458/70‑1 (MQ=255)
                                                           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGACTTTGGcttct       >  1:3374360/1‑71 (MQ=255)
                                                            agagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGt                              <  1:1031158/47‑1 (MQ=255)
                                                            agagCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgc                <  1:905995/61‑2 (MQ=255)
                                                             gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTa                         >  1:131459/1‑51 (MQ=255)
                                                             gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                        <  1:330008/52‑1 (MQ=255)
                                                             gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt              <  1:100857/62‑1 (MQ=255)
                                                               gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt              <  1:2674380/60‑1 (MQ=255)
                                                               gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTCTACCAGTGAGTTTGGcttcttctt   <  1:1939964/71‑1 (MQ=255)
                                                                cATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCatt  <  1:3717555/71‑1 (MQ=255)
                                                                  tCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTa                         <  1:2780335/46‑1 (MQ=255)
                                                                  |                                                                    
GCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTT  >  minE/2919450‑2919584

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: