Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,493 T→A 22.0% Y68F (TAT→TTT)  argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4930TA22.0% 187.8 / 51.1 91Y68F (TAT→TTT) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (26/45);  new base A (11/9);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAC  >  minE/2919425‑2919562
                                                                    |                                                                     
tAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCaaat                                                                     >  1:3334071/1‑71 (MQ=255)
tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCatat                                                                     >  1:1032572/1‑71 (MQ=255)
tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAAGCGAGCACCCTGGTCatat                                                                     >  1:3624586/1‑71 (MQ=255)
      aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCa                                                               <  1:1679702/71‑1 (MQ=255)
       ctctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCaa                                                              >  1:3308221/1‑71 (MQ=255)
        tctTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGc                                                                >  1:1662501/1‑68 (MQ=255)
        tctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCaa                                                              >  1:2704415/1‑70 (MQ=255)
        tctTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAAc                                                             >  1:788008/1‑71 (MQ=255)
          tttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAAc                                                             <  1:962910/69‑1 (MQ=255)
          tttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCATGGTCATATGCGGCAAc                                                             <  1:1954740/69‑1 (MQ=255)
           ttATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGc                                                                <  1:2668362/65‑1 (MQ=255)
            tATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTc                                                          <  1:339524/70‑1 (MQ=255)
            tATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTc                                                          >  1:575732/1‑70 (MQ=255)
             aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCg                                                         <  1:417330/70‑1 (MQ=255)
             aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCg                                                         >  1:3808433/1‑70 (MQ=255)
              tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                                       <  1:1746299/71‑1 (MQ=255)
              tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                                                        <  1:3506041/70‑1 (MQ=255)
              tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                                                        <  1:773765/70‑1 (MQ=255)
               gACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGc                                                                   <  1:503415/58‑1 (MQ=255)
               gACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGAACCCTTGACATATGc                                                                   <  1:1341441/58‑1 (MQ=255)
                cccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                                      <  1:3064838/69‑1 (MQ=255)
                cccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAag                                                     <  1:3312819/70‑1 (MQ=255)
                aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaaa                                                      <  1:244577/70‑1 (MQ=255)
                aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAag                                                     <  1:1603931/71‑1 (MQ=255)
                aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAag                                                     <  1:2414506/71‑1 (MQ=255)
                 ccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                                      >  1:1550738/1‑69 (MQ=255)
                 ccAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGaaa                                                      >  1:3400372/1‑69 (MQ=255)
                  cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCCCCCTGGTCATATGCGGCAACTTCGAAagag                                                   >  1:2090553/1‑71 (MQ=255)
                  cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAagag                                                   >  1:2074159/1‑71 (MQ=255)
                    aTCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCa                                                 <  1:3100397/71‑1 (MQ=255)
                     tCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCAt                                                >  1:1266996/1‑71 (MQ=255)
                       tGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGc                                                                   >  1:2709100/1‑50 (MQ=255)
                       tGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCg                                              >  1:1136939/1‑71 (MQ=255)
                         gCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATc                                               <  1:988330/68‑1 (MQ=255)
                          cTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGt                                           <  1:1150435/71‑1 (MQ=255)
                             ccGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                         >  1:1021955/1‑70 (MQ=255)
                             ccGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAc                                         >  1:1352634/1‑70 (MQ=255)
                              cGCTTGGGCCGAGATAAGTAACGCGAGCACGCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                         <  1:1755589/69‑1 (MQ=255)
                              cGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                       <  1:47219/71‑1 (MQ=255)
                              cGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAaga                                                    >  1:1256399/1‑58 (MQ=255)
                               gCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgag                                      <  1:3292804/71‑1 (MQ=255)
                                cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATg                                                                    <  1:2390389/40‑1 (MQ=255)
                                cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATc                                               <  1:3930069/61‑1 (MQ=255)
                                cTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagt                                     >  1:175560/1‑71 (MQ=255)
                                  tGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                    >  1:764412/1‑70 (MQ=255)
                                   gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcg                                   <  1:2916072/70‑1 (MQ=255)
                                   gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                    <  1:733628/69‑1 (MQ=255)
                                   gggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                    <  1:3227791/69‑1 (MQ=255)
                                    ggCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcgag                                 <  1:509753/71‑1 (MQ=255)
                                      ccGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACtt                                                           >  1:3289058/1‑43 (MQ=255)
                                       cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                                       <  1:2694615/46‑1 (MQ=255)
                                       cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                                       <  1:519303/46‑1 (MQ=255)
                                       cGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCt                              >  1:676457/1‑71 (MQ=255)
                                        gagaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGt                                           >  1:3365593/1‑57 (MQ=255)
                                          gaTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCtttt                           >  1:2066507/1‑71 (MQ=255)
                                           aTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACt                                                            <  1:1050373/37‑1 (MQ=255)
                                           aTAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCttttt                          <  1:1716666/71‑1 (MQ=255)
                                             aaGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCa                                                 <  1:2993297/46‑1 (MQ=255)
                                               gTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGaa                      <  1:2915227/71‑1 (MQ=255)
                                                tAACGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                                       <  1:2984119/37‑1 (MQ=255)
                                                 aaCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAcgagtcg                                   <  1:1874776/56‑1 (MQ=255)
                                                  aCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGAt                   <  1:3790563/71‑1 (MQ=255)
                                                  aCGCGAGCACCCTGGTCATATGCGGCAACTTCGaa                                                       <  1:3311641/35‑1 (MQ=255)
                                                  aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTc                         <  1:635119/65‑1 (MQ=255)
                                                  aCGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt                   <  1:902519/71‑1 (MQ=255)
                                                  aCGCGAGCACCCTGGACATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt                   <  1:3738397/71‑1 (MQ=255)
                                                   cgcgcgcgCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCGTTTTCGAAGAt                   >  1:1109394/1‑70 (MQ=255)
                                                   cgcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGt                                           <  1:278775/46‑1 (MQ=255)
                                                    gcgcgcgcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGAt                   <  1:590861/69‑1 (MQ=255)
                                                    gcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGAt                   >  1:1052411/1‑69 (MQ=255)
                                                    gcgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGa                 <  1:360040/71‑1 (MQ=255)
                                                     cgAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAg                <  1:2008898/71‑1 (MQ=255)
                                                      gAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAg                     >  1:1141367/1‑65 (MQ=255)
                                                        gcgcccTGGTCAAATGCGGGAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAgcgc             >  1:45921/4‑71 (MQ=255)
                                                        gcgcccTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAgcgc             >  1:3311355/4‑71 (MQ=255)
                                                         cACCCTGGTCATATGCTGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGcgcg            <  1:2405417/71‑1 (MQ=255)
                                                         cACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAc                                         <  1:764474/42‑1 (MQ=255)
                                                         cACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAgc               <  1:2483171/68‑1 (MQ=255)
                                                          gcccTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGa           <  1:2332387/70‑1 (MQ=255)
                                                          gcccTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGa           <  1:3713242/70‑1 (MQ=255)
                                                           cccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGTGTCGAGTCTTTTTCGAAGATGAGCGCGAt          >  1:3537270/1‑71 (MQ=255)
                                                           cccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGcgcg            >  1:3466483/1‑69 (MQ=255)
                                                           cccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGAt          >  1:499683/1‑71 (MQ=255)
                                                            ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                       >  1:616707/1‑41 (MQ=255)
                                                            ccTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCttt                            <  1:274598/52‑1 (MQ=255)
                                                              tGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcg                                        <  1:2148317/38‑1 (MQ=255)
                                                              tGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTa                                          >  1:1809294/1‑36 (MQ=38)
                                                              tGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGaa                      >  1:3894423/1‑56 (MQ=255)
                                                                gTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgag                                      >  1:91274/1‑38 (MQ=255)
                                                                gTCATATGCGGCAACTTCGAAAGAGCATCGGGTAcga                                       >  1:874462/1‑37 (MQ=255)
                                                                gTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGttt      <  1:3354898/70‑1 (MQ=255)
                                                                gTCAAATGCGGCAACTTCGAAAGAGCAACGGGTAcgagtcgagt                                >  1:2425440/1‑44 (MQ=255)
                                                                 tCATATGCGGCAACTTCGAAAGAGCATCGGGTAcgagtc                                    <  1:1912777/39‑1 (MQ=255)
                                                                 tCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGAt          >  1:2341495/1‑65 (MQ=255)
                                                                   atatGCGGCAACTTCGAAAGAGCATCGGGTAcgagtcga                                  >  1:3707273/1‑39 (MQ=255)
                                                                   atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGt        <  1:2381121/65‑1 (MQ=255)
                                                                   atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc  >  1:1518154/1‑71 (MQ=255)
                                                                   atatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc  >  1:2691091/1‑71 (MQ=255)
                                                                    tatGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGa                 <  1:2539659/55‑1 (MQ=255)
                                                                    |                                                                     
TAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAC  >  minE/2919425‑2919562

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: