Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 29,154 G→C 7.7% V220V (GTG→GTC carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE29,1540GC 7.7% 229.1 / 7.8 78V220V (GTG→GTCcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (28/44);  new base C (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.26e-01

TCTCTCTCCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATC  >  minE/29085‑29220
                                                                     |                                                                  
tctctctcCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAgtgg                                                                   >  1:1365814/1‑71 (MQ=255)
    tctcCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCg                                                               >  1:4109938/1‑71 (MQ=255)
        cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                            >  1:437485/1‑70 (MQ=255)
        cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                            >  1:361581/1‑70 (MQ=255)
        cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                           <  1:1854292/71‑1 (MQ=255)
        cGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGAt                                                           <  1:3390540/71‑1 (MQ=255)
        cGACCAAAGAGTGGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGa                                                            >  1:281821/1‑70 (MQ=255)
            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)
                   ttGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAgtg                                                                    <  1:60717/51‑1 (MQ=255)
                   ttGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAgtg                                                                    <  1:1689759/51‑1 (MQ=255)
                    tGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAgtg                                                                    <  1:380649/50‑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)
                           tgatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAgtgg                                                                   >  1:386825/1‑44 (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:2476285/1‑71 (MQ=255)
                            gatgaGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGcatcat                                       >  1:2212175/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:2103628/1‑71 (MQ=255)
                                 gTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCt                                  >  1:1875249/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: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: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: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:1451196/71‑1 (MQ=255)
                                                                 aaGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATc  >  1:3774596/1‑71 (MQ=255)
                                                                  aGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTc                 <  1:2943821/55‑1 (MQ=255)
                                                                    tggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat         <  1:3791995/61‑1 (MQ=255)
                                                                    tggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat         <  1:3520752/61‑1 (MQ=255)
                                                                    tggtgCGTGATAAAAACCACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat         <  1:1179256/61‑1 (MQ=255)
                                                                     ggtgCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTcgatgcgat         <  1:2983477/60‑1 (MQ=255)
                                                                     |                                                                  
TCTCTCTCCGACCAAAGAGTTGCTGATTGATGAGTCGCTGATCGGCTGGAAAGAGTACGAGATGGAAGTGGTGCGTGATAAAAACGACAACTGCATCATCGTCTGCTCTATCGAAAACTTCGATGCGATGGGCATC  >  minE/29085‑29220

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: