Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,240 T→C 12.9% K152K (AAA→AAG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2400TC12.9% 248.8 / 24.7 93K152K (AAA→AAGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (47/34);  new base C (6/6);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGT  >  minE/2919172‑2919308
                                                                    |                                                                    
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)
         gagttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                                            <  1:1098310/67‑1 (MQ=255)
         gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                           <  1:4103696/71‑1 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            <  1:2921634/70‑1 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            <  1:132593/70‑1 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            >  1:481889/1‑70 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            <  1:598495/70‑1 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGAGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            >  1:2397820/1‑70 (MQ=255)
         gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                            >  1:93207/1‑70 (MQ=255)
           aTTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGc                                                             >  1:1354108/1‑67 (MQ=255)
            ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                        >  1:1127789/1‑71 (MQ=255)
            ttGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATg                                                        >  1:110867/1‑71 (MQ=255)
              gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                               >  1:2797244/1‑62 (MQ=255)
              gCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCgg                                                               >  1:635454/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)
                     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:1267869/71‑1 (MQ=255)
                           cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                         <  1:1368559/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:3266232/42‑1 (MQ=255)
                              gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTg                                                                   <  1:915770/42‑1 (MQ=255)
                              gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                         >  1:41487/1‑68 (MQ=255)
                                   aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCaaa                                                          >  1:3445416/1‑46 (MQ=255)
                                   aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                  <  1:331130/70‑1 (MQ=255)
                                   aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                  <  1:3590668/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:4013411/1‑70 (MQ=255)
                                        cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                             >  1:3527105/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:2004264/71‑1 (MQ=255)
                                             aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                       >  1:2004543/1‑71 (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:2751769/71‑1 (MQ=255)
                                                  cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                  >  1:789520/1‑71 (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)
                                                         tCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt               >  1:2744971/1‑67 (MQ=255)
                                                         tCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa           >  1:2685177/1‑71 (MQ=255)
                                                          cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                       >  1:485610/1‑42 (MQ=255)
                                                          cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                       >  1:3043588/1‑42 (MQ=255)
                                                          cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                       >  1:733461/1‑42 (MQ=255)
                                                           gTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                       >  1:956752/1‑41 (MQ=255)
                                                           gTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                  >  1:2624806/1‑62 (MQ=255)
                                                            ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc        <  1:3418106/71‑1 (MQ=255)
                                                            ttGAACGCTTTGCCGGGCAAATGCTCCTCCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc        <  1:2043968/71‑1 (MQ=255)
                                                             tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa           >  1:1535395/1‑67 (MQ=255)
                                                             tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa           >  1:2506129/1‑67 (MQ=255)
                                                              gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa       <  1:2427816/70‑1 (MQ=255)
                                                              gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt      >  1:2388981/1‑71 (MQ=255)
                                                              gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt      >  1:1358819/1‑71 (MQ=255)
                                                               aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa           <  1:2860077/65‑1 (MQ=255)
                                                                aCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc        >  1:2311023/1‑67 (MQ=255)
                                                                  gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTTGGATGGAACTCATTGGt  <  1:932480/71‑1 (MQ=255)
                                                                   cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa       <  1:2773059/65‑1 (MQ=255)
                                                                   cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg   >  1:3077225/1‑69 (MQ=255)
                                                                   cTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTg    >  1:1958010/1‑68 (MQ=255)
                                                                    tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt  >  1:3385894/1‑69 (MQ=255)
                                                                    |                                                                    
TCGAGCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGT  >  minE/2919172‑2919308

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: