Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,429 C→T 6.8% S89S (TCG→TCA argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,4290CT 6.8% 236.4 / 4.9 73S89S (TCG→TCAargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (41/27);  new base T (2/3);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATG  >  minE/2919360‑2919496
                                                                     |                                                                   
ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                     >  1:597122/1‑70 (MQ=255)
ctcCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                     >  1:4052687/1‑70 (MQ=255)
   cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                     <  1:685844/67‑1 (MQ=255)
   cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                                                   >  1:2786790/1‑69 (MQ=255)
   cTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                                                   <  1:1358783/69‑1 (MQ=255)
     gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                                                               >  1:732057/1‑71 (MQ=255)
     gACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                                                               >  1:1004185/1‑71 (MQ=255)
      aCCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                                                                <  1:209558/69‑1 (MQ=255)
        cATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActct                                                               <  1:1241767/68‑1 (MQ=255)
         aTAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGGGGGCAGTGTCTTTAATCGACTCTTTAt                                                           <  1:3481903/71‑1 (MQ=255)
          tAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGActc                                                                <  1:3852490/65‑1 (MQ=255)
          tAGCCGCGATACTGAATACCGTCATACATGCGACAAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTAt                                                           >  1:355954/1‑70 (MQ=255)
           aGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                                         >  1:1128770/1‑71 (MQ=255)
            gCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGAc                                                        >  1:2593185/1‑71 (MQ=255)
              cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                                                 >  1:1069099/1‑60 (MQ=255)
              cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                                                 >  1:3056088/1‑60 (MQ=255)
              cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGAct                                                                 >  1:1255088/1‑60 (MQ=255)
              cgcgATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCa                                                      <  1:1120459/71‑1 (MQ=255)
                   tACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCtt                                                              <  1:4109202/58‑1 (MQ=255)
                    aatGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                     >  1:2421525/3‑50 (MQ=255)
                    aCTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                                   >  1:1156875/1‑68 (MQ=255)
                     cTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTgg                                                <  1:413946/70‑1 (MQ=255)
                       gAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTg                                             >  1:1953984/1‑71 (MQ=255)
                        aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGa                                                                   >  1:1581934/1‑48 (MQ=255)
                        aaTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGc                                            >  1:3935998/1‑71 (MQ=255)
                         aTACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATc                                                                     <  1:4037096/45‑1 (MQ=255)
                          tgccGTCATACATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCg                                          <  1:482992/69‑1 (MQ=255)
                          tACCGTCATACATGCGACCAAGAACGCGGGCAGTGTCTTTAATCGa                                                                   >  1:1483838/1‑46 (MQ=255)
                            ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGc                                         <  1:3018203/70‑1 (MQ=255)
                            ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                        <  1:1804524/71‑1 (MQ=255)
                            ccGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCt                                        <  1:2335345/71‑1 (MQ=255)
                              gTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCttt                                                             <  1:254613/48‑1 (MQ=255)
                               tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                                   <  1:3181627/57‑1 (MQ=255)
                               tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                                   <  1:1763417/57‑1 (MQ=255)
                               tCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATc                                                   <  1:414439/57‑1 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGcct                                          >  1:2081716/1‑63 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc                                   >  1:593064/1‑71 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc                                   >  1:3359434/1‑71 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc                                   >  1:3858920/1‑71 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc                                   >  1:2231304/1‑71 (MQ=255)
                                 aTACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGc                                   >  1:465663/1‑71 (MQ=255)
                                  tACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGcc                                  >  1:1527442/1‑71 (MQ=255)
                                    cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                                            >  1:2087736/1‑43 (MQ=255)
                                    cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTa                                                            >  1:2082932/1‑43 (MQ=255)
                                    cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGa                                                         >  1:3082744/1‑46 (MQ=255)
                                    cATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGc                                               >  1:1975978/1‑56 (MQ=255)
                                    cATCCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGccct                                 <  1:1830262/70‑2 (MQ=255)
                                       gCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCgag                               >  1:175127/1‑69 (MQ=255)
                                        cGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTg                                                 >  1:2184810/1‑50 (MQ=255)
                                          aCCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCa                                                      >  1:3903239/1‑43 (MQ=255)
                                            cAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                         <  1:3181648/70‑1 (MQ=255)
                                             aaGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                       >  1:2889833/1‑71 (MQ=255)
                                              aGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCga                                <  1:1348478/61‑1 (MQ=255)
                                                 aCGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAAcgc                    <  1:2266016/70‑1 (MQ=255)
                                                  cgcgGGCAGTGTCTTTAATCGACTCTTTATGACCa                                                      <  1:3883851/35‑1 (MQ=255)
                                                   gcgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                 <  1:3278124/71‑1 (MQ=255)
                                                    cgGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGAt                             <  1:1477739/58‑1 (MQ=255)
                                                       gCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCAcc             >  1:1359381/1‑71 (MQ=255)
                                                           tgtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGc                                         <  1:3730096/39‑1 (MQ=255)
                                                            gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAg                 >  1:3588813/1‑62 (MQ=255)
                                                            gtCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCa               >  1:3581998/1‑64 (MQ=255)
                                                             tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCtt                                       >  1:2397848/1‑39 (MQ=255)
                                                             tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCTAGATAAGTAACGCGAGCACCCTGGt        <  1:1427368/70‑1 (MQ=255)
                                                             tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                         <  1:1909759/53‑1 (MQ=255)
                                                             tCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTaa                       <  1:1498458/55‑1 (MQ=255)
                                                              cTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                           >  1:2875433/1‑50 (MQ=255)
                                                              cTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                           >  1:2176606/1‑50 (MQ=255)
                                                              cTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATaa                           >  1:1111547/1‑50 (MQ=255)
                                                              cTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTgg         <  1:2384274/68‑1 (MQ=255)
                                                               tttAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCa      >  1:3471144/1‑70 (MQ=255)
                                                               cttAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCa      >  1:777738/2‑70 (MQ=255)
                                                                 tAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTaa                       >  1:2203709/1‑51 (MQ=255)
                                                                 tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGATATAAGTAACGCGAGCACCCTGGTCatat   >  1:1293599/1‑71 (MQ=255)
                                                                 tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCatat   >  1:1165112/1‑71 (MQ=255)
                                                                 tAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCata    >  1:150592/1‑70 (MQ=255)
                                                                  aaTTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCaaat   <  1:2411471/70‑1 (MQ=255)
                                                                  aaTCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATg  >  1:968957/1‑71 (MQ=255)
                                                                    tCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGt                         <  1:639508/46‑1 (MQ=255)
                                                                     cGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCt           <  1:1521693/59‑1 (MQ=255)
                                                                     |                                                                   
CTCCTGACCATAGCCGCGATACTGAATACCGTCATACATGCGACCAAGCACGCGGGCAGTGTCTTTAATCGACTCTTTATGACCAATCTGGCTGCCGCTTGGGCCGAGATAAGTAACGCGAGCACCCTGGTCATATG  >  minE/2919360‑2919496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: