Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,267 G→C 7.1% L143L (CTC→CTG) ‡ argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,2670GC 7.1% 254.5 / 7.9 84L143L (CTC→CTG) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (42/36);  new base C (4/2);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACG  >  minE/2919199‑2919330
                                                                    |                                                               
cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                               <  1:1267869/71‑1 (MQ=255)
cgcgcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                               <  1:1368559/71‑1 (MQ=255)
   gcgTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                               >  1:41487/1‑68 (MQ=255)
        aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                                        <  1:3590668/70‑1 (MQ=255)
        aCCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATcc                                                        <  1:331130/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)
            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)
                  aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                             <  1:2004264/71‑1 (MQ=255)
                  aCCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTg                                             >  1:2004543/1‑71 (MQ=255)
                    cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                           <  1:85451/71‑1 (MQ=255)
                      gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg                                         >  1:1666852/1‑71 (MQ=255)
                      gCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTg                                         <  1:1871423/71‑1 (MQ=255)
                       cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                        <  1:2751769/71‑1 (MQ=255)
                       cGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                        >  1:789520/1‑71 (MQ=255)
                        gTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTgg                                        >  1:3600264/1‑70 (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: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)
                                          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)
                                                ggCAAATGCTCCTGCATGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc              <  1:545322/71‑1 (MQ=255)
                                                 gCAAATGCTCCTGCATTGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc              <  1:2475523/71‑1 (MQ=255)
                                                 gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc              <  1:3169777/71‑1 (MQ=255)
                                                 gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc              <  1:2361766/71‑1 (MQ=255)
                                                  cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc              >  1:2085531/1‑70 (MQ=255)
                                                   aggtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCa            >  1:98435/4‑71 (MQ=255)
                                                   aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTg                          >  1:3569019/1‑57 (MQ=255)
                                                   aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc             <  1:1389860/70‑1 (MQ=255)
                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                   <  1:582538/47‑1 (MQ=255)
                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                   <  1:1632614/47‑1 (MQ=255)
                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                   <  1:2956864/47‑1 (MQ=255)
                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa            >  1:1553601/1‑70 (MQ=255)
                                                    aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa            <  1:2202348/70‑1 (MQ=255)
                                                     aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                         <  1:228071/56‑1 (MQ=255)
                                                     aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTgg                         <  1:425009/56‑1 (MQ=255)
                                                      tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa            >  1:2689462/1‑68 (MQ=255)
                                                      tGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAc         <  1:1137237/71‑1 (MQ=255)
                                                      tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAgcc                   >  1:1982011/1‑59 (MQ=255)
                                                       gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACt                               <  1:1621294/48‑1 (MQ=255)
                                                       gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCATGCCATTCCATAcc        >  1:647500/1‑71 (MQ=255)
                                                           cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCgg      <  1:54855/69‑1 (MQ=255)
                                                           cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGc     <  1:2539398/70‑1 (MQ=255)
                                                           cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCgg      <  1:1283027/69‑1 (MQ=255)
                                                            tGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                   <  1:3690060/39‑1 (MQ=255)
                                                             gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCAc   >  1:3675509/1‑70 (MQ=255)
                                                             gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg  >  1:2395215/1‑71 (MQ=255)
                                                             gCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGtt                           >  1:1988023/1‑46 (MQ=12)
                                                              cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                 >  1:1828206/1‑39 (MQ=255)
                                                              cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                 >  1:321127/1‑39 (MQ=255)
                                                              cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                 >  1:4136648/1‑39 (MQ=255)
                                                              cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCAc   <  1:1229525/69‑1 (MQ=255)
                                                                  gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGaa                                 >  1:716041/1‑35 (MQ=255)
                                                                  gTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                   >  1:2418342/1‑49 (MQ=255)
                                                                    |                                                               
CGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACG  >  minE/2919199‑2919330

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: