Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,228 T→C 17.8% E156E (GAA→GAG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2280TC17.8% 191.7 / 35.8 84E156E (GAA→GAGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (39/30);  new base C (8/7);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

AGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGG  >  minE/2919160‑2919297
                                                                    |                                                                     
aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                                                     >  1:73158/1‑71 (MQ=255)
aGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTc                                                                     >  1:2368877/1‑71 (MQ=255)
    ccgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGtt                                                                  >  1:929025/1‑70 (MQ=255)
     cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                                                                <  1:2055691/71‑1 (MQ=255)
     cgcAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGa                                                                <  1:1296743/71‑1 (MQ=255)
         gCTTCCAGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGc                                                            >  1:866859/1‑71 (MQ=255)
            tCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                                         >  1:2638871/1‑71 (MQ=255)
             ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg                                                        >  1:2298655/3‑71 (MQ=255)
             ccaGCATCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTg                                                        >  1:573383/3‑71 (MQ=255)
             cGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCttt                                                         <  1:3383346/70‑1 (MQ=255)
              gAGCATCGAATTGCCCATGTTGTTAGGCGTGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                       <  1:1053910/71‑1 (MQ=255)
              gAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                       >  1:3029943/1‑71 (MQ=255)
               aGCATCGAATTGCCCATGTTGTTAAGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGcc                                                      <  1:2289205/71‑1 (MQ=255)
                gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg                                                     <  1:902106/71‑1 (MQ=255)
                 cATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGc                                                       <  1:3393915/68‑1 (MQ=255)
                   tCGAGTTGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGttaa                                                                >  1:3008743/1‑55 (MQ=255)
                   tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                                                 <  1:3627368/56‑1 (MQ=255)
                   tCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTg                                                                 <  1:1649990/56‑1 (MQ=255)
                     gagttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                                 <  1:1098310/67‑1 (MQ=255)
                     gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                <  1:4103696/71‑1 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGaa                                                               >  1:547197/1‑56 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 >  1:481889/1‑70 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 <  1:598495/70‑1 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 <  1:2921634/70‑1 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 <  1:132593/70‑1 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 >  1:2397820/1‑70 (MQ=255)
                     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                 >  1:93207/1‑70 (MQ=255)
                       gttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGtta                                                                 <  1:3533291/51‑2 (MQ=255)
                       aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGc                                                  >  1:1354108/1‑67 (MQ=255)
                        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                             >  1:110867/1‑71 (MQ=255)
                        ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                             >  1:1127789/1‑71 (MQ=255)
                          gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                    >  1:635454/1‑62 (MQ=255)
                          gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                    >  1:2797244/1‑62 (MQ=255)
                              aTGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg                                       <  1:2449036/71‑1 (MQ=255)
                              aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGc                                            <  1:3455239/66‑1 (MQ=255)
                              aTGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                       <  1:2386096/71‑1 (MQ=255)
                               tgttgtTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTg                                       >  1:2568003/1‑70 (MQ=255)
                               tgttgtTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                       >  1:1100509/1‑70 (MQ=255)
                                gttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCt                                        >  1:398805/1‑68 (MQ=255)
                                 ttgttgCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGttaa                                                                >  1:3163357/1‑41 (MQ=16)
                                 ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                          >  1:1663463/1‑65 (MQ=255)
                                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                     >  1:364041/1‑69 (MQ=255)
                                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                     >  1:359370/1‑69 (MQ=255)
                                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATg                                   >  1:2263890/1‑71 (MQ=255)
                                   gttGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATgg                                  >  1:94299/1‑71 (MQ=255)
                                     tACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGt                                 >  1:2504116/1‑70 (MQ=255)
                                      aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                     >  1:2431715/1‑65 (MQ=255)
                                       cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                              <  1:1368559/71‑1 (MQ=255)
                                       cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                              <  1:1267869/71‑1 (MQ=255)
                                       cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtcat                              <  1:649083/71‑4 (MQ=255)
                                       cgcgcaTCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGtcat                              <  1:2559838/71‑4 (MQ=255)
                                        gcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                       <  1:2697759/61‑1 (MQ=255)
                                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                                                        <  1:915770/42‑1 (MQ=255)
                                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                                                        <  1:3266232/42‑1 (MQ=255)
                                          gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                              >  1:41487/1‑68 (MQ=255)
                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                               >  1:3445416/1‑46 (MQ=255)
                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                       <  1:3590668/70‑1 (MQ=255)
                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                       <  1:331130/70‑1 (MQ=255)
                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                      >  1:789079/1‑71 (MQ=255)
                                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                      >  1:2325574/1‑71 (MQ=255)
                                                 cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGcc                    <  1:1040683/71‑1 (MQ=255)
                                                   gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                             >  1:3395957/1‑44 (MQ=255)
                                                   gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                  >  1:1286198/1‑71 (MQ=255)
                                                    cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                  >  1:3527105/1‑70 (MQ=255)
                                                    cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                  >  1:4013411/1‑70 (MQ=255)
                                                     gtAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                <  1:2285979/70‑1 (MQ=255)
                                                       aGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc              >  1:2015197/1‑71 (MQ=255)
                                                         aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                             >  1:2276215/1‑38 (MQ=255)
                                                         aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg            >  1:2004543/1‑71 (MQ=255)
                                                         aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg            <  1:2004264/71‑1 (MQ=255)
                                                           cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg          <  1:85451/71‑1 (MQ=255)
                                                           cAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATg                                   <  1:3059441/46‑1 (MQ=255)
                                                             gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCAt                                    <  1:2488407/43‑1 (MQ=255)
                                                             gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg        >  1:1666852/1‑71 (MQ=255)
                                                             gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg        <  1:1871423/71‑1 (MQ=255)
                                                             gCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGt                                 >  1:1190934/1‑46 (MQ=255)
                                                              cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg       >  1:789520/1‑71 (MQ=255)
                                                              cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg       <  1:2751769/71‑1 (MQ=255)
                                                               gTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg       >  1:3600264/1‑70 (MQ=255)
                                                                tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTg                                       <  1:680515/37‑1 (MQ=255)
                                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCGTGCa                                     >  1:1390525/1‑37 (MQ=255)
                                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg   >  1:990945/1‑71 (MQ=255)
                                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg   <  1:3168188/71‑1 (MQ=255)
                                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTCGGATg   >  1:2055453/1‑71 (MQ=255)
                                                                   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                          <  1:497683/47‑1 (MQ=255)
                                                                   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaag                           <  1:307900/46‑1 (MQ=255)
                                                                   tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg  <  1:2113160/71‑1 (MQ=255)
                                                                    |                                                                     
AGCGCCGCAGCTTCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGG  >  minE/2919160‑2919297

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: