Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,270 A→C 14.2% 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,2700AC14.2% 180.6 / 26.4 84L142L (CTT→CTGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base C (4/8);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919205‑2919339
                                                                 |                                                                     
tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc                                                                  <  1:984147/71‑1 (MQ=255)
    cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                                                    <  1:1332913/65‑1 (MQ=255)
      gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:2729199/1‑71 (MQ=255)
      gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:3375487/1‑71 (MQ=255)
      gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:1443836/1‑71 (MQ=255)
       cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            >  1:215778/1‑70 (MQ=255)
       cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                            <  1:2524156/70‑1 (MQ=255)
         tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                         >  1:3820175/1‑71 (MQ=255)
         tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                         >  1:1566546/1‑71 (MQ=255)
          aGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                                                        >  1:345191/1‑71 (MQ=255)
           gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:488704/71‑1 (MQ=255)
           gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:1845784/71‑1 (MQ=255)
           gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                       <  1:678766/71‑1 (MQ=255)
           gACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       >  1:2324041/1‑71 (MQ=255)
            aCCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                                          <  1:1536491/67‑1 (MQ=255)
             ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                                     <  1:417966/71‑1 (MQ=255)
              cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                                      >  1:3417191/1‑53 (MQ=255)
              cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                                    >  1:898756/1‑71 (MQ=255)
                   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:1516533/71‑1 (MQ=255)
                      tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            <  1:291035/71‑1 (MQ=255)
                      tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            <  1:2940543/71‑1 (MQ=255)
                         cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                         >  1:464517/1‑71 (MQ=255)
                         cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                         >  1:3500525/1‑71 (MQ=255)
                          gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                          <  1:3462692/69‑1 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                            >  1:16463/1‑66 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                       <  1:1014133/71‑1 (MQ=255)
                           ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                       <  1:1671315/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)
                               aCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTgg                                                     >  1:705940/1‑52 (MQ=255)
                                cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                                   >  1:1315962/1‑38 (MQ=255)
                                 gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                 <  1:2727425/71‑1 (MQ=255)
                                 gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                 <  1:2405156/71‑1 (MQ=255)
                                 gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                                 <  1:56103/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:1081914/1‑71 (MQ=255)
                                      gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                            >  1:356765/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:2438395/71‑1 (MQ=255)
                                              aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                    <  1:1589273/71‑1 (MQ=255)
                                               gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       <  1:2163477/34‑1 (MQ=16)
                                               gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                       <  1:2551466/34‑1 (MQ=16)
                                               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:4085/56‑1 (MQ=255)
                                                   tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                              <  1:3008816/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:254901/1‑71 (MQ=255)
                                                       gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg           >  1:1567680/1‑71 (MQ=255)
                                                        cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc          <  1:1445464/71‑1 (MQ=255)
                                                        cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc          <  1:985564/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)
                                                                 |                                                                     
TCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCGTAT  >  minE/2919205‑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: