Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,251 A→G 10.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,2510AG10.7% 227.9 / 16.2 84L149L (TTG→CTG) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (45/30);  new base G (4/5);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCA  >  minE/2919181‑2919320
                                                                      |                                                                     
gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                                       <  1:4103696/71‑1 (MQ=255)
   ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                                    >  1:110867/1‑71 (MQ=255)
   ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                                    >  1:1127789/1‑71 (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)
                     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:789079/1‑71 (MQ=255)
                          aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                                             >  1:2325574/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:2004264/71‑1 (MQ=255)
                                    aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                   >  1:2004543/1‑71 (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:1871423/71‑1 (MQ=255)
                                        gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg                               >  1:1666852/1‑71 (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:3043588/1‑42 (MQ=255)
                                                 cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                   >  1:733461/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:1358819/1‑71 (MQ=255)
                                                     gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                  >  1:2388981/1‑71 (MQ=255)
                                                      aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                       <  1:2860077/65‑1 (MQ=255)
                                                       aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                    >  1:2311023/1‑67 (MQ=255)
                                                         gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTTGGATGGAACTCATTGGt              <  1:932480/71‑1 (MQ=255)
                                                          cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                   <  1:2773059/65‑1 (MQ=255)
                                                          cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg               >  1:3077225/1‑69 (MQ=255)
                                                          cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTg                >  1:1958010/1‑68 (MQ=255)
                                                           tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt              >  1:3385894/1‑69 (MQ=255)
                                                            ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg           >  1:1611566/1‑71 (MQ=255)
                                                             tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa            <  1:974414/69‑1 (MQ=255)
                                                             tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg          >  1:2211170/1‑71 (MQ=255)
                                                               ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc        <  1:757283/71‑1 (MQ=255)
                                                                cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa            >  1:3480953/1‑66 (MQ=255)
                                                                cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa            >  1:2964642/1‑66 (MQ=255)
                                                                  ggCAAATGCTCCTGCATGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    <  1:545322/71‑1 (MQ=255)
                                                                   gCAAATGCTCCTGCATTGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    <  1:2475523/71‑1 (MQ=255)
                                                                   gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    <  1:2361766/71‑1 (MQ=255)
                                                                   gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    <  1:3169777/71‑1 (MQ=255)
                                                                    cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    >  1:2085531/1‑70 (MQ=255)
                                                                     aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTg                >  1:3569019/1‑57 (MQ=255)
                                                                     aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc   <  1:1389860/70‑1 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                         <  1:1632614/47‑1 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                         <  1:2956864/47‑1 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                         <  1:582538/47‑1 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa  >  1:1553601/1‑70 (MQ=255)
                                                                      aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa  <  1:2202348/70‑1 (MQ=255)
                                                                      |                                                                     
GAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCA  >  minE/2919181‑2919320

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: