Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,267 G→C 13.9% 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,2670GC13.9% 222.0 / 22.2 86L143L (CTC→CTG) ‡argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (4/8);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCG  >  minE/2919198‑2919336
                                                                     |                                                                     
aCGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTg                                                                       <  1:3118263/70‑1 (MQ=255)
  gcgcgTCACCTGCATAGACAAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTgag                                                                     >  1:3389424/1‑70 (MQ=255)
       tCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATc                                                               <  1:984147/71‑1 (MQ=255)
           cTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaaga                                                                 <  1:1332913/65‑1 (MQ=255)
             gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                         >  1:3375487/1‑71 (MQ=255)
             gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                         >  1:1443836/1‑71 (MQ=255)
             gCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                         >  1:2729199/1‑71 (MQ=255)
              cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                         <  1:2524156/70‑1 (MQ=255)
              cATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcag                                                         >  1:215778/1‑70 (MQ=255)
                tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                      >  1:3820175/1‑71 (MQ=255)
                tAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCcagcag                                                      >  1:1566546/1‑71 (MQ=255)
                 aGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCagcagc                                                     >  1:345191/1‑71 (MQ=255)
                  gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                    <  1:678766/71‑1 (MQ=255)
                  gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                    <  1:488704/71‑1 (MQ=255)
                  gACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCt                                                    <  1:1845784/71‑1 (MQ=255)
                  gACCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                    >  1:2324041/1‑71 (MQ=255)
                   aCCAGCGTCATCTCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCcagca                                                       <  1:1536491/67‑1 (MQ=255)
                    ccAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGc                                                  <  1:417966/71‑1 (MQ=255)
                     cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGagaa                                                                   >  1:3417191/1‑53 (MQ=255)
                     cAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCg                                                 >  1:898756/1‑71 (MQ=255)
                          tcctttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                            <  1:1672665/68‑1 (MQ=255)
                          tCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGa                                            <  1:1495109/71‑1 (MQ=255)
                           cATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGAt                                           >  1:3709877/1‑71 (MQ=255)
                            atctCGTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTg                                          <  1:2544862/68‑1 (MQ=255)
                            aTTTCGTTGAACGCTTTGCCTGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                          <  1:1517941/71‑1 (MQ=255)
                             tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGTGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         <  1:2800/71‑1 (MQ=255)
                             tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         <  1:291035/71‑1 (MQ=255)
                             tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         <  1:1516533/71‑1 (MQ=255)
                             tttCGTTGAACGCTTTGCCGGGCAAATGCTCCTGAATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         <  1:2940543/71‑1 (MQ=255)
                                cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                      >  1:3500525/1‑71 (MQ=255)
                                cGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAAc                                      >  1:464517/1‑71 (MQ=255)
                                 gTTAAACGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                       <  1:3462692/69‑1 (MQ=255)
                                  ttGAACGCTTTGCCGGGCAAATGGTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         >  1:16463/1‑66 (MQ=255)
                                  ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                    <  1:1671315/71‑1 (MQ=255)
                                  ttGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                    <  1:1014133/71‑1 (MQ=255)
                                   tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGTATGGAACTCa                                   >  1:242989/1‑71 (MQ=255)
                                   tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTc                                    <  1:856089/70‑1 (MQ=255)
                                   tGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCa                                   >  1:632586/1‑71 (MQ=255)
                                    gAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCa                                   <  1:2928175/70‑1 (MQ=255)
                                     aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATgg                                         <  1:3392923/63‑1 (MQ=255)
                                     aaCGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAt                                  <  1:935183/70‑1 (MQ=255)
                                     aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                  <  1:3478967/70‑3 (MQ=255)
                                     aaCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTcgt                                  <  1:1989517/70‑3 (MQ=255)
                                      aCGCCTTGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTgg                                                  >  1:705940/1‑52 (MQ=255)
                                       cGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGAt                                                                >  1:1315962/1‑38 (MQ=255)
                                        gCTTTGCCGGGCAATTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                              <  1:2727425/71‑1 (MQ=255)
                                        gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                              <  1:56103/71‑1 (MQ=255)
                                        gCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                              <  1:2405156/71‑1 (MQ=255)
                                          tttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGt                              >  1:25113/1‑69 (MQ=255)
                                          cttGCCCGGCAGGTGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTCGAACTCGTTGGt                              >  1:260038/2‑69 (MQ=255)
                                           ttGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                           >  1:1152492/1‑71 (MQ=255)
                                            tGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAgg                          >  1:983456/1‑71 (MQ=255)
                                             gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                         >  1:1081914/1‑71 (MQ=255)
                                             gCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGc                         >  1:356765/1‑71 (MQ=255)
                                               cGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGcc                        <  1:1194635/70‑1 (MQ=255)
                                                gggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAt                      <  1:2248350/71‑1 (MQ=255)
                                                 ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                     <  1:3875150/71‑1 (MQ=255)
                                                 ggCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCAtt                     <  1:3517193/71‑1 (MQ=255)
                                                  gCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                    >  1:2506764/1‑71 (MQ=255)
                                                   cAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTc                    >  1:3647491/1‑70 (MQ=255)
                                                   cAAATGCTCCTGCATGGTGAGAAGATCCGCAAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                   <  1:2330912/71‑1 (MQ=255)
                                                    aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTcc                   >  1:1394035/1‑70 (MQ=255)
                                                    aaaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCa                  >  1:3419366/1‑71 (MQ=255)
                                                     aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                        >  1:301951/1‑48 (MQ=255)
                                                     aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGa                                        >  1:2164908/1‑48 (MQ=255)
                                                     aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCAtt                                 >  1:2429364/1‑55 (MQ=255)
                                                     aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                 <  1:1589273/71‑1 (MQ=255)
                                                     aaTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCAt                 <  1:2438395/71‑1 (MQ=255)
                                                      gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                    <  1:2551466/34‑1 (MQ=16)
                                                      gtGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCt                                                    <  1:2163477/34‑1 (MQ=16)
                                                      aTGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATa                <  1:738670/71‑1 (MQ=255)
                                                       tGCTCCTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGaa                                       >  1:278211/1‑47 (MQ=16)
                                                        gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATg                                          >  1:3463126/1‑43 (MQ=255)
                                                        gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc              >  1:1154139/1‑71 (MQ=255)
                                                        gCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATAcc              >  1:2131374/1‑71 (MQ=255)
                                                          tCCTGCATGGTGAGAAGATCCGCCAGCAGGTGCGTGGGATGGAACTCATTGGTCAgg                          >  1:1842506/1‑57 (MQ=255)
                                                          tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                           <  1:3008816/56‑1 (MQ=255)
                                                          tCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAg                           <  1:4085/56‑1 (MQ=255)
                                                           ccTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATTGAACTCATTGGTCAg                           <  1:3099690/55‑1 (MQ=255)
                                                            cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCa          >  1:2599004/1‑71 (MQ=255)
                                                            cTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGGCAGGCCATTCCATACCgg            <  1:2504979/69‑1 (MQ=255)
                                                            cTGCATGGTCATCAGGTCCGCCAGCAGCTGGGTCGGGTGGAACTCGTTGGTCAGCCCGTTCCACACCGGc           <  1:2359793/70‑1 (MQ=255)
                                                              gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg        >  1:254901/1‑71 (MQ=255)
                                                              gCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACg        >  1:1567680/1‑71 (MQ=255)
                                                               cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc       <  1:1445464/71‑1 (MQ=255)
                                                               cATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc       <  1:985564/71‑1 (MQ=255)
                                                                aTGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGc       >  1:1852288/1‑70 (MQ=255)
                                                                    tGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCg  <  1:236448/71‑1 (MQ=255)
                                                                     |                                                                     
ACGCGCGTCACCTGCATAGACCAGCGTCATTTCGTTGAACGCTTTGCCGGGCAAATGCTCCTGCATGGTGAGAAGATCCGCCAGCAGCTGCGTGGGATGGAACTCATTGGTCAGGCCATTCCATACCGGCACGCTAGCG  >  minE/2919198‑2919336

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: