Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,213 A→G 23.0% Y161Y (TAT→TAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2130AG23.0% 209.8 / 63.7 113Y161Y (TAT→TACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (39/48);  new base G (16/10);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 115 total aligned reads displayed.

CAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAG  >  minE/2919144‑2919281
                                                                     |                                                                    
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)
   aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCGATGTTGTTACGCGCGTCACCTGCATAg                                                                   <  1:1724833/70‑1 (MQ=255)
   aTCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGa                                                                  >  1:992121/1‑71 (MQ=255)
      cAGCCCGGTCAGCGCGGCCGCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGAc                                                                 >  1:2948535/1‑69 (MQ=255)
        gACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGc                                                             >  1:1582406/1‑71 (MQ=255)
         acaggTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCg                                                            >  1:18140/4‑71 (MQ=255)
         aCCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCa                                                               <  1:383810/68‑1 (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)
                      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: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)
                             cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                         <  1:839445/70‑1 (MQ=255)
                             cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                         <  1:299123/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)
                                 cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                                                            <  1:2930009/47‑1 (MQ=255)
                                 cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCg                                                            <  1:2374444/47‑1 (MQ=255)
                                 cATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGcccg                                    >  1:3841291/1‑69 (MQ=255)
                                 cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                    <  1:886526/71‑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:2583324/70‑1 (MQ=255)
                                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                 <  1:58346/70‑1 (MQ=255)
                                      aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                                      <  1:1063702/64‑1 (MQ=255)
                                      aaTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                               >  1:3741343/1‑71 (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:3546383/55‑1 (MQ=255)
                                       aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAAc                                              <  1:2687687/55‑1 (MQ=255)
                                       aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAAc                                              <  1:668902/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)
                                             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:2860423/1‑70 (MQ=255)
                                               tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                       >  1:793078/1‑70 (MQ=255)
                                                gttgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                     >  1:2923270/1‑71 (MQ=255)
                                                 ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                     <  1:1330796/70‑1 (MQ=255)
                                                 ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                     <  1:789498/70‑1 (MQ=255)
                                                 ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCa                     <  1:2089736/70‑1 (MQ=255)
                                                 ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCAt                    >  1:1159080/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)
                                                   gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATgg                  >  1:2865924/1‑71 (MQ=255)
                                                   gttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                   <  1:2745274/70‑1 (MQ=255)
                                                    ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                  <  1:358280/70‑1 (MQ=255)
                                                    ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                  <  1:1224030/70‑1 (MQ=255)
                                                    ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                 <  1:1745529/71‑1 (MQ=255)
                                                      gcgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                 >  1:3906515/2‑53 (MQ=255)
                                                      gcgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                 >  1:428561/2‑53 (MQ=255)
                                                      aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                 >  1:291980/1‑69 (MQ=255)
                                                      aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTg                <  1:3118263/70‑1 (MQ=255)
                                                       cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtca               <  1:414890/70‑3 (MQ=255)
                                                        gcgcgTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag              >  1:3389424/1‑70 (MQ=255)
                                                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                 <  1:2678815/49‑1 (MQ=255)
                                                             tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc        <  1:984147/71‑1 (MQ=255)
                                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCtt                                          >  1:331443/1‑35 (MQ=255)
                                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCtt                                          >  1:3142143/1‑35 (MQ=255)
                                                                 cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga          <  1:1332913/65‑1 (MQ=255)
                                                                   gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATgg                  <  1:1371044/55‑1 (MQ=255)
                                                                   gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag  >  1:2729199/1‑71 (MQ=255)
                                                                   gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag  >  1:3375487/1‑71 (MQ=255)
                                                                    cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag  >  1:215778/1‑70 (MQ=255)
                                                                     |                                                                    
CAAATCCAGACCGGTAAGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAG  >  minE/2919144‑2919281

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: