Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,198 A→G 20.3% R166R (CGT→CGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1980AG20.3% 231.4 / 52.8 113R166R (CGT→CGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (40/50);  new base G (16/7);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 3.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 121 total aligned reads displayed.

TGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTG  >  minE/2919129‑2919267
                                                                     |                                                                     
tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                                      <  1:553622/71‑1 (MQ=255)
    gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                                                  <  1:2175231/71‑1 (MQ=255)
    gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                                                  <  1:2179825/71‑1 (MQ=255)
       aacAGGCGCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATc                                                               >  1:242975/3‑71 (MQ=255)
       aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGt                                                                >  1:1056823/1‑70 (MQ=255)
       aCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTc                                                               <  1:1261596/71‑1 (MQ=255)
           gACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCt                                                           >  1:3557828/1‑71 (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)
             cGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc                                                         <  1:1517398/71‑1 (MQ=255)
              gCAGATCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGcg                                                        >  1:2437121/1‑70 (MQ=255)
                agaTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTa                                                      <  1:1072696/69‑1 (MQ=255)
                aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                                        >  1:2616502/1‑69 (MQ=255)
                aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                                      <  1:3293749/71‑1 (MQ=255)
                 aaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAg                                                     <  1:3413603/71‑1 (MQ=255)
                  aTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTgcgcgca                                                                 <  1:3216249/58‑2 (MQ=25)
                  aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCGATGTTGTTACGCGCGTCACCTGCATAg                                                     <  1:1724833/70‑1 (MQ=255)
                  aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                    >  1:2904147/1‑71 (MQ=255)
                  aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                    >  1:992121/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)
                          cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                                      <  1:1661539/45‑1 (MQ=255)
                          cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                                             >  1:3549354/1‑70 (MQ=255)
                            gtcaGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCAt                                          >  1:371926/4‑71 (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: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)
                                         cTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGc                                                         <  1:3952579/43‑1 (MQ=255)
                                         cTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt                             <  1:1992691/71‑1 (MQ=255)
                                         cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt                                     <  1:361106/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:299123/70‑1 (MQ=255)
                                            cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                           <  1:839445/70‑1 (MQ=255)
                                            cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCt                             >  1:3850122/1‑68 (MQ=255)
                                             gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                          <  1:1436368/70‑1 (MQ=255)
                                             gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                          <  1:3146009/70‑1 (MQ=255)
                                              aGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                        <  1:3287714/71‑1 (MQ=255)
                                                cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                                              <  1:2374444/47‑1 (MQ=255)
                                                cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGcccg                      >  1:3841291/1‑69 (MQ=255)
                                                cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg                       <  1:694927/70‑1 (MQ=255)
                                                 aTCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                        >  1:1172768/1‑68 (MQ=255)
                                                  tCGAGTTGCCCATGTTGTTGCGCGCATAGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCctt                           >  1:1901390/1‑64 (MQ=255)
                                                    gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                  <  1:2893689/71‑1 (MQ=255)
                                                    gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa                  <  1:3745365/71‑1 (MQ=255)
                                                     aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                        <  1:1063702/64‑1 (MQ=255)
                                                     aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                 >  1:3741343/1‑71 (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:2687687/55‑1 (MQ=255)
                                                        tGCCCATGTTGTTGCACGCATCGCCCGCGTAGACCAGCGTCATATCGTTAAACGc                              >  1:1316265/1‑55 (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)
                                                             aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc              <  1:2935722/66‑1 (MQ=255)
                                                              tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                       <  1:243703/40‑1 (MQ=255)
                                                              tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc            >  1:2627365/1‑67 (MQ=255)
                                                              tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg         >  1:793078/1‑70 (MQ=255)
                                                               gttgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa       >  1:2923270/1‑71 (MQ=255)
                                                                ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa       <  1:789498/70‑1 (MQ=255)
                                                                ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa       <  1:2089736/70‑1 (MQ=255)
                                                                ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa       <  1:1330796/70‑1 (MQ=255)
                                                                ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCAt      >  1:1159080/1‑71 (MQ=255)
                                                                ttgttACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAAt      >  1:1735195/1‑71 (MQ=255)
                                                                 tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt      <  1:465540/70‑1 (MQ=255)
                                                                  gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATgg    >  1:2865924/1‑71 (MQ=255)
                                                                  gttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg     <  1:2745274/70‑1 (MQ=255)
                                                                   ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg    <  1:1224030/70‑1 (MQ=255)
                                                                   ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg    <  1:358280/70‑1 (MQ=255)
                                                                   ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt   <  1:1745529/71‑1 (MQ=255)
                                                                     aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt                                    >  1:861328/1‑36 (MQ=255)
                                                                     aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt   >  1:291980/1‑69 (MQ=255)
                                                                     aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTg  <  1:3118263/70‑1 (MQ=255)
                                                                     |                                                                     
TGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTG  >  minE/2919129‑2919267

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: