Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,204 G→A 23.5% D164D (GAC→GAT argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2040GA23.5% 218.0 / 64.0 115D164D (GAC→GATargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (41/47);  new base A (17/10);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 121 total aligned reads displayed.

ACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAG  >  minE/2919136‑2919269
                                                                    |                                                                 
aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                                                                  >  1:1056823/1‑70 (MQ=255)
aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTc                                                                 <  1:1261596/71‑1 (MQ=255)
    gACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCt                                                             >  1:282821/1‑71 (MQ=255)
     gcGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGccc                                                             >  1:359579/2‑69 (MQ=255)
     ccGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                                            >  1:652588/2‑71 (MQ=255)
     aCGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTg                                                            >  1:281883/1‑71 (MQ=255)
      cGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGccc                                                             <  1:1487388/69‑2 (MQ=255)
       gCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGcg                                                          >  1:2437121/1‑70 (MQ=255)
         agaTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTa                                                        <  1:1072696/69‑1 (MQ=255)
         aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                                          >  1:2616502/1‑69 (MQ=255)
          aaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAg                                                       <  1:3413603/71‑1 (MQ=255)
           aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCGATGTTGTTACGCGCGTCACCTGCATAg                                                       <  1:1724833/70‑1 (MQ=255)
           aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                      >  1:992121/1‑71 (MQ=255)
           aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                      >  1:2904147/1‑71 (MQ=255)
              cAGCCCGGTCAGCGCGGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAc                                                     >  1:2948535/1‑69 (MQ=255)
              cAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                                                >  1:1661717/1‑58 (MQ=255)
                gACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGc                                                 >  1:1582406/1‑71 (MQ=255)
                 aCCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                                                   <  1:383810/68‑1 (MQ=255)
                  ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                                               >  1:3004126/1‑71 (MQ=255)
                  ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                                               <  1:205712/71‑1 (MQ=255)
                   cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                                               >  1:3549354/1‑70 (MQ=255)
                     gtcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAt                                            >  1:371926/4‑71 (MQ=255)
                       aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGAc                                                     >  1:2726544/1‑60 (MQ=255)
                       aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTc                                              >  1:3090568/1‑67 (MQ=255)
                       aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAttt                                          >  1:721565/1‑71 (MQ=255)
                       aaGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAtt                                           >  1:3856891/1‑70 (MQ=255)
                        aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                         >  1:3935842/1‑71 (MQ=255)
                        aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                         <  1:910884/71‑1 (MQ=255)
                        aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                         >  1:740140/1‑71 (MQ=255)
                        aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGAATAGACCAGCGTCATTTc                                         >  1:2800822/1‑71 (MQ=255)
                          cgccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGGCACCTGCATAGACCAGCGTCATTTCGt                                       <  1:1808402/71‑1 (MQ=255)
                           gccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt                                      <  1:199087/71‑1 (MQ=255)
                            ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                     <  1:1578336/71‑1 (MQ=255)
                            ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                     <  1:478987/71‑1 (MQ=255)
                             cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTAGGCGCGTCACCt                                                             <  1:1501341/46‑1 (MQ=255)
                              gcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGttaac                                   >  1:3249707/1‑68 (MQ=255)
                                aGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACg                                 >  1:910126/1‑71 (MQ=255)
                                 gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc                                >  1:2102341/1‑71 (MQ=255)
                                  cTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt                               <  1:1992691/71‑1 (MQ=255)
                                  cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt                                       <  1:361106/63‑1 (MQ=255)
                                  cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt                                       <  1:365812/63‑1 (MQ=255)
                                   ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt                              >  1:1117797/1‑71 (MQ=255)
                                   ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt                              >  1:1581253/1‑71 (MQ=255)
                                    tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                                                >  1:195632/1‑36 (MQ=255)
                                    tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                                          >  1:2467971/1‑42 (MQ=255)
                                    tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                     >  1:2810286/1‑63 (MQ=255)
                                     ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg                            >  1:2974639/3‑71 (MQ=255)
                                     cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                                                <  1:1986271/35‑1 (MQ=255)
                                     cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                             <  1:839445/70‑1 (MQ=255)
                                     cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                             <  1:299123/70‑1 (MQ=255)
                                     cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt                               >  1:3850122/1‑68 (MQ=255)
                                     cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt                               >  1:3835298/1‑68 (MQ=255)
                                      gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                            <  1:3146009/70‑1 (MQ=255)
                                      gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                            <  1:1436368/70‑1 (MQ=255)
                                       aGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                          <  1:3287714/71‑1 (MQ=255)
                                         cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                                                <  1:2930009/47‑1 (MQ=255)
                                         cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                                                <  1:2374444/47‑1 (MQ=255)
                                         cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg                         <  1:694927/70‑1 (MQ=255)
                                          aTCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                          >  1:1172768/1‑68 (MQ=255)
                                           tCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTc                                              >  1:1320609/1‑47 (MQ=255)
                                           tCGAGTTGCCCATGTTGTTGCGCGCATAGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCctt                             >  1:1901390/1‑64 (MQ=255)
                                             gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                    <  1:2893689/71‑1 (MQ=255)
                                             gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                    <  1:3745365/71‑1 (MQ=255)
                                             gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                     <  1:2583324/70‑1 (MQ=255)
                                              aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                          <  1:1063702/64‑1 (MQ=255)
                                              aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                   <  1:3235175/71‑1 (MQ=255)
                                               gttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTaaa                                   >  1:82953/2‑54 (MQ=255)
                                               aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCAtt                                           <  1:890348/46‑1 (MQ=255)
                                               aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAAc                                  <  1:668902/55‑1 (MQ=255)
                                               aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAAc                                  <  1:3546383/55‑1 (MQ=255)
                                               aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAAc                                  <  1:2687687/55‑1 (MQ=255)
                                                  gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                             >  1:2686168/1‑57 (MQ=255)
                                                  gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCt               >  1:554028/1‑71 (MQ=255)
                                                   cccATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGc                      <  1:3755735/63‑1 (MQ=255)
                                                    ccATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTcc             <  1:3782032/71‑1 (MQ=255)
                                                     cATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCt            >  1:3347436/1‑71 (MQ=255)
                                                      aTGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg           <  1:1265581/71‑1 (MQ=255)
                                                       tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                         <  1:243703/40‑1 (MQ=255)
                                                       tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg           >  1:793078/1‑70 (MQ=255)
                                                       tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg           >  1:2860423/1‑70 (MQ=255)
                                                        gttgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa         >  1:2923270/1‑71 (MQ=255)
                                                         ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa         <  1:2089736/70‑1 (MQ=255)
                                                         ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa         <  1:1330796/70‑1 (MQ=255)
                                                         ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa         <  1:789498/70‑1 (MQ=255)
                                                         ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCAt        >  1:1159080/1‑71 (MQ=255)
                                                         ttgttACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAAt        >  1:1735195/1‑71 (MQ=255)
                                                          tgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg       >  1:2304940/1‑71 (MQ=255)
                                                          tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt        <  1:465540/70‑1 (MQ=255)
                                                          tgttACGCGCGTCACCTGCATAGACAAGCGTCAttt                                          >  1:1518977/1‑36 (MQ=255)
                                                            ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg      <  1:1224030/70‑1 (MQ=255)
                                                            ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg      <  1:358280/70‑1 (MQ=255)
                                                            ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt     <  1:1745529/71‑1 (MQ=255)
                                                              gcgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                     >  1:428561/2‑53 (MQ=255)
                                                              gcgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                     >  1:3906515/2‑53 (MQ=255)
                                                              aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt                                      >  1:861328/1‑36 (MQ=255)
                                                              aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt     >  1:291980/1‑69 (MQ=255)
                                                               cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtca   <  1:414890/70‑3 (MQ=255)
                                                                gcgcgTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag  >  1:3389424/1‑70 (MQ=255)
                                                                  gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                     <  1:2678815/49‑1 (MQ=255)
                                                                  gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                     <  1:124757/49‑1 (MQ=255)
                                                                    |                                                                 
ACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAG  >  minE/2919136‑2919269

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: