Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,246 G→C 11.2% P150P (CCC→CCG argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2460GC11.2% 244.3 / 18.1 89P150P (CCC→CCGargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (49/30);  new base C (4/6);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCC  >  minE/2919176‑2919314
                                                                      |                                                                    
gCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCg                                                                      <  1:902106/71‑1 (MQ=255)
     gagttGCCCATGTTGTTGCGCGCATCGCCCGCGTAGACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCa                                                                  <  1:1098310/67‑1 (MQ=255)
     gAATTGCCCATGTTGTTACGCTCGTCACCTGCATAGACCAGGGTCATTTCGTTGAACGCTTTGCCGGGCaa                                                                 <  1:4103696/71‑1 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                                  <  1:598495/70‑1 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                                  <  1:132593/70‑1 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                                  >  1:481889/1‑70 (MQ=255)
     gAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCa                                                                  <  1:2921634/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: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)
                 ttgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTc                                                           >  1:1663463/1‑65 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                      >  1:359370/1‑69 (MQ=255)
                  tgttACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCa                                                      >  1:364041/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)
                          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:2325574/1‑71 (MQ=255)
                               aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCg                                       >  1:789079/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: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:3168188/71‑1 (MQ=255)
                                                  aTTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                    >  1:990945/1‑71 (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:733461/1‑42 (MQ=255)
                                                      cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                             >  1:3043588/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:2506129/1‑67 (MQ=255)
                                                         tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                 >  1:1535395/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)
                                                                 ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg     >  1:1611566/1‑71 (MQ=255)
                                                                  tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa      <  1:974414/69‑1 (MQ=255)
                                                                  tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg    >  1:2211170/1‑71 (MQ=255)
                                                                    ccGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc  <  1:757283/71‑1 (MQ=255)
                                                                     cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa      >  1:2964642/1‑66 (MQ=255)
                                                                     cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCa      >  1:3480953/1‑66 (MQ=255)
                                                                      |                                                                    
GCATCGAATTGCCCATGTTGTTACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCC  >  minE/2919176‑2919314

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: