Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 195,592 G→A 18.0% T79T (ACC→ACT yaeF ← predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE195,5920GA18.0% 170.7 / 28.7 72T79T (ACC→ACTyaeFpredicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (36/23);  new base A (10/3);  total (46/26)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CGTTATTATCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTA  >  minE/195522‑195659
                                                                      |                                                                    
cGTTATTATCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAg                                                                      >  1:2819186/1‑71 (MQ=255)
    attatCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGtt                                                                   >  1:1684319/1‑70 (MQ=255)
    attatCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTa                                                                  >  1:3915835/1‑71 (MQ=255)
    attatCACCCAGAAAGATCGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTa                                                                  >  1:3205393/1‑71 (MQ=255)
        tCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTAccc                                                               >  1:403885/1‑70 (MQ=255)
          aCCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTAcccc                                                              >  1:2360157/1‑69 (MQ=255)
          aCCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTAcccc                                                              >  1:3240370/1‑69 (MQ=255)
          aCCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTAcccc                                                              >  1:1954847/1‑69 (MQ=255)
            ccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAg                                                          >  1:1957409/1‑70 (MQ=255)
            ccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAgg                                                         >  1:3546783/1‑71 (MQ=255)
            ccAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAgg                                                         >  1:348189/1‑71 (MQ=255)
               gAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAg                                                                      <  1:1953243/56‑1 (MQ=255)
                 aaGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAg                                                                      >  1:4015320/1‑54 (MQ=255)
                 aaGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTg                                                     <  1:2465691/70‑1 (MQ=255)
                  aGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgag                                                   <  1:1771730/71‑1 (MQ=255)
                    aTTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCTTgagga                                                 <  1:3802648/71‑1 (MQ=255)
                     ttGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAg                                                                      <  1:35619/50‑1 (MQ=255)
                     ttGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgag                                                   >  1:4128260/1‑68 (MQ=255)
                     ttGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgag                                                   >  1:483203/1‑68 (MQ=255)
                     ttGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCTTgagga                                                 <  1:2156269/70‑1 (MQ=255)
                      tGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTa                                                                  >  1:1951418/1‑53 (MQ=255)
                      tGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgaggag                                                >  1:3589599/1‑70 (MQ=255)
                      tGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgaggag                                                >  1:2709755/1‑70 (MQ=255)
                         aaCGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTa                                                                  >  1:1243353/1‑50 (MQ=255)
                            gTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGc                                          >  1:3421789/1‑70 (MQ=255)
                            gTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCTTGAGGAGAACAGc                                          >  1:741014/1‑70 (MQ=255)
                               aCTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCaaa                                       <  1:2065138/70‑1 (MQ=255)
                                cTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCaaa                                       <  1:650737/69‑1 (MQ=255)
                                 tCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTAcccc                                                              <  1:3644695/46‑1 (MQ=255)
                                   acaGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGcc                                  <  1:3431942/71‑1 (MQ=255)
                                     agagGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACa                                            >  1:4042712/1‑59 (MQ=255)
                                     agagGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACa                                            >  1:3814706/1‑59 (MQ=255)
                                       agGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTgagga                                                 >  1:421198/1‑52 (MQ=255)
                                       agGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCtt                                                      >  1:2750809/1‑47 (MQ=255)
                                       agGAAGTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCtt                                                      >  1:4073991/1‑47 (MQ=255)
                                        ggAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCt                                                       >  1:4005189/1‑45 (MQ=255)
                                         gAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                            >  1:4118584/1‑71 (MQ=255)
                                         gAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGtt                             >  1:3376521/1‑70 (MQ=255)
                                           aGTGCTGAAGACGCTGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          <  1:2248324/71‑1 (MQ=255)
                                           aGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          <  1:1613945/71‑1 (MQ=255)
                                           aGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          <  1:569385/71‑1 (MQ=255)
                                            gTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTg                                                     <  1:353237/43‑1 (MQ=255)
                                            gTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTg                                                     <  1:2625311/43‑1 (MQ=255)
                                            gTGCTGAAGACGCGGATTCCAAATGAAGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTgag                         <  1:1064089/71‑1 (MQ=255)
                                             tGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAgg                                                         >  1:480434/1‑38 (MQ=255)
                                             tGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGcc                                  <  1:776555/61‑1 (MQ=255)
                                              gCTGAAGACGCGGATTCCAAATGAGGTTACCCCaa                                                           >  1:2513833/1‑35 (MQ=255)
                                              gCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGaa                                              >  1:1910507/1‑48 (MQ=255)
                                              gCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                            >  1:2943710/1‑66 (MQ=255)
                                              gCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTgag                         >  1:1252817/1‑69 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTg                                                     >  1:245518/1‑40 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTg                                                     <  1:2657398/40‑1 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          >  1:3408482/1‑67 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          >  1:4128998/1‑67 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATc                      <  1:125595/71‑1 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATc                      <  1:2012137/71‑1 (MQ=255)
                                               cTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATc                      <  1:1431243/71‑1 (MQ=255)
                                                  aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGAt                       >  1:3462184/1‑67 (MQ=255)
                                                  aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGAt                       <  1:1780794/67‑1 (MQ=255)
                                                  aaGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCa                     <  1:2789065/69‑1 (MQ=255)
                                                    gACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          >  1:3058229/1‑62 (MQ=255)
                                                    gACGCGGATTCCAAATGAAGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTg                 >  1:944634/1‑71 (MQ=255)
                                                             tCCAAATGAGGTTACCCCATGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGAt                       <  1:802704/56‑1 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATc                      >  1:3919163/1‑55 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATc                      >  1:930825/1‑55 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCa                     >  1:2984758/1‑56 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTc       <  1:1813397/70‑1 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCt      <  1:3689429/71‑1 (MQ=255)
                                                               cAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCt      >  1:3430828/1‑71 (MQ=255)
                                                               cAAATGAAGTTACCCCAAGGCTTGAGGAGAACAGCAAAt                                      >  1:1585332/1‑39 (MQ=255)
                                                               cAAATGAAGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCg                                 >  1:138421/1‑44 (MQ=255)
                                                                 aaTGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTgag                         <  1:1762329/50‑1 (MQ=255)
                                                                   tGAGGTTACCCCTAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCtttt   <  1:1843147/71‑1 (MQ=255)
                                                                   tGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                            >  1:3628787/1‑45 (MQ=255)
                                                                   tGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTa  >  1:3271512/1‑71 (MQ=255)
                                                                     aGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGttt                            >  1:1024316/1‑43 (MQ=255)
                                                                      ggTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTga                          >  1:1169448/1‑44 (MQ=255)
                                                                      |                                                                    
CGTTATTATCACCCAGAAAGATTGCAACGTGACTCACAGAGGAAGTGCTGAAGACGCGGATTCCAAATGAGGTTACCCCAAGGCTTGAGGAGAACAGCAAATCGCCGGGTTTGAGATCAGGTGCCGTTATTTCTTTTA  >  minE/195522‑195659

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: