Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,163 T→C 6.7% D223D (GAT→GAC carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,1630TC 6.7% 253.5 / 9.2 90D223D (GAT→GACcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (30/54);  new base C (4/2);  total (34/56)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGG  >  minE/29093‑29228
                                                                      |                                                                 
cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                                   <  1:1854292/71‑1 (MQ=255)
cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                                   <  1:3390540/71‑1 (MQ=255)
    cAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaa                                                               <  1:2083262/71‑1 (MQ=255)
    cAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaa                                                               <  1:3348318/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:2476285/1‑71 (MQ=255)
                    gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                               >  1:2212175/1‑71 (MQ=255)
                    gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                               >  1:3982518/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:1019702/2‑67 (MQ=255)
                      agaGTCGCTGCTGGGCTGGAAAGAGTACGAGATGGAAGTCGTTCGCGACAAGAACGACAACTGcatc                                                 >  1:384156/2‑67 (MQ=255)
                       gaGTCGCTGATCGGCTGGAACGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCg                                             <  1:1504168/70‑1 (MQ=255)
                         gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                          >  1:1875249/1‑71 (MQ=255)
                         gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                          >  1:2103628/1‑71 (MQ=255)
                         gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                          >  1:1069807/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:422497/1‑41 (MQ=255)
                                   ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATaaaaa                                                              >  1:2972137/1‑41 (MQ=255)
                                   ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                 <  1:3962498/70‑1 (MQ=255)
                                   ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                 <  1:2626890/70‑1 (MQ=255)
                                   ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                                 <  1:1222667/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:2310295/1‑71 (MQ=255)
                                             agTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgat                      >  1:1266611/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:1538002/71‑1 (MQ=255)
                                                  gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                 <  1:3920045/71‑1 (MQ=255)
                                                  gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                 <  1:3349105/71‑1 (MQ=255)
                                                   agaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgatg                >  1:1206870/1‑71 (MQ=255)
                                                    gaTTGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat                 <  1:194416/69‑1 (MQ=255)
                                                    gaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg               <  1:3567857/71‑1 (MQ=255)
                                                    gaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATgg               <  1:2166276/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:3774596/1‑71 (MQ=255)
                                                         aaGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATc          <  1:1451196/71‑1 (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:2054986/1‑65 (MQ=255)
                                                                     aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAcc    >  1:3919952/1‑65 (MQ=255)
                                                                     aTAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAcc    >  1:2878055/1‑65 (MQ=255)
                                                                      |                                                                 
CGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACACCGG  >  minE/29093‑29228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: