Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,270 A→C 7.7% L142L (CTT→CTG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2700AC 7.7% 197.6 / 13.1 78L142L (CTT→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (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.66e-01

ACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACG  >  minE/2919207‑2919330
                                                               |                                                            
aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                                        <  1:331130/70‑1 (MQ=255)
aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                                        <  1:3590668/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:789520/1‑71 (MQ=255)
               cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                        <  1:2751769/71‑1 (MQ=255)
                gTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                        >  1:3600264/1‑70 (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:485610/1‑42 (MQ=255)
                       cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                             >  1:3043588/1‑42 (MQ=255)
                       cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                             >  1:733461/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)
                            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:2964642/1‑66 (MQ=255)
                                      cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa                      >  1:3480953/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: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: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:321127/1‑39 (MQ=255)
                                                      cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                 >  1:4136648/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)
                                                               |                                                            
ACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACG  >  minE/2919207‑2919330

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: