Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,567 G→A 8.3% L43L (CTC→CTT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5670GA 8.3% 226.2 / 7.2 72L43L (CTC→CTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (25/41);  new base A (2/4);  total (27/45)
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

CTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTT  >  minE/2919503‑2919635
                                                                |                                                                    
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)
                             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:246912/68‑1 (MQ=255)
                                 tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                  <  1:2086482/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:976836/71‑1 (MQ=255)
                                             gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                   <  1:3638689/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:95183/71‑1 (MQ=255)
                                                  aTGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTg              <  1:1849651/71‑1 (MQ=255)
                                                    gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTaa                 >  1:3958553/1‑66 (MQ=255)
                                                      ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag           >  1:1348143/1‑70 (MQ=255)
                                                      ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag           >  1:498776/1‑70 (MQ=255)
                                                       tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                          >  1:1958947/1‑54 (MQ=255)
                                                       tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                          >  1:594043/1‑54 (MQ=255)
                                                       tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag           >  1:3199911/1‑69 (MQ=255)
                                                             aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca         <  1:2387988/65‑1 (MQ=255)
                                                             aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca         <  1:2509406/65‑1 (MQ=255)
                                                             aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca         <  1:1825216/65‑1 (MQ=255)
                                                             aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca         <  1:2736414/65‑1 (MQ=255)
                                                             aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagca         <  1:1192689/65‑1 (MQ=255)
                                                              gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt  <  1:2475526/71‑1 (MQ=255)
                                                              gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt  <  1:1268207/71‑1 (MQ=255)
                                                              gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtc  <  1:1120670/71‑2 (MQ=255)
                                                               tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCgg                     <  1:3791958/51‑1 (MQ=255)
                                                               tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                   <  1:280414/53‑1 (MQ=255)
                                                                gAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAg                             >  1:1109434/1‑42 (MQ=255)
                                                                |                                                                    
CTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTT  >  minE/2919503‑2919635

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: