Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,567 G→A 12.8% 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,5670GA12.8% 298.4 / 27.7 117L43L (CTC→CTTargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (44/58);  new base A (6/9);  total (50/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 119 total aligned reads displayed.

GGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTG  >  minE/2919498‑2919636
                                                                     |                                                                     
ggCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGa                                                                      <  1:3349754/71‑1 (MQ=255)
 gCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                                     <  1:5144288/71‑1 (MQ=255)
  cAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAgc                                                                    <  1:1253813/71‑2 (MQ=255)
    aCTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                                  <  1:3504061/71‑1 (MQ=255)
     cTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg                                                                 <  1:2213201/71‑1 (MQ=255)
     cTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTg                                                                 <  1:701491/71‑1 (MQ=255)
       tCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttctg                                                               >  1:1113318/1‑68 (MQ=255)
         gAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGctt                                                             <  1:176865/71‑1 (MQ=255)
          aaaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   <  1:3507232/64‑1 (MQ=255)
          aaaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   <  1:1300597/64‑1 (MQ=255)
          aaaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   <  1:5007830/64‑1 (MQ=255)
            agagCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                           >  1:1687932/1‑70 (MQ=255)
              agCATCTGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                         <  1:5378828/70‑1 (MQ=255)
               gCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttc                                                            >  1:1738901/1‑66 (MQ=255)
                  tCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGt                                                                    <  1:3143502/55‑1 (MQ=255)
                  tCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGt                                                                    <  1:5598032/55‑1 (MQ=255)
                   cGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGttt                                                                  <  1:4079618/56‑1 (MQ=255)
                   cGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   <  1:1294892/55‑1 (MQ=255)
                    gggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttct                                                           >  1:163657/1‑62 (MQ=255)
                    gggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCg                                                  >  1:3633395/1‑71 (MQ=255)
                     ggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGc                                                 >  1:2542096/1‑71 (MQ=255)
                         cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                             <  1:2151520/71‑1 (MQ=255)
                            gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          <  1:4606374/71‑1 (MQ=255)
                            gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          <  1:314772/71‑1 (MQ=255)
                            gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          <  1:4875748/71‑1 (MQ=255)
                             tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                                     <  1:1699371/43‑1 (MQ=255)
                             tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAg                                                                     <  1:2563633/43‑1 (MQ=255)
                             tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttc                                                         <  1:4788046/55‑1 (MQ=255)
                               gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttcat                                                       >  1:1111453/1‑53 (MQ=255)
                                agtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       <  1:4496425/70‑1 (MQ=255)
                                 gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTgg                                                                >  1:201510/1‑44 (MQ=255)
                                 gtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTgg                                                                >  1:42092/1‑44 (MQ=255)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttc                                                                 <  1:3112486/42‑2 (MQ=25)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGcttc                                                                 <  1:395573/42‑2 (MQ=25)
                                  tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGc                                    <  1:892347/71‑1 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTTCCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                              <  1:3720964/59‑1 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtt                                  >  1:1066549/1‑71 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGct                                                              >  1:4191883/1‑43 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTAcc                                                   <  1:4653121/54‑1 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                              <  1:3109575/59‑1 (MQ=255)
                                    tttttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCttt                                              <  1:836030/59‑1 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   >  1:251679/1‑36 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGc                                                               >  1:1683283/1‑40 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGc                                                 >  1:4436992/1‑54 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       >  1:2838773/1‑64 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       >  1:3489964/1‑64 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                 <  1:3048147/70‑1 (MQ=255)
                                      tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:2952393/69‑1 (MQ=255)
                                       ttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttg                                <  1:4069147/70‑2 (MQ=255)
                                       ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                       <  1:4276065/47‑1 (MQ=255)
                                       ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                 >  1:5548960/1‑69 (MQ=255)
                                       ttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTc                                <  1:4869782/70‑1 (MQ=255)
                                        tCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAg                              >  1:67537/1‑71 (MQ=255)
                                        tCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                       >  1:3784716/1‑46 (MQ=255)
                                          gAAGATGAGCGGGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCt                            <  1:3950715/71‑1 (MQ=255)
                                          gAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCt                            <  1:1862618/71‑1 (MQ=255)
                                          gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                               <  1:1361473/68‑1 (MQ=255)
                                          gAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCt                            <  1:3297626/71‑1 (MQ=255)
                                           aaGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                           <  1:5082994/71‑1 (MQ=255)
                                           aaGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCATTCAGCtt                           <  1:4292735/71‑1 (MQ=255)
                                            aGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCtgt                          >  1:3697866/1‑69 (MQ=255)
                                            aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCtt                           <  1:3693964/70‑1 (MQ=255)
                                            aGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTc                          >  1:3863955/1‑71 (MQ=255)
                                              aTGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttg                                <  1:3547569/63‑2 (MQ=255)
                                               tGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgcg                       >  1:1786043/1‑71 (MQ=255)
                                               tGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTcgcg                       >  1:3365949/1‑71 (MQ=255)
                                                 tgcgcGATTTTTTTACCAGTGAGTTTGGCTTCTTCttt                                                      <  1:5185375/37‑1 (MQ=38)
                                                 aGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                             >  1:4260061/1‑47 (MQ=255)
                                                 aGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTACTTTACCGCtttt                                             >  1:652379/1‑47 (MQ=255)
                                                  gcgcGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCt                                   >  1:1133370/1‑56 (MQ=255)
                                                   cgcgATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCtttga                               >  1:2565545/1‑57 (MQ=255)
                                                   cgcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTAGCGGCTa                   <  1:1875674/71‑1 (MQ=255)
                                                    gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTaa                  <  1:1695396/71‑1 (MQ=255)
                                                    gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTaa                  <  1:5586578/71‑1 (MQ=255)
                                                    gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTa                   <  1:615408/70‑1 (MQ=255)
                                                        tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCg                                                  >  1:5546672/1‑35 (MQ=255)
                                                        tGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgc              <  1:915484/71‑1 (MQ=255)
                                                         gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGc                                                 >  1:3147086/1‑35 (MQ=255)
                                                         gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCa                               >  1:913339/1‑53 (MQ=255)
                                                         gTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACt                >  1:3403086/1‑68 (MQ=255)
                                                          tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                             <  1:486361/38‑1 (MQ=255)
                                                          tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCtttt                                             <  1:5010164/38‑1 (MQ=255)
                                                           ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       >  1:3951347/1‑43 (MQ=255)
                                                           ttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       >  1:514493/1‑43 (MQ=255)
                                                           ttttACCAGTGAGATTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttta                                >  1:157661/1‑49 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          >  1:2322836/1‑39 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          >  1:266719/1‑39 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          >  1:2908816/1‑39 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCttt                                 >  1:1046939/1‑48 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGc                             >  1:1053741/1‑52 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcagc           <  1:532089/70‑1 (MQ=255)
                                                            tttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgcag            >  1:4765830/1‑69 (MQ=255)
                                                                   gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATc                                       >  1:3231310/1‑35 (MQ=255)
                                                                   gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGtt   >  1:3736610/1‑71 (MQ=255)
                                                                    tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGcagcag         <  1:167875/64‑1 (MQ=255)
                                                                    tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGcagcag         <  1:1345453/64‑1 (MQ=255)
                                                                    tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:287974/71‑1 (MQ=255)
                                                                    tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:589972/71‑1 (MQ=255)
                                                                    tGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:2613338/71‑1 (MQ=255)
                                                                     gAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:451904/70‑1 (MQ=255)
                                                                     |                                                                     
GGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTG  >  minE/2919498‑2919636

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: