Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,252 A→G 24.5% H148H (CAT→CAC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2520AG24.5% 165.9 / 54.5 98H148H (CAT→CACargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (32/42);  new base G (12/12);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 103 total aligned reads displayed.

ATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAT  >  minE/2919183‑2919321
                                                                     |                                                                     
aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAAt                                                                      >  1:1202170/1‑71 (MQ=255)
aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAAt                                                                      >  1:3432057/1‑71 (MQ=255)
  tGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                                    >  1:663767/1‑71 (MQ=255)
   gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                                    <  1:3273038/70‑1 (MQ=255)
   gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                                    <  1:181489/70‑1 (MQ=255)
   gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                                                    <  1:3308099/70‑1 (MQ=255)
   gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCt                                                                   >  1:5246231/1‑71 (MQ=255)
    cccATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  <  1:5137429/71‑1 (MQ=255)
      cATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCt                                                                <  1:1725906/71‑1 (MQ=255)
        tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  <  1:5407096/67‑1 (MQ=255)
        tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  <  1:4580676/67‑1 (MQ=255)
        tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGc                                                              <  1:2808378/71‑1 (MQ=255)
          ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                                       >  1:3554266/1‑60 (MQ=255)
          ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  >  1:1514884/1‑65 (MQ=255)
          ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                                  >  1:4593094/1‑65 (MQ=255)
          ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                             <  1:5141844/70‑1 (MQ=255)
           tgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                                           >  1:911300/1‑71 (MQ=255)
           tgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                                           >  1:3016110/1‑71 (MQ=255)
           tgttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                                           >  1:2632796/1‑71 (MQ=255)
            gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCt                                                                   >  1:647392/1‑62 (MQ=255)
             ttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATTGt                                                         <  1:2013748/71‑1 (MQ=255)
              tACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                                           >  1:4202588/1‑68 (MQ=255)
              tACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                                           >  1:2044921/1‑68 (MQ=255)
              tACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                                                         >  1:3598803/1‑70 (MQ=255)
               gcgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtca                                                       >  1:650508/2‑69 (MQ=255)
                cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtca                                                       <  1:846672/70‑3 (MQ=255)
                  cgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGc                                                              <  1:4367059/61‑1 (MQ=255)
                     gTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                      >  1:5313960/1‑66 (MQ=255)
                     gTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                 <  1:4850782/71‑1 (MQ=255)
                       cACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                                               <  1:3114670/55‑1 (MQ=255)
                       cACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                               >  1:2351368/1‑71 (MQ=255)
                        aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                                              >  1:1206329/1‑71 (MQ=255)
                           tGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGc                                             >  1:1455516/1‑69 (MQ=255)
                            gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGcc                                            <  1:5354942/69‑1 (MQ=255)
                            gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGcc                                            <  1:5039536/69‑1 (MQ=255)
                             cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCt                                                                   <  1:4560080/45‑1 (MQ=255)
                             cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCt                                                                   <  1:3144619/45‑1 (MQ=255)
                              gtAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                        <  1:3256820/70‑1 (MQ=255)
                              aTAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagca                                        <  1:2980714/71‑1 (MQ=255)
                               tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                       <  1:200232/71‑1 (MQ=255)
                                 gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                             >  1:914893/1‑47 (MQ=255)
                                 gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     <  1:5519083/71‑1 (MQ=255)
                                 gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     <  1:2168190/71‑1 (MQ=255)
                                 gACCAGCGCCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                     >  1:5551481/1‑71 (MQ=255)
                                  aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                    <  1:3663917/71‑1 (MQ=255)
                                  aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                    >  1:4465203/1‑71 (MQ=255)
                                  aCCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTTCATGGTCATCAGGTCCGCCAGCAGCt                                     <  1:79489/70‑1 (MQ=255)
                                   ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                       >  1:3137687/1‑67 (MQ=255)
                                   ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                   <  1:3548393/71‑1 (MQ=255)
                                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                     >  1:3230000/1‑68 (MQ=255)
                                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                   <  1:2206133/70‑1 (MQ=255)
                                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                  >  1:3087521/1‑71 (MQ=255)
                                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                  <  1:2776115/71‑1 (MQ=255)
                                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTACTGCATGGTGAGAAGATCCGCcagcag                                       >  1:5283254/1‑66 (MQ=255)
                                      gCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGtc                                <  1:3738147/71‑2 (MQ=255)
                                       cGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGtc                                <  1:1594569/70‑2 (MQ=255)
                                        gTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGt                                                         >  1:3611139/1‑44 (MQ=16)
                                        gTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTcgg                              <  1:2875392/71‑1 (MQ=255)
                                         tCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTcggg                             <  1:988855/71‑2 (MQ=255)
                                         tCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTcggg                             >  1:958237/1‑70 (MQ=255)
                                           aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                 >  1:4524402/1‑49 (MQ=255)
                                           aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                 >  1:5114721/1‑49 (MQ=255)
                                           aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                            >  1:1165251/1‑70 (MQ=255)
                                           aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                            >  1:2937142/1‑70 (MQ=255)
                                           aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                            <  1:3775002/70‑1 (MQ=255)
                                            tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                      <  1:1933938/43‑1 (MQ=255)
                                            tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                      <  1:3301538/43‑1 (MQ=255)
                                            tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                           >  1:1521858/1‑70 (MQ=255)
                                             ttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:98710/70‑1 (MQ=255)
                                             ttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          >  1:4612318/1‑70 (MQ=255)
                                             ttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          >  1:2694170/1‑70 (MQ=255)
                                              tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTcggg                             >  1:1377731/1‑65 (MQ=255)
                                              tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTcggg                             >  1:203654/1‑65 (MQ=255)
                                              tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                        >  1:4193708/1‑71 (MQ=255)
                                                gTTGTACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                      <  1:5298801/71‑1 (MQ=255)
                                                gTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                      <  1:2503699/71‑1 (MQ=255)
                                                gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTgg                          <  1:2222422/67‑1 (MQ=255)
                                                 ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                      <  1:4187329/70‑1 (MQ=255)
                                                  tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagc                                         <  1:2175456/50‑1 (MQ=255)
                                                  tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                          <  1:5487741/65‑1 (MQ=255)
                                                  tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTcg                    >  1:3286920/1‑70 (MQ=255)
                                                  tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                    >  1:94922/1‑71 (MQ=255)
                                                   gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   >  1:641879/1‑71 (MQ=255)
                                                   gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   >  1:3541285/1‑71 (MQ=255)
                                                   aaaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTgg                          <  1:2672590/63‑1 (MQ=255)
                                                   aaaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcg                    <  1:1719064/69‑2 (MQ=255)
                                                    aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                   <  1:893408/70‑3 (MQ=255)
                                                      cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                >  1:1824122/1‑71 (MQ=255)
                                                       gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt               <  1:3591007/71‑1 (MQ=255)
                                                            gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCTAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg           <  1:3382133/70‑1 (MQ=255)
                                                            gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg           <  1:5174129/70‑1 (MQ=255)
                                                             ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                      <  1:4193588/58‑1 (MQ=255)
                                                             ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc          <  1:3721670/70‑1 (MQ=255)
                                                             ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc          <  1:4659488/70‑1 (MQ=255)
                                                                 gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc     <  1:71432/71‑1 (MQ=255)
                                                                 gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc     <  1:1168233/71‑1 (MQ=255)
                                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   >  1:724246/1‑54 (MQ=255)
                                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   >  1:682287/1‑54 (MQ=255)
                                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                   >  1:5669422/1‑54 (MQ=255)
                                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt  <  1:2454188/71‑1 (MQ=255)
                                                                     |                                                                     
ATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAT  >  minE/2919183‑2919321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: