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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2940AG 8.8% 219.6 / 10.4 79H134H (CAT→CACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (2/5);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCT  >  minE/2919224‑2919364
                                                                      |                                                                      
tcctttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                        <  1:1672665/68‑1 (MQ=255)
tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                                                        <  1:1495109/71‑1 (MQ=255)
 cATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                                                                       >  1:3709877/1‑71 (MQ=255)
  atctCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTg                                                                      <  1:2544862/68‑1 (MQ=255)
  aTTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                      <  1:1517941/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGTGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     <  1:2800/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     <  1:291035/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     <  1:1516533/71‑1 (MQ=255)
   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     <  1:2940543/71‑1 (MQ=255)
      cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                                  >  1:3500525/1‑71 (MQ=255)
      cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                                                  >  1:464517/1‑71 (MQ=255)
       gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                                   <  1:3462692/69‑1 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     >  1:16463/1‑66 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                                <  1:1671315/71‑1 (MQ=255)
        ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                                <  1:1014133/71‑1 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGTATGGAACTCa                                                               >  1:242989/1‑71 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                                                <  1:856089/70‑1 (MQ=255)
         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                                               >  1:632586/1‑71 (MQ=255)
          gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCa                                                               <  1:2928175/70‑1 (MQ=255)
           aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                                                     <  1:3392923/63‑1 (MQ=255)
           aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                              <  1:935183/70‑1 (MQ=255)
           aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                              <  1:1989517/70‑3 (MQ=255)
           aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                                              <  1:3478967/70‑3 (MQ=255)
              gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                          <  1:2727425/71‑1 (MQ=255)
              gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                          <  1:56103/71‑1 (MQ=255)
              gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                          <  1:2405156/71‑1 (MQ=255)
                tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                                          >  1:25113/1‑69 (MQ=255)
                cttGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTCGAACTCGTTGGt                                                          >  1:260038/2‑69 (MQ=255)
                 ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                       >  1:1152492/1‑71 (MQ=255)
                  tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                                                      >  1:983456/1‑71 (MQ=255)
                   gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                     >  1:356765/1‑71 (MQ=255)
                   gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                                                     >  1:1081914/1‑71 (MQ=255)
                     cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc                                                    <  1:1194635/70‑1 (MQ=255)
                      gggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt                                                  <  1:2248350/71‑1 (MQ=255)
                       ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                                 <  1:3517193/71‑1 (MQ=255)
                       ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                                                 <  1:3875150/71‑1 (MQ=255)
                        gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                                >  1:2506764/1‑71 (MQ=255)
                         cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                                                >  1:3647491/1‑70 (MQ=255)
                         cAAATGCTCCTGCATGGTGAGAAGATCCGCAAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                               <  1:2330912/71‑1 (MQ=255)
                          aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                                               >  1:1394035/1‑70 (MQ=255)
                          aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                                              >  1:3419366/1‑71 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                                    >  1:2164908/1‑48 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                                                    >  1:301951/1‑48 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                                             >  1:2429364/1‑55 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                             <  1:1589273/71‑1 (MQ=255)
                           aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                                             <  1:2438395/71‑1 (MQ=255)
                            aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATa                                            <  1:738670/71‑1 (MQ=255)
                             tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                                                   >  1:278211/1‑47 (MQ=16)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                                                      >  1:3463126/1‑43 (MQ=255)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                          >  1:1154139/1‑71 (MQ=255)
                              gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc                                          >  1:2131374/1‑71 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGGTGCGTGGGATGGAACTCATTGGTCAgg                                                      >  1:1842506/1‑57 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                       <  1:3008816/56‑1 (MQ=255)
                                tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                                                       <  1:4085/56‑1 (MQ=255)
                                 ccTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCATTGGTCAg                                                       <  1:3099690/55‑1 (MQ=255)
                                  cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa                                      >  1:2599004/1‑71 (MQ=255)
                                  cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGGCAGGCCATTCCATACCgg                                        <  1:2504979/69‑1 (MQ=255)
                                  cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGc                                       <  1:2359793/70‑1 (MQ=255)
                                    gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                                    >  1:1567680/1‑71 (MQ=255)
                                    gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg                                    >  1:254901/1‑71 (MQ=255)
                                     cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                                   <  1:985564/71‑1 (MQ=255)
                                     cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                                   <  1:1445464/71‑1 (MQ=255)
                                      aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc                                   >  1:1852288/1‑70 (MQ=255)
                                          tGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCg                              <  1:236448/71‑1 (MQ=255)
                                            agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                                              <  1:2142807/37‑1 (MQ=255)
                                            agaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTa                            <  1:3870150/71‑1 (MQ=255)
                                             gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                           >  1:2488977/1‑71 (MQ=255)
                                             gaagaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                           <  1:1442393/71‑1 (MQ=255)
                                               agaTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTc                         >  1:2220404/1‑71 (MQ=255)
                                                 aTCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTc                                                         <  1:627250/37‑1 (MQ=255)
                                                  tCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGc                      >  1:2926612/1‑71 (MQ=255)
                                                    cGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCa                    <  1:2358888/71‑1 (MQ=255)
                                                       cagcagCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGcc                     <  1:1370370/67‑1 (MQ=255)
                                                          cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt              >  1:2596451/1‑71 (MQ=255)
                                                          cagcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt              >  1:2651682/1‑71 (MQ=255)
                                                           agcTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTc             <  1:568277/71‑1 (MQ=255)
                                                             cTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGttt              <  1:2716573/68‑1 (MQ=255)
                                                                cGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACaa        >  1:3275705/1‑71 (MQ=255)
                                                                agTGGGATGGAACTCATTGGTCAGGCCATTCCATACAGGCACGCTAGCGt                             >  1:3782737/2‑50 (MQ=255)
                                                                 gTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAt                           >  1:3475825/1‑51 (MQ=255)
                                                                   gggATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAAt       >  1:1176296/1‑69 (MQ=255)
                                                                      aTGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt  >  1:263066/1‑71 (MQ=255)
                                                                      aTGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCt  >  1:934276/1‑71 (MQ=255)
                                                                      |                                                                      
TCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTATTCCGCCAGTGTTTCGACAATCTCCT  >  minE/2919224‑2919364

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: