Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,564 A→G 21.2% 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,5640AG21.2% 195.6 / 49.0 99T44T (ACT→ACCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (30/48);  new base G (7/14);  total (37/62)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTG  >  minE/2919500‑2919633
                                                                |                                                                     
cAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                                  <  1:573463/70‑1 (MQ=255)
cAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGt                                                                 <  1:2060975/71‑1 (MQ=255)
cAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgc                                                                 <  1:1477191/71‑2 (MQ=255)
 aaCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                <  1:2673147/71‑1 (MQ=255)
 aaCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgct                                                                <  1:3128023/71‑3 (MQ=255)
 aaCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgct                                                                >  1:261907/1‑69 (MQ=255)
 aaCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTATCGAAGATGA‑‑GCGATGTTTTTACCGGTAAGcttc                                                              >  1:1687322/1‑70 (MQ=255)
  aCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                               >  1:3420171/1‑71 (MQ=255)
  aCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                               <  1:2970049/71‑1 (MQ=255)
  aCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                               <  1:823246/71‑1 (MQ=255)
      cGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTgg                                                             <  1:2693425/69‑1 (MQ=255)
      cGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                           >  1:3577608/1‑71 (MQ=255)
      cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                            <  1:963201/70‑4 (MQ=255)
      cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                            >  1:1750826/1‑67 (MQ=255)
      cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                            >  1:3855338/1‑67 (MQ=255)
      cGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                            <  1:2668452/70‑4 (MQ=255)
         aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGt                                                                     <  1:1426296/58‑1 (MQ=255)
         aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                        >  1:1460845/1‑71 (MQ=255)
           gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                      <  1:1332095/71‑1 (MQ=255)
           gagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                      <  1:2143960/71‑1 (MQ=255)
            agCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttct                                                     >  1:2330690/1‑71 (MQ=255)
            agCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttct                                                     <  1:3049360/71‑1 (MQ=255)
             gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                    >  1:2547126/1‑71 (MQ=255)
              cATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                    <  1:2522935/70‑1 (MQ=255)
                     tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                            <  1:2555359/71‑1 (MQ=255)
                     tACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                            <  1:3063819/71‑1 (MQ=255)
                       cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                           <  1:1964499/70‑1 (MQ=255)
                       cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                           >  1:1144197/1‑70 (MQ=255)
                        gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGt                                                                 <  1:1101406/47‑1 (MQ=255)
                        gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                    <  1:2133161/60‑1 (MQ=255)
                          gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                       <  1:3491272/71‑1 (MQ=255)
                          gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                       <  1:156901/71‑1 (MQ=255)
                          gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                       <  1:458861/71‑1 (MQ=255)
                          gtcgagtcTTTTTCGAAGATGAGCCCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                        >  1:1430430/1‑70 (MQ=255)
                           tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTa                                      <  1:2460278/71‑1 (MQ=255)
                            cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                               >  1:3115051/1‑45 (MQ=255)
                            cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                            <  1:932425/64‑1 (MQ=255)
                             gagtcTTTTTGGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATc                                    >  1:2406720/1‑71 (MQ=255)
                             gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGATTCTGTACTTCCTTGCCATTTTTTTTATc                                    >  1:1382560/1‑71 (MQ=255)
                             gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTAt                                     <  1:1185952/70‑1 (MQ=255)
                             gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                    <  1:917204/71‑1 (MQ=255)
                               gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCg                                   >  1:737502/1‑70 (MQ=255)
                               gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCgg                                  <  1:1258608/71‑1 (MQ=255)
                                tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGta                                                                    <  1:345880/36‑2 (MQ=37)
                                tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                 <  1:3193547/71‑1 (MQ=255)
                                tcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                      >  1:3147016/1‑50 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTCCCAGTGAGttt                                                               >  1:2351469/1‑39 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                               >  1:2958424/1‑39 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                >  1:2538324/1‑38 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                    <  1:1690350/50‑1 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                        <  1:2084985/46‑1 (MQ=255)
                                  tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                         >  1:1812235/1‑61 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                               <  1:3306119/69‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                               <  1:2106005/69‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                <  1:260631/68‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                <  1:2574720/68‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                <  1:3732892/68‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                           <  1:1047821/41‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                              <  1:1167246/70‑1 (MQ=255)
                                    tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                              <  1:1762441/70‑1 (MQ=255)
                                     ttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttg                             <  1:875383/70‑2 (MQ=255)
                                     ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                             <  1:918403/70‑1 (MQ=255)
                                      tCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                           <  1:963369/39‑1 (MQ=255)
                                        gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                           <  1:2989852/53‑1 (MQ=255)
                                        gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                    >  1:2414015/1‑60 (MQ=255)
                                          aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                  >  1:996348/1‑44 (MQ=255)
                                          aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                  >  1:1467276/1‑44 (MQ=255)
                                           gATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                      >  1:3382425/1‑39 (MQ=255)
                                           gATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                        >  1:3680665/1‑37 (MQ=255)
                                           gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTc                                         >  1:3413156/1‑52 (MQ=255)
                                           gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                      >  1:1783812/1‑71 (MQ=255)
                                           gATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcg                      <  1:1160/71‑1 (MQ=255)
                                            aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                               <  1:3525389/61‑1 (MQ=255)
                                            aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                     <  1:1868592/71‑1 (MQ=255)
                                            aTGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgc                     >  1:104555/1‑71 (MQ=255)
                                             tGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                 <  1:1810735/58‑1 (MQ=255)
                                             tGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                          >  1:3398013/1‑49 (MQ=255)
                                              gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                           >  1:2509847/1‑47 (MQ=255)
                                              gAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                           >  1:79207/1‑47 (MQ=255)
                                               aGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGc                  >  1:89350/1‑71 (MQ=255)
                                                gcgcGATGTTTTTACCAGTGTGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                 <  1:2880690/71‑1 (MQ=255)
                                                gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                 <  1:2324083/71‑1 (MQ=255)
                                                gcgcGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                 <  1:2816680/71‑1 (MQ=255)
                                                 cgcgATGTTTTTAGCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                <  1:1296647/71‑1 (MQ=255)
                                                 cgcgATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGcaa                <  1:3578651/71‑3 (MQ=255)
                                                  gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTa                                                  <  1:307647/36‑1 (MQ=255)
                                                   cgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                             >  1:726269/1‑56 (MQ=255)
                                                     aTGTTTTTACCAGTGAGGTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                <  1:728554/67‑1 (MQ=255)
                                                      tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                   >  1:3055440/1‑47 (MQ=255)
                                                      tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCg                                   >  1:194072/1‑47 (MQ=255)
                                                       gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                              >  1:1636867/1‑51 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                        <  1:451192/40‑1 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                        <  1:792723/40‑1 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCt                                        <  1:2845223/40‑1 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                            <  1:1537903/52‑1 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag         <  1:2830672/71‑1 (MQ=255)
                                                        tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgca          <  1:2776959/70‑1 (MQ=255)
                                                          tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag         >  1:455139/1‑69 (MQ=255)
                                                               cAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTg  <  1:1053806/71‑1 (MQ=255)
                                                                |                                                                     
CAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTG  >  minE/2919500‑2919633

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: