Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,252 A→G 10.7% H148H (CAT→CAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2520AG10.7% 228.3 / 17.4 84H148H (CAT→CACargIornithine 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.75e-01

TTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCA  >  minE/2919184‑2919320
                                                                    |                                                                    
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:1368559/71‑1 (MQ=255)
               cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                     <  1:1267869/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:331130/70‑1 (MQ=255)
                       aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                              <  1:3590668/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: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:3168188/71‑1 (MQ=255)
                                          aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                          >  1:990945/1‑71 (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: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:3169777/71‑1 (MQ=255)
                                                                gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc    <  1:2361766/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:582538/47‑1 (MQ=255)
                                                                   aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                         <  1:1632614/47‑1 (MQ=255)
                                                                   aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                         <  1:2956864/47‑1 (MQ=255)
                                                                   aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa  >  1:1553601/1‑70 (MQ=255)
                                                                   aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa  <  1:2202348/70‑1 (MQ=255)
                                                                    aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg               <  1:425009/56‑1 (MQ=255)
                                                                    aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg               <  1:228071/56‑1 (MQ=255)
                                                                    |                                                                    
TTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCA  >  minE/2919184‑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: