Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,183 A→G 15.6% N171N (AAT→AAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,1830AG15.6% 227.2 / 31.7 96N171N (AAT→AACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (36/45);  new base G (10/5);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 104 total aligned reads displayed.

CCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAA  >  minE/2919113‑2919252
                                                                      |                                                                     
ccGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGaa                                                                       >  1:3638961/1‑71 (MQ=255)
 cGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAAt                                                                      >  1:1093903/1‑71 (MQ=255)
     cAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                                  <  1:258419/71‑1 (MQ=255)
      aGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGcc                                                                  <  1:3182796/70‑1 (MQ=255)
         aCGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtg                                                              >  1:3831688/1‑71 (MQ=255)
           gCTTGTGGCGCGACCAGACGCAAATCCGGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                                                            <  1:750023/71‑1 (MQ=255)
           gCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                                                            <  1:2695455/71‑1 (MQ=255)
            cTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgtt                                                            <  1:225363/70‑1 (MQ=255)
               gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                          >  1:1885405/1‑69 (MQ=255)
               gtgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttgt                                                          <  1:2005970/69‑1 (MQ=255)
                tgGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                       <  1:553622/71‑1 (MQ=255)
                  gcgcGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGGATCGAATTGCCCATGTTGTTAcgc                                                     >  1:3063692/1‑71 (MQ=255)
                   cgcgTCCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAtgttg                                                           <  1:1486235/64‑1 (MQ=255)
                    gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCAt                                                               <  1:2133395/59‑1 (MQ=255)
                    gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                                   <  1:2179825/71‑1 (MQ=255)
                    gcgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAcgcgc                                                   <  1:2175231/71‑1 (MQ=255)
                     cgACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATgttgtt                                                         >  1:783572/1‑64 (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)
                               cAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCAt                                        <  1:2048631/71‑1 (MQ=255)
                                agaTCCAGCCCGGTCAGCGCCGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTa                                       <  1:1072696/69‑1 (MQ=255)
                                aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                         >  1:2616502/1‑69 (MQ=255)
                                aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                       <  1:3293749/71‑1 (MQ=255)
                                aaaTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATa                                       <  1:1711077/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)
                                     cAGCCCGGTCAGCGCGGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAc                                    >  1:2948535/1‑69 (MQ=255)
                                     cAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCa                                               >  1:1661717/1‑58 (MQ=255)
                                       gACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGc                                >  1:1582406/1‑71 (MQ=255)
                                        acaggTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCg                               >  1:18140/4‑71 (MQ=255)
                                        aCCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                                  <  1:383810/68‑1 (MQ=255)
                                         ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                              <  1:205712/71‑1 (MQ=255)
                                         ccGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGt                              >  1:3004126/1‑71 (MQ=255)
                                          cGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTAc                                                       <  1:1661539/45‑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:740140/1‑71 (MQ=255)
                                               aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                        <  1:910884/71‑1 (MQ=255)
                                               aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                        >  1:3935842/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:365812/63‑1 (MQ=255)
                                                         cTTCGAGCAGCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGt                      <  1:361106/63‑1 (MQ=255)
                                                          ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt             >  1:1581253/1‑71 (MQ=255)
                                                          ttCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCtt             >  1:1117797/1‑71 (MQ=255)
                                                           tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCa                                         >  1:2467971/1‑42 (MQ=255)
                                                           tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                    >  1:2810286/1‑63 (MQ=255)
                                                           tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                    >  1:1166324/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:3835298/1‑68 (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)
                                                                cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                               <  1:2930009/47‑1 (MQ=255)
                                                                cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGcccg       >  1:3841291/1‑69 (MQ=255)
                                                                cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg       <  1:886526/71‑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:762208/71‑1 (MQ=255)
                                                                    gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaa   <  1:3745365/71‑1 (MQ=255)
                                                                    gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa    <  1:58346/70‑1 (MQ=255)
                                                                    gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa    <  1:2583324/70‑1 (MQ=255)
                                                                     aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc         <  1:1063702/64‑1 (MQ=255)
                                                                     aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa  <  1:3235175/71‑1 (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)
                                                                      |                                                                     
CCGGCCAGCACGCTTGTGGCGCGACCAGACGCAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAA  >  minE/2919113‑2919252

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: