Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,234 G→A 14.4% F154F (TTC→TTT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2340GA14.4% 240.2 / 27.2 90F154F (TTC→TTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (44/33);  new base A (7/6);  total (51/39)
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.65e-01

CGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAT  >  minE/2919165‑2919304
                                                                     |                                                                      
cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                                                                       <  1:2055691/71‑1 (MQ=255)
cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                                                                       <  1:1296743/71‑1 (MQ=255)
    gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc                                                                   >  1:866859/1‑71 (MQ=255)
       tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                                                >  1:2638871/1‑71 (MQ=255)
        ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg                                                               >  1:2298655/3‑71 (MQ=255)
        ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg                                                               >  1:573383/3‑71 (MQ=255)
        cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                                                <  1:3383346/70‑1 (MQ=255)
         gAGCATCGAATTGCCCATGTTGTTAGGCGTGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                              <  1:1053910/71‑1 (MQ=255)
         gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                              >  1:3029943/1‑71 (MQ=255)
          aGCATCGAATTGCCCATGTTGTTAAGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                                                             <  1:2289205/71‑1 (MQ=255)
           gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg                                                            <  1:902106/71‑1 (MQ=255)
            cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                              <  1:3393915/68‑1 (MQ=255)
              tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                                                        <  1:1649990/56‑1 (MQ=255)
              tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                                                        <  1:3627368/56‑1 (MQ=255)
                gagttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                                        <  1:1098310/67‑1 (MQ=255)
                gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                       <  1:4103696/71‑1 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGaa                                                                      >  1:547197/1‑56 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        <  1:2921634/70‑1 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        <  1:132593/70‑1 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        <  1:598495/70‑1 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        >  1:481889/1‑70 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        >  1:2397820/1‑70 (MQ=255)
                gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                        >  1:93207/1‑70 (MQ=255)
                  aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGc                                                         >  1:1354108/1‑67 (MQ=255)
                   ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                    >  1:1127789/1‑71 (MQ=255)
                   ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                    >  1:110867/1‑71 (MQ=255)
                     gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                           >  1:635454/1‑62 (MQ=255)
                     gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                           >  1:2797244/1‑62 (MQ=255)
                         aTGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg                                              <  1:2449036/71‑1 (MQ=255)
                         aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                   <  1:3455239/66‑1 (MQ=255)
                         aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                              <  1:2386096/71‑1 (MQ=255)
                          tgttgtTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg                                              >  1:2568003/1‑70 (MQ=255)
                          tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                              >  1:1100509/1‑70 (MQ=255)
                           gttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCt                                               >  1:398805/1‑68 (MQ=255)
                            ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                 >  1:1663463/1‑65 (MQ=255)
                             tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                            >  1:359370/1‑69 (MQ=255)
                             tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                            >  1:364041/1‑69 (MQ=255)
                             tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                          >  1:2263890/1‑71 (MQ=255)
                              gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATgg                                         >  1:94299/1‑71 (MQ=255)
                                tACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                                        >  1:2504116/1‑70 (MQ=255)
                                 aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                            >  1:2431715/1‑65 (MQ=255)
                                  cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                     <  1:1267869/71‑1 (MQ=255)
                                  cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                     <  1:1368559/71‑1 (MQ=255)
                                  cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtcat                                     <  1:2559838/71‑4 (MQ=255)
                                  cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtcat                                     <  1:649083/71‑4 (MQ=255)
                                   gcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                              <  1:2697759/61‑1 (MQ=255)
                                     gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                                                               <  1:915770/42‑1 (MQ=255)
                                     gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                                                               <  1:3266232/42‑1 (MQ=255)
                                     gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                     >  1:41487/1‑68 (MQ=255)
                                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                      >  1:3445416/1‑46 (MQ=255)
                                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                              <  1:3590668/70‑1 (MQ=255)
                                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                              <  1:331130/70‑1 (MQ=255)
                                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                             >  1:2325574/1‑71 (MQ=255)
                                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                             >  1:789079/1‑71 (MQ=255)
                                            cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGcc                           <  1:1040683/71‑1 (MQ=255)
                                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                    >  1:3395957/1‑44 (MQ=255)
                                              gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                         >  1:1286198/1‑71 (MQ=255)
                                               cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                         >  1:4013411/1‑70 (MQ=255)
                                               cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                         >  1:3527105/1‑70 (MQ=255)
                                                gtAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                       <  1:2285979/70‑1 (MQ=255)
                                                  aGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                     >  1:2015197/1‑71 (MQ=255)
                                                    aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                    >  1:2276215/1‑38 (MQ=255)
                                                    aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                   >  1:2004543/1‑71 (MQ=255)
                                                    aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                   <  1:2004264/71‑1 (MQ=255)
                                                      cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                 <  1:85451/71‑1 (MQ=255)
                                                      cAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                          <  1:3059441/46‑1 (MQ=255)
                                                        gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                           <  1:2488407/43‑1 (MQ=255)
                                                        gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg               >  1:1666852/1‑71 (MQ=255)
                                                        gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg               <  1:1871423/71‑1 (MQ=255)
                                                        gCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGt                                        >  1:1190934/1‑46 (MQ=255)
                                                         cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg              <  1:2751769/71‑1 (MQ=255)
                                                         cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg              >  1:789520/1‑71 (MQ=255)
                                                          gTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg              >  1:3600264/1‑70 (MQ=255)
                                                           tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                              <  1:680515/37‑1 (MQ=255)
                                                             aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCGTGCa                                            >  1:1390525/1‑37 (MQ=255)
                                                             aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg          >  1:990945/1‑71 (MQ=255)
                                                             aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg          <  1:3168188/71‑1 (MQ=255)
                                                             aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTCGGATg          >  1:2055453/1‑71 (MQ=255)
                                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                 <  1:497683/47‑1 (MQ=255)
                                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaag                                  <  1:307900/46‑1 (MQ=255)
                                                              tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg         <  1:2113160/71‑1 (MQ=255)
                                                                tCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt           >  1:2744971/1‑67 (MQ=255)
                                                                tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa       >  1:2685177/1‑71 (MQ=255)
                                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                   >  1:733461/1‑42 (MQ=255)
                                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                   >  1:3043588/1‑42 (MQ=255)
                                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                   >  1:485610/1‑42 (MQ=255)
                                                                  gTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                   >  1:956752/1‑41 (MQ=255)
                                                                  gTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg              >  1:2624806/1‑62 (MQ=255)
                                                                   ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc    <  1:3418106/71‑1 (MQ=255)
                                                                   ttGAACGCTTTGCCGGGCAAATGCTCCTCCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc    <  1:2043968/71‑1 (MQ=255)
                                                                    tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa       >  1:1535395/1‑67 (MQ=255)
                                                                    tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa       >  1:2506129/1‑67 (MQ=255)
                                                                     gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa   <  1:2427816/70‑1 (MQ=255)
                                                                     gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt  >  1:2388981/1‑71 (MQ=255)
                                                                     gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt  >  1:1358819/1‑71 (MQ=255)
                                                                     |                                                                      
CGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAT  >  minE/2919165‑2919304

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