Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,441 A→C 12.8% G85G (GGT→GGG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4410AC12.8% 197.8 / 27.2 94G85G (GGT→GGGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base C (7/5);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAG  >  minE/2919373‑2919511
                                                                    |                                                                      
ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                                                        >  1:2516378/1‑69 (MQ=255)
ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                                                       >  1:287660/1‑70 (MQ=255)
ccGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                                                       >  1:3877993/1‑70 (MQ=255)
 cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAcc                                                                      >  1:3140073/1‑70 (MQ=255)
   cgATACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAAt                                                                   >  1:3343308/1‑71 (MQ=255)
   cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                                                   <  1:1072865/71‑1 (MQ=255)
   cgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAAt                                                                   <  1:1501430/71‑1 (MQ=255)
     aTACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                                                                 <  1:3140774/71‑1 (MQ=255)
      tACTTAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCt                                                                 <  1:3182861/70‑1 (MQ=255)
      tACTGAATGCCGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTg                                                                >  1:1151001/1‑71 (MQ=255)
           aaaaCCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                                            <  1:2421608/67‑1 (MQ=255)
           aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg                                                               >  1:1196747/1‑67 (MQ=255)
           aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                                            <  1:1588042/70‑1 (MQ=255)
           aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                                            <  1:3400262/70‑1 (MQ=255)
           aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                                            <  1:3576235/70‑1 (MQ=255)
             tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                                                         <  1:1014431/71‑1 (MQ=255)
             tACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCg                                                         >  1:2875305/1‑71 (MQ=255)
               ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                                       <  1:530444/71‑1 (MQ=255)
                cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                                                      >  1:390611/1‑71 (MQ=255)
                cGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                                                      <  1:3216134/71‑1 (MQ=255)
                  tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTgg                                                    <  1:3916074/71‑1 (MQ=255)
                  tCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCgg                                                    >  1:3047201/1‑71 (MQ=255)
                     tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                                                 <  1:3711769/71‑1 (MQ=255)
                       cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGcccta                                               <  1:1101366/71‑3 (MQ=255)
                        aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                                                 >  1:3325441/1‑68 (MQ=255)
                        aTGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                                                 >  1:2753508/1‑68 (MQ=255)
                          gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                             >  1:1671216/1‑70 (MQ=255)
                          gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                             >  1:2474233/1‑70 (MQ=255)
                          gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgaga                                             >  1:286668/1‑70 (MQ=255)
                           cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                                            >  1:628717/1‑70 (MQ=255)
                            gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                          <  1:1545066/71‑1 (MQ=255)
                            gACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                                          <  1:3385590/71‑1 (MQ=255)
                              ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                                                          <  1:1916688/53‑1 (MQ=255)
                              ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                         >  1:3253438/1‑70 (MQ=255)
                              ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                         >  1:1435759/1‑70 (MQ=255)
                              ccAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAg                                         >  1:2606645/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)
                                     cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                                                  <  1:3360299/38‑1 (MQ=255)
                                     cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                 <  1:1212140/71‑1 (MQ=255)
                                      gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                >  1:2676850/1‑71 (MQ=255)
                                      gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                                >  1:3674753/1‑71 (MQ=255)
                                           cAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCAccc                           >  1:1171111/1‑71 (MQ=255)
                                            ggtgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCt                                                       <  1:1502732/41‑1 (MQ=37)
                                             gtgtCTTTAATCGACTCTTTATGACCAATCTGGCTGcc                                                          >  1:80866/1‑38 (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)
                                                tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCga                                               >  1:1111400/1‑46 (MQ=255)
                                                tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                 >  1:135462/1‑60 (MQ=255)
                                                tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                 >  1:140120/1‑60 (MQ=255)
                                                tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGa                                 >  1:3227485/1‑60 (MQ=255)
                                                  tttAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                                       <  1:1706149/36‑1 (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)
                                                         gACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgag                                              <  1:1462822/38‑1 (MQ=255)
                                                          aCTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                                            <  1:889712/39‑1 (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: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)
                                                                    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)
                                                                    |                                                                      
CCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATGCGGCAACTTCGAAAG  >  minE/2919373‑2919511

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: