Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,157 G→T 7.6% V221V (GTG→GTT carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,1570GT 7.6% 258.8 / 4.4 78V221V (GTG→GTTcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (28/44);  new base T (4/2);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCTCCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAC  >  minE/29089‑29225
                                                                    |                                                                    
tctcCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCg                                                                    >  1:4109938/1‑71 (MQ=255)
    cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                                 >  1:361581/1‑70 (MQ=255)
    cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                                 >  1:437485/1‑70 (MQ=255)
    cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                                <  1:3390540/71‑1 (MQ=255)
    cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                                <  1:1854292/71‑1 (MQ=255)
    cGACCAAAGAGTGGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                                 >  1:281821/1‑70 (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:1340825/71‑1 (MQ=255)
          aaGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAAc                                                          <  1:1370709/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:3982518/1‑71 (MQ=255)
                        gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                            >  1:2212175/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: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:1222667/70‑1 (MQ=255)
                                       ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                              <  1:2626890/70‑1 (MQ=255)
                                       ggCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCg                              <  1:3962498/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:893095/71‑1 (MQ=255)
                                              aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                      <  1:280784/71‑1 (MQ=255)
                                              aaGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                      <  1:2481305/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:3349105/71‑1 (MQ=255)
                                                      gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat              <  1:1538002/71‑1 (MQ=255)
                                                      gagaTGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat              <  1:3920045/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: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:1955758/71‑1 (MQ=255)
                                                                  gtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCacac  <  1:2972940/71‑1 (MQ=255)
                                                                  gtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCacac  <  1:1055805/71‑1 (MQ=255)
                                                                    |                                                                    
TCTCCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATCCACAC  >  minE/29089‑29225

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: