Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,404 C→T 8.5% V538V (GTC→GTT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,4040CT 8.5% 210.8 / 7.4 70V538V (GTC→GTTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/4);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATG  >  minE/386337‑386473
                                                                   |                                                                     
ctgcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTc                                                                       <  1:3246299/68‑1 (MQ=255)
   cTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCg                                                                 >  1:2941614/1‑71 (MQ=255)
   cTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCg                                                                 >  1:2749834/1‑71 (MQ=255)
    tAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCg                                                                 >  1:3860843/1‑70 (MQ=255)
    tAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGc                                                                >  1:912283/1‑71 (MQ=255)
    tAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGc                                                                >  1:501940/1‑71 (MQ=255)
    tAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGc                                                                >  1:3865320/1‑71 (MQ=255)
      aCCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTg                                                              <  1:1074402/71‑1 (MQ=255)
       ccTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTgg                                                             >  1:3942557/1‑71 (MQ=255)
       ccTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTgg                                                             >  1:3909612/1‑71 (MQ=255)
          gctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGggcg                                                          >  1:2721759/1‑71 (MQ=255)
           ctgGCAACAGCGCTGTGCGTACTGGCGTGGGGGGACTTCCTCCATCAGGGCGTGGTCGATCCGCTGggcg                                                          <  1:3777253/70‑1 (MQ=255)
              gCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCAt                                                      <  1:3093336/71‑1 (MQ=255)
              gCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCAt                                                      <  1:22034/71‑1 (MQ=255)
               cAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCAtt                                                     <  1:3042952/71‑1 (MQ=255)
                aaCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTa                                                    <  1:2616200/71‑1 (MQ=255)
                 aCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTaa                                                   <  1:1059878/71‑1 (MQ=255)
                 aCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTaa                                                   <  1:2511833/71‑1 (MQ=255)
                 aCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTaa                                                   <  1:2394814/71‑1 (MQ=255)
                     cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCt                                                               <  1:3911863/55‑1 (MQ=255)
                      gcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActc                                              <  1:1992996/71‑1 (MQ=255)
                          tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGggc                                                           <  1:2940070/54‑1 (MQ=255)
                          tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                            <  1:218731/69‑1 (MQ=255)
                          tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                            <  1:1298947/69‑1 (MQ=255)
                            cGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTgg                                         <  1:3802474/70‑1 (MQ=255)
                            cGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTgg                                         <  1:2854286/70‑1 (MQ=255)
                              tACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCg                                      >  1:3546732/1‑71 (MQ=255)
                               aCTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGAGGCATTaa                                                   >  1:1781670/1‑57 (MQ=255)
                                cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGa                                                                     <  1:3566853/38‑1 (MQ=255)
                                cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGa                                                                     <  1:129059/38‑1 (MQ=255)
                                cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCg                                                                 <  1:2688352/42‑1 (MQ=255)
                                  ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                  >  1:1572174/1‑71 (MQ=255)
                                  ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                  >  1:1556499/1‑71 (MQ=255)
                                   gCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGGGGCATTAACACTCTGTGGCCGCTGtt                                 >  1:3860946/1‑71 (MQ=255)
                                   gCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGtt                                 <  1:2578321/71‑1 (MQ=255)
                                   gCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGtt                                 <  1:1008083/71‑1 (MQ=255)
                                   gCGGGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                            <  1:2022494/60‑1 (MQ=255)
                                     gTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTg                               >  1:3422852/1‑71 (MQ=255)
                                      tGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                              <  1:1384757/71‑1 (MQ=255)
                                      tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCt                                                               <  1:224996/38‑1 (MQ=12)
                                        gggCTACTTCCTCCATCATGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGt                             <  1:697887/70‑1 (MQ=255)
                                        gggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                            >  1:633637/1‑71 (MQ=255)
                                        gggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCCGTGGCCGCTGtt                                 <  1:42016/66‑1 (MQ=255)
                                         ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtgtg                                          >  1:90581/1‑56 (MQ=255)
                                          gCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAtt                          <  1:1400524/71‑1 (MQ=255)
                                            tacttactCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGtt                                 >  1:2705529/1‑62 (MQ=255)
                                            tACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGc                        >  1:1419834/1‑71 (MQ=255)
                                             aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCGGCTGTTTGGTATTGc                        >  1:149171/1‑70 (MQ=255)
                                             aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGc                        >  1:3195728/1‑70 (MQ=255)
                                             aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGc                        >  1:141739/1‑70 (MQ=255)
                                              cTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGcc                       >  1:1986492/1‑70 (MQ=255)
                                               ttCCTCCATCAGGGCGTGGTCGATCCGCCGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCaa                     >  1:2960936/1‑71 (MQ=255)
                                                tcctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCaa                     >  1:2481852/1‑70 (MQ=255)
                                                    ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                    >  1:2955634/1‑67 (MQ=255)
                                                     cATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGAt               <  1:155986/71‑1 (MQ=255)
                                                      aTCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTAtct                         >  1:3356265/1‑58 (MQ=255)
                                                      aTCAGGGCGTGGTCGATCCGCTTGGCGGCATTAACAct                                               <  1:1214192/38‑1 (MQ=255)
                                                      aTCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATg              <  1:665327/71‑1 (MQ=255)
                                                       tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATg              <  1:1210251/70‑1 (MQ=255)
                                                       tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGc             >  1:2737037/1‑71 (MQ=255)
                                                       tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                    >  1:2183047/1‑64 (MQ=255)
                                                       tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGc             <  1:2140766/71‑1 (MQ=255)
                                                           ggCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGc         <  1:452322/71‑1 (MQ=255)
                                                            gCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                              <  1:3584432/49‑1 (MQ=255)
                                                            gCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                              <  1:2442816/49‑1 (MQ=255)
                                                              gTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCt            <  1:1683315/65‑1 (MQ=19)
                                                               tGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCaa                     <  1:2982150/55‑1 (MQ=16)
                                                               tGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAg       <  1:2971277/69‑1 (MQ=255)
                                                                 gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGa    <  1:1116822/70‑1 (MQ=255)
                                                                  tCGATCCGCTGGGCGGCATTAACACTCTGTGGGCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATg  >  1:205137/1‑71 (MQ=255)
                                                                  tCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATg  >  1:1063233/1‑71 (MQ=255)
                                                                   cGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCt                                    <  1:773005/36‑1 (MQ=255)
                                                                   |                                                                     
CTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATG  >  minE/386337‑386473

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: