Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,251 A→G 20.7% L149L (TTG→CTG)  argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2510AG20.7% 175.2 / 41.8 87L149L (TTG→CTG) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (8/10);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAT  >  minE/2919181‑2919321
                                                                      |                                                                      
gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                                        <  1:2893689/71‑1 (MQ=255)
gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                                        <  1:3745365/71‑1 (MQ=255)
gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                                        <  1:762208/71‑1 (MQ=255)
 aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                                       >  1:3741343/1‑71 (MQ=255)
 aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                                       <  1:3235175/71‑1 (MQ=255)
     gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCt                                                                   >  1:554028/1‑71 (MQ=255)
       ccATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTcc                                                                 <  1:3782032/71‑1 (MQ=255)
        cATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCt                                                                >  1:3347436/1‑71 (MQ=255)
         aTGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg                                                               <  1:1265581/71‑1 (MQ=255)
         aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                                    <  1:2935722/66‑1 (MQ=255)
          tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  >  1:2627365/1‑67 (MQ=255)
          tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                                               >  1:793078/1‑70 (MQ=255)
          tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                                               >  1:2860423/1‑70 (MQ=255)
           gttgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                                                             >  1:2923270/1‑71 (MQ=255)
            ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                                                             <  1:1330796/70‑1 (MQ=255)
            ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                                                             <  1:789498/70‑1 (MQ=255)
            ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                                                             <  1:2089736/70‑1 (MQ=255)
            ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCAt                                                            >  1:1159080/1‑71 (MQ=255)
            ttgttACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAAt                                                            >  1:1735195/1‑71 (MQ=255)
             tgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                                           >  1:2304940/1‑71 (MQ=255)
             tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                                            <  1:465540/70‑1 (MQ=255)
              gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATgg                                                          >  1:2865924/1‑71 (MQ=255)
              gttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                                           <  1:2745274/70‑1 (MQ=255)
               ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                                                          <  1:358280/70‑1 (MQ=255)
               ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                                                          <  1:1224030/70‑1 (MQ=255)
               ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                                                         <  1:1745529/71‑1 (MQ=255)
                 aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                                                         >  1:291980/1‑69 (MQ=255)
                 aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTg                                                        <  1:3118263/70‑1 (MQ=255)
                  cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtca                                                       <  1:414890/70‑3 (MQ=255)
                   gcgcgTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                      >  1:3389424/1‑70 (MQ=255)
                        tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc                                                <  1:984147/71‑1 (MQ=255)
                            cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                                  <  1:1332913/65‑1 (MQ=255)
                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                                                          <  1:1371044/55‑1 (MQ=255)
                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                          >  1:1443836/1‑71 (MQ=255)
                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                          >  1:3375487/1‑71 (MQ=255)
                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                          >  1:2729199/1‑71 (MQ=255)
                               cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                          >  1:215778/1‑70 (MQ=255)
                               cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                          <  1:2524156/70‑1 (MQ=255)
                                 tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                       >  1:3820175/1‑71 (MQ=255)
                                 tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                       >  1:1566546/1‑71 (MQ=255)
                                  aGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                                      >  1:345191/1‑71 (MQ=255)
                                   gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     <  1:488704/71‑1 (MQ=255)
                                   gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     <  1:678766/71‑1 (MQ=255)
                                   gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     <  1:1845784/71‑1 (MQ=255)
                                   gACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                     >  1:2324041/1‑71 (MQ=255)
                                    aCCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                        <  1:1536491/67‑1 (MQ=255)
                                     ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                   <  1:417966/71‑1 (MQ=255)
                                      cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                    >  1:3417191/1‑53 (MQ=255)
                                      cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                  >  1:898756/1‑71 (MQ=255)
                                           tcctttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                             <  1:1672665/68‑1 (MQ=255)
                                           tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                             <  1:1495109/71‑1 (MQ=255)
                                            cATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                            >  1:3709877/1‑71 (MQ=255)
                                             atctCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTg                           <  1:2544862/68‑1 (MQ=255)
                                             aTTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                           <  1:1517941/71‑1 (MQ=255)
                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                                           <  1:3224539/38‑1 (MQ=255)
                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGTGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:2800/71‑1 (MQ=255)
                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:291035/71‑1 (MQ=255)
                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:1516533/71‑1 (MQ=255)
                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:2940543/71‑1 (MQ=255)
                                               ttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                             >  1:2885270/1‑35 (MQ=255)
                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                       >  1:464517/1‑71 (MQ=255)
                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                       >  1:3500525/1‑71 (MQ=255)
                                                  gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                        <  1:3462692/69‑1 (MQ=255)
                                                   ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          >  1:16463/1‑66 (MQ=255)
                                                   ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                     <  1:1014133/71‑1 (MQ=255)
                                                   ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                     <  1:1671315/71‑1 (MQ=255)
                                                    tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGTATGGAACTCa                    >  1:242989/1‑71 (MQ=255)
                                                    tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                     <  1:856089/70‑1 (MQ=255)
                                                    tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                    >  1:632586/1‑71 (MQ=255)
                                                     gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCa                    <  1:2928175/70‑1 (MQ=255)
                                                      aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:3392923/63‑1 (MQ=255)
                                                      aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   <  1:935183/70‑1 (MQ=255)
                                                      aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                   <  1:3478967/70‑3 (MQ=255)
                                                      aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                   <  1:1989517/70‑3 (MQ=255)
                                                       aCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTgg                                   >  1:705940/1‑52 (MQ=255)
                                                        cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                 >  1:1315962/1‑38 (MQ=255)
                                                         gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt               <  1:2727425/71‑1 (MQ=255)
                                                         gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt               <  1:2405156/71‑1 (MQ=255)
                                                         gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt               <  1:56103/71‑1 (MQ=255)
                                                           tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt               >  1:25113/1‑69 (MQ=255)
                                                           cttGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTCGAACTCGTTGGt               >  1:260038/2‑69 (MQ=255)
                                                            ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg            >  1:1152492/1‑71 (MQ=255)
                                                             tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg           >  1:983456/1‑71 (MQ=255)
                                                              gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc          >  1:1081914/1‑71 (MQ=255)
                                                              gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc          >  1:356765/1‑71 (MQ=255)
                                                                cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc         <  1:1194635/70‑1 (MQ=255)
                                                                 gggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt       <  1:2248350/71‑1 (MQ=255)
                                                                  ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt      <  1:3875150/71‑1 (MQ=255)
                                                                  ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt      <  1:3517193/71‑1 (MQ=255)
                                                                   gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc     >  1:2506764/1‑71 (MQ=255)
                                                                    cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc     >  1:3647491/1‑70 (MQ=255)
                                                                    cAAATGCTCCTGCATGGTGAGAAGATCCGCAAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc    <  1:2330912/71‑1 (MQ=255)
                                                                     aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc    >  1:1394035/1‑70 (MQ=255)
                                                                     aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa   >  1:3419366/1‑71 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                         >  1:2164908/1‑48 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                         >  1:301951/1‑48 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                  >  1:2429364/1‑55 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt  <  1:1589273/71‑1 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt  <  1:2438395/71‑1 (MQ=255)
                                                                      |                                                                      
GAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAT  >  minE/2919181‑2919321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: