Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,273 A→G 7.7% D141D (GAT→GAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2730AG 7.7% 217.5 / 8.7 78D141D (GAT→GACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base G (4/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919207‑2919339
                                                                  |                                                                  
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)
    gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:1286198/1‑71 (MQ=255)
     cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:3527105/1‑70 (MQ=255)
     cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:4013411/1‑70 (MQ=255)
      gtAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                                          <  1:2285979/70‑1 (MQ=255)
        aGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                                                        >  1:2015197/1‑71 (MQ=255)
          aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                                      >  1:2004543/1‑71 (MQ=255)
          aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                                      <  1:2004264/71‑1 (MQ=255)
            cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                                    <  1:85451/71‑1 (MQ=255)
              gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg                                                  >  1:1666852/1‑71 (MQ=255)
              gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg                                                  <  1:1871423/71‑1 (MQ=255)
               cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                                 <  1:2751769/71‑1 (MQ=255)
               cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                                 >  1:789520/1‑71 (MQ=255)
                gTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                                 >  1:3600264/1‑70 (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)
                    tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            <  1:2113160/71‑1 (MQ=255)
                      tCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                                              >  1:2744971/1‑67 (MQ=255)
                      tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                          >  1:2685177/1‑71 (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)
                            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)
                                           aggtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCa                     >  1:98435/4‑71 (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:2202348/70‑1 (MQ=255)
                                            aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                     >  1:1553601/1‑70 (MQ=255)
                                             aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                                  <  1:228071/56‑1 (MQ=255)
                                             aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                                  <  1:425009/56‑1 (MQ=255)
                                              tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                     >  1:2689462/1‑68 (MQ=255)
                                              tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc                  <  1:1137237/71‑1 (MQ=255)
                                              tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAgcc                            >  1:1982011/1‑59 (MQ=255)
                                               gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                                        <  1:1621294/48‑1 (MQ=255)
                                               gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCATGCCATTCCATAcc                 >  1:647500/1‑71 (MQ=255)
                                                   cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCgg               <  1:54855/69‑1 (MQ=255)
                                                   cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGc              <  1:2539398/70‑1 (MQ=255)
                                                   cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCgg               <  1:1283027/69‑1 (MQ=255)
                                                    tGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            <  1:3690060/39‑1 (MQ=255)
                                                     gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCAc            >  1:3675509/1‑70 (MQ=255)
                                                     gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg           >  1:2395215/1‑71 (MQ=255)
                                                     gCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGtt                                    >  1:1988023/1‑46 (MQ=12)
                                                      cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                          >  1:4136648/1‑39 (MQ=255)
                                                      cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                          >  1:321127/1‑39 (MQ=255)
                                                      cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                          >  1:1828206/1‑39 (MQ=255)
                                                      cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCAc            <  1:1229525/69‑1 (MQ=255)
                                                          gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                          >  1:716041/1‑35 (MQ=255)
                                                          gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                            >  1:2418342/1‑49 (MQ=255)
                                                             agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                     >  1:3191303/1‑53 (MQ=255)
                                                              tcagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt  <  1:3958675/69‑1 (MQ=255)
                                                                  aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt                         <  1:2828547/44‑1 (MQ=255)
                                                                  |                                                                  
ACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919207‑2919339

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: