Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,564 A→G 10.7% T44T (ACT→ACC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5640AG10.7% 177.0 / 11.2 65T44T (ACT→ACCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base G (3/4);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCA  >  minE/2919498‑2919628
                                                                  |                                                                
ggCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGAAGTTTTTACCAGt                                                                <  1:3453768/69‑1 (MQ=255)
     cTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg                                                         <  1:457845/71‑1 (MQ=255)
     cTTCGAAAGAGCAACGGGTACGATTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttc                                                         <  1:2919292/71‑2 (MQ=255)
      ttCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                             <  1:2729500/66‑1 (MQ=255)
      ttCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg                                                         >  1:3770343/1‑70 (MQ=255)
        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                       >  1:1068263/1‑67 (MQ=255)
        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                       >  1:392637/1‑67 (MQ=255)
        cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttct                                                        <  1:1149638/69‑3 (MQ=255)
         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGctt                                                     <  1:952633/71‑1 (MQ=255)
         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                      <  1:1519431/70‑1 (MQ=255)
           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                   <  1:467653/71‑1 (MQ=255)
           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                   <  1:3823534/71‑1 (MQ=255)
           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttc                                                    <  1:3578458/70‑1 (MQ=255)
           aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGACTTTGGcttct                                                   >  1:3374360/1‑71 (MQ=255)
            agagCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgc                                                            <  1:905995/61‑2 (MQ=255)
             gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                          <  1:100857/62‑1 (MQ=255)
               gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                          <  1:2674380/60‑1 (MQ=255)
               gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTCTACCAGTGAGTTTGGcttcttctt                                               <  1:1939964/71‑1 (MQ=255)
                cATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCatt                                              <  1:3717555/71‑1 (MQ=255)
                         cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                      <  1:3155227/70‑1 (MQ=255)
                          gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGt                                                            <  1:2125763/47‑1 (MQ=255)
                          gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                       <  1:3854496/68‑1 (MQ=255)
                            gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                  <  1:1232633/71‑1 (MQ=255)
                            gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                   >  1:255071/1‑70 (MQ=255)
                             tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTa                                 >  1:3780291/1‑71 (MQ=255)
                               gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATc                               <  1:501393/71‑1 (MQ=255)
                               gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                               <  1:4054159/71‑1 (MQ=255)
                                agtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTg                                                               <  1:899249/38‑1 (MQ=255)
                                 gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTCTGGCTTCTTCTTTAcc                                           >  1:2137065/1‑57 (MQ=255)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGt                                                                <  1:2239758/35‑1 (MQ=255)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                             <  1:814002/38‑1 (MQ=255)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTAt                                >  1:743782/1‑67 (MQ=255)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGc                            >  1:2023203/1‑71 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                               >  1:246622/1‑50 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttctt                                                  >  1:579742/1‑47 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                          >  1:1809560/1‑71 (MQ=255)
                                    tttttCGAAGATGAGCACGATGTTTTTACCAGTGAGttt                                                          <  1:1955299/39‑1 (MQ=255)
                                    tttttCGAAGATGAGAGCGATGTTTTTACCAGTGAGTTTGGc                                                       >  1:1444758/1‑42 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCGGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTAt                                >  1:736587/1‑63 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                         <  1:3601709/70‑1 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                           <  1:2086482/68‑1 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                           <  1:246912/68‑1 (MQ=255)
                                        tCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                       <  1:1310845/70‑1 (MQ=255)
                                         cGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGc                     >  1:3073001/1‑71 (MQ=255)
                                           aaGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                    >  1:1383695/1‑54 (MQ=255)
                                           aaGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                    >  1:1046343/1‑54 (MQ=255)
                                             gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                 >  1:3248903/1‑71 (MQ=255)
                                              gtGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCgg                             <  1:179334/57‑1 (MQ=255)
                                              aTGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                               <  1:826621/40‑1 (MQ=255)
                                              aTGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                 <  1:3021790/38‑1 (MQ=255)
                                              aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                <  1:1640929/71‑1 (MQ=255)
                                                 aGCGCGATGTTTTTACCAGTGAGTATGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGc             >  1:1765994/1‑71 (MQ=255)
                                                  gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt            <  1:3638689/71‑1 (MQ=255)
                                                  gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt            <  1:976836/71‑1 (MQ=255)
                                                   cgcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa           <  1:3166401/71‑1 (MQ=255)
                                                    gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                             >  1:343917/1‑36 (MQ=255)
                                                    gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                             >  1:100231/1‑36 (MQ=255)
                                                       aTGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTg       <  1:1849651/71‑1 (MQ=255)
                                                       aTGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTg       <  1:95183/71‑1 (MQ=255)
                                                         gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTaa          >  1:3958553/1‑66 (MQ=255)
                                                           ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag    >  1:498776/1‑70 (MQ=255)
                                                           ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag    >  1:1348143/1‑70 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                   >  1:594043/1‑54 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                   >  1:1958947/1‑54 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag    >  1:3199911/1‑69 (MQ=255)
                                                                  aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca  <  1:2736414/65‑1 (MQ=255)
                                                                  aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca  <  1:2509406/65‑1 (MQ=255)
                                                                  aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca  <  1:2387988/65‑1 (MQ=255)
                                                                  aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca  <  1:1825216/65‑1 (MQ=255)
                                                                  aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca  <  1:1192689/65‑1 (MQ=255)
                                                                  |                                                                
GGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCA  >  minE/2919498‑2919628

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: