Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,471 A→G 14.6% T75T (ACT→ACC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4710AG14.6% 232.5 / 28.5 96T75T (ACT→ACCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (44/38);  new base G (8/6);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTT  >  minE/2919401‑2919538
                                                                      |                                                                   
gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                                                     <  1:3385590/71‑1 (MQ=255)
gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                                                     <  1:1545066/71‑1 (MQ=255)
  ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                    >  1:2606645/1‑70 (MQ=255)
  ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                    >  1:3253438/1‑70 (MQ=255)
  ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                    >  1:1435759/1‑70 (MQ=255)
  ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                                                    >  1:3744153/1‑70 (MQ=255)
   cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTa                                                                  <  1:1338284/71‑1 (MQ=255)
   cAAGCACGCGGGCAGTGTCTTTAATAGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                                                                   >  1:3808790/1‑70 (MQ=255)
     aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                                >  1:2019871/1‑71 (MQ=255)
     aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                                <  1:800302/71‑1 (MQ=255)
       cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                              <  1:198903/71‑1 (MQ=255)
       cACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcg                                                               >  1:1024561/1‑70 (MQ=255)
         cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                            <  1:1212140/71‑1 (MQ=255)
          gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                           >  1:3674753/1‑71 (MQ=255)
          gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                           >  1:2676850/1‑71 (MQ=255)
               cAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCAccc                                                      >  1:1171111/1‑71 (MQ=255)
                 gtgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGccc                                                      <  1:1796656/69‑1 (MQ=255)
                  tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                                 >  1:2262576/1‑57 (MQ=255)
                  tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                   >  1:3772954/1‑71 (MQ=255)
                  tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                   >  1:305225/1‑71 (MQ=255)
                  tgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAAc                                                                >  1:1536508/1‑58 (MQ=255)
                   gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg                                                   <  1:1451349/70‑1 (MQ=255)
                   gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGt                                                  >  1:2564401/1‑71 (MQ=255)
                   gtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAAcgc                                                              <  1:1223077/59‑1 (MQ=255)
                    tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                            >  1:3227485/1‑60 (MQ=255)
                    tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                            >  1:140120/1‑60 (MQ=255)
                    tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                            >  1:135462/1‑60 (MQ=255)
                      tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                                                                   <  1:3029960/51‑1 (MQ=255)
                      tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                              <  1:2679631/56‑1 (MQ=255)
                      tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                                           >  1:456162/1‑59 (MQ=255)
                      tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                >  1:3480390/1‑70 (MQ=255)
                      tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa                                                >  1:910340/1‑70 (MQ=255)
                        tAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCaaat                                             >  1:3334071/1‑71 (MQ=255)
                        tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCatat                                             >  1:1032572/1‑71 (MQ=255)
                        tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAAGCGAGCACCCTGGTCatat                                             >  1:3624586/1‑71 (MQ=255)
                         aaTTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCt                                                     >  1:2810014/1‑62 (MQ=255)
                          aTCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAc                                                                <  1:3024510/50‑1 (MQ=255)
                           tCGACTCTTTATGACCAATCTGGCTGACGCTTGGGCCGAGATAAGTAACGCGAGCACCCt                                                     >  1:1648351/1‑60 (MQ=255)
                              aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                                 <  1:1961293/45‑1 (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)
                                     aTGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                                                                 >  1:2501674/1‑38 (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:773765/70‑1 (MQ=255)
                                      tGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGa                                <  1:3506041/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)
                                        aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                              <  1:1027469/38‑1 (MQ=255)
                                        aCCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                                            >  1:965550/1‑40 (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)
                                          cAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                                                              >  1:657330/1‑36 (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)
                                              cTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCa                                                         >  1:2226734/1‑37 (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)
                                                                      |                                                                   
GACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTT  >  minE/2919401‑2919538

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: