Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,166 A→G 7.3% K224K (AAA→AAG carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,1660AG 7.3% 229.7 / 9.2 82K224K (AAA→AAGcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (27/49);  new base G (4/2);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGGTGACTCCA  >  minE/29097‑29236
                                                                     |                                                                      
cAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaa                                                                       <  1:3348318/71‑1 (MQ=255)
cAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaa                                                                       <  1:2083262/71‑1 (MQ=255)
cAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaa                                                                        >  1:3231441/1‑70 (MQ=255)
  aaGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaa                                                                       <  1:566721/69‑1 (MQ=255)
  aaGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAAc                                                                     <  1:1370709/71‑1 (MQ=255)
  aaGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAAc                                                                     <  1:1340825/71‑1 (MQ=255)
    gagTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGa                                                                   >  1:6795/1‑71 (MQ=255)
      gTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACa                                                                 <  1:1411020/71‑1 (MQ=255)
          cTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGAc                                                                  <  1:3319799/66‑1 (MQ=255)
           tgattgatGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACa                                                                 <  1:2229300/66‑1 (MQ=255)
           tgattgatGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACa                                                                 <  1:163292/66‑1 (MQ=255)
             attgatGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTTATAAAAACGACAACTGc                                                            <  1:3742618/69‑1 (MQ=255)
             attgatGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACaa                                                                >  1:1420460/1‑65 (MQ=255)
               tgatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAAc                                                                     <  1:3317612/58‑1 (MQ=255)
               tgatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGAc                                                                  >  1:2471959/1‑61 (MQ=255)
               tgatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGGGGTGCGTGATAAAAAc                                                                     <  1:1353721/58‑1 (MQ=255)
                gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                                       >  1:2212175/1‑71 (MQ=255)
                gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                                       >  1:3982518/1‑71 (MQ=255)
                gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                                       >  1:2476285/1‑71 (MQ=255)
                gaagaGTCGCTGCTGGGCTGGAAAGAGTACGAGATGGAAGTCGTTCGCGACAAGAACGACAACTGca                                                           <  1:1682420/64‑1 (MQ=255)
                gaagaGTCGCTGCTGGGCTGGAAAGAGTACGAGATGGAAGTCGTTCGCGACAAGAACGACAAATGca                                                           <  1:3391729/64‑1 (MQ=255)
                 atgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGa                                                                   >  1:657205/1‑58 (MQ=255)
                  tgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCg                                                     <  1:1808405/71‑1 (MQ=255)
                  agaGTCGCTGCTGGGCTGGAAAGAGTACGAGATGGAAGTCGTTCGCGACAAGAACGACAACTGcatc                                                         >  1:384156/2‑67 (MQ=255)
                  agaGTCGCTGCTGGGCTGGAAAGAGTACGAGATGGAAGTCGTTCGCGACAAGAACGACAACTGcatc                                                         >  1:1019702/2‑67 (MQ=255)
                   gaGTCGCTGATCGGCTGGAACGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCg                                                     <  1:1504168/70‑1 (MQ=255)
                     gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                                  >  1:1875249/1‑71 (MQ=255)
                     gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                                  >  1:1069807/1‑71 (MQ=255)
                     gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                                  >  1:2103628/1‑71 (MQ=255)
                        gCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaa                                                                        <  1:3196388/46‑1 (MQ=255)
                        gCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACaa                                                                >  1:1969373/1‑54 (MQ=255)
                            aTCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGAc                                                                  <  1:192789/48‑1 (MQ=255)
                            aTCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTa                                            <  1:702284/70‑1 (MQ=255)
                            aTCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTAt                                           >  1:2058862/1‑71 (MQ=255)
                               ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaaa                                                                      >  1:2972137/1‑41 (MQ=255)
                               ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaaa                                                                      >  1:422497/1‑41 (MQ=255)
                               ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                         <  1:3962498/70‑1 (MQ=255)
                               ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                         <  1:1222667/70‑1 (MQ=255)
                               ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                         <  1:2626890/70‑1 (MQ=255)
                                gCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGaa                                       <  1:706362/71‑1 (MQ=255)
                                 cTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTc                                                   >  1:939821/1‑58 (MQ=255)
                                  tGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaa                                                                        <  1:1245005/36‑1 (MQ=255)
                                   ggAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGct                                               <  1:2608360/60‑1 (MQ=255)
                                      aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGt                                                    >  1:3152400/1‑52 (MQ=255)
                                      aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTc                                                   >  1:213725/1‑53 (MQ=255)
                                      aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                                 <  1:2481305/71‑1 (MQ=255)
                                      aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                                 <  1:893095/71‑1 (MQ=255)
                                      aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                                 <  1:280784/71‑1 (MQ=255)
                                       agagTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcg                                <  1:1946519/71‑1 (MQ=255)
                                       agagTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                                 >  1:3073476/1‑70 (MQ=255)
                                         agTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgat                              >  1:1266611/1‑71 (MQ=255)
                                         agTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgat                              >  1:2310295/1‑71 (MQ=255)
                                           tccGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgc                            <  1:1511192/69‑1 (MQ=255)
                                            aCGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgc                            <  1:983501/70‑1 (MQ=255)
                                            aCGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgc                            <  1:2066933/70‑1 (MQ=255)
                                             cGAGATTGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcga                          <  1:1265662/71‑1 (MQ=255)
                                             cGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcga                          <  1:570091/71‑1 (MQ=255)
                                              gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:3349105/71‑1 (MQ=255)
                                              gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:3920045/71‑1 (MQ=255)
                                              gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:1538002/71‑1 (MQ=255)
                                               agaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgatg                        >  1:1206870/1‑71 (MQ=255)
                                                gaTTGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:194416/69‑1 (MQ=255)
                                                gaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg                       <  1:2166276/71‑1 (MQ=255)
                                                gaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg                       <  1:3567857/71‑1 (MQ=255)
                                                gaTGGAAGTCGTTCGCGACAAGAACGACAACTGCATCATCATCTGCTCGATCGAGAACCTggatccga                          >  1:2434487/1‑68 (MQ=255)
                                                gaTGGAAGTCGTTCGCGACAAGAACGACAACTGCATCATCATCTGCTCGATCAAGAACCTggatccga                          >  1:792767/1‑68 (MQ=255)
                                                 aTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg                       <  1:2171051/70‑1 (MQ=255)
                                                   ggAAGTGGTGCTTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:560068/66‑1 (MQ=255)
                                                   ggAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCa                    >  1:2935347/1‑71 (MQ=255)
                                                     aaGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATc                  <  1:1451196/71‑1 (MQ=255)
                                                     aaGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATc                  >  1:3774596/1‑71 (MQ=255)
                                                      aGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                                 <  1:2943821/55‑1 (MQ=255)
                                                        tggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:3520752/61‑1 (MQ=255)
                                                        tggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:3791995/61‑1 (MQ=255)
                                                        tggtgCGTGATAAAAACCACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:1179256/61‑1 (MQ=255)
                                                         ggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                         <  1:2983477/60‑1 (MQ=255)
                                                          gtgCGTGATAAAAACGACAACTGCATCATCGTCTTCTCTATCGAAAACTTcgatgcgat                         <  1:659271/59‑1 (MQ=255)
                                                          gtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCacac             <  1:1055805/71‑1 (MQ=255)
                                                          gtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCacac             <  1:2972940/71‑1 (MQ=255)
                                                          gtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCacac             <  1:1955758/71‑1 (MQ=255)
                                                             cGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACt                                   <  1:3742816/46‑1 (MQ=255)
                                                             cGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACt                                   <  1:2589048/46‑1 (MQ=255)
                                                             cGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgatg                        >  1:4003132/1‑57 (MQ=255)
                                                             cGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg                       >  1:1295749/1‑58 (MQ=255)
                                                               tGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTctatgc                            <  1:2934962/51‑1 (MQ=255)
                                                               tGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCgg          <  1:3371505/69‑1 (MQ=255)
                                                                 aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGaaaa                                     >  1:1472897/1‑40 (MQ=255)
                                                                 aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAcc            >  1:3919952/1‑65 (MQ=255)
                                                                 aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAcc            >  1:2054986/1‑65 (MQ=255)
                                                                 aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAcc            >  1:2878055/1‑65 (MQ=255)
                                                                   aaaaaCGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGGTGACTc    >  1:3786070/1‑71 (MQ=255)
                                                                    aaaaCGACAACTGCATCATCGTCTGCTCTATCGaa                                       <  1:153407/35‑1 (MQ=255)
                                                                     aaaCGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg                       >  1:3641017/1‑50 (MQ=255)
                                                                     aaaCGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGGTGACTCCa  <  1:484221/71‑1 (MQ=255)
                                                                     |                                                                      
CAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGGTGACTCCA  >  minE/29097‑29236

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: