Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,136,538 T→C 39.7% K129K (AAA→AAG yniB ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,136,5380TC39.7% 40.6 / 75.6 68K129K (AAA→AAGyniBpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/24);  new base C (13/14);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGG  >  minE/1136472‑1136608
                                                                  |                                                                      
gAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTgg                                                                    >  1:4333386/1‑71 (MQ=255)
gAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTgg                                                                    <  1:830604/71‑1 (MQ=255)
  aaaTAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTGACGCGTTAGCCCCTCTTCCCCTTTGGCt                                                                  <  1:2514627/71‑1 (MQ=255)
    aTAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCtct                                                                <  1:1839950/71‑3 (MQ=255)
    aTAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCATCTTCCCCCTTGGCttt                                                                <  1:3430574/71‑1 (MQ=255)
       gTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCttttt                                                              <  1:1808952/70‑1 (MQ=255)
        tATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCttttt                                                              <  1:3788931/69‑1 (MQ=255)
          tGGTGCGGAACAACAATACGGGACTCAATTTGTTCCCGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTc                                                             >  1:2274411/1‑68 (MQ=255)
          tGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAg                                                          <  1:331186/71‑1 (MQ=255)
           ggTGCGGAACAACAATACGGGACTCAATTTGTTCACGGGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAg                                                          <  1:1303004/70‑1 (MQ=255)
           ggTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAg                                                          <  1:1880778/70‑1 (MQ=255)
               cGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCa                                                           >  1:2128084/1‑65 (MQ=255)
                ggAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCTCCCTTGGCTTTTTCCAGAATGa                                                     <  1:68959/70‑1 (MQ=255)
                  aacaacAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAg                                                    >  1:3182386/1‑69 (MQ=255)
                   acaacaATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTg                                                 <  1:442104/71‑1 (MQ=255)
                       caaTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCtttt                                                               >  1:4681897/1‑53 (MQ=255)
                       caaTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCt                                             >  1:2045915/1‑71 (MQ=255)
                       caaTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCt                                             >  1:4212213/1‑71 (MQ=255)
                        aaTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGaa                                                        >  1:4642585/1‑59 (MQ=255)
                         aTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTc                                           <  1:1004719/71‑1 (MQ=255)
                         aTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTc                                           <  1:3517158/71‑1 (MQ=255)
                         aTACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTc                                           <  1:2106881/71‑1 (MQ=255)
                          tACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTc                                           <  1:3407409/70‑1 (MQ=255)
                          tACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGaa                                                        >  1:363458/1‑57 (MQ=255)
                          tACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGaa                                                        >  1:2176401/1‑57 (MQ=255)
                          tACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGaa                                                        >  1:383360/1‑57 (MQ=255)
                          tACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCa                                          <  1:2036884/71‑1 (MQ=255)
                             gggACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCaa                                         <  1:4247959/69‑1 (MQ=255)
                             gggACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATCAGTTGATCTTCaa                                         <  1:2479995/69‑1 (MQ=255)
                              ggACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGAt                                               <  1:4637335/62‑1 (MQ=255)
                                aCTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCt                                    <  1:2921203/71‑1 (MQ=255)
                                  tCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGc                                                                   <  1:1080943/38‑1 (MQ=255)
                                  tCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGc                                                                   <  1:330124/38‑1 (MQ=255)
                                  tCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAg                                                          <  1:2762621/47‑1 (MQ=255)
                                  tCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTg                                                 <  1:3483676/56‑1 (MQ=255)
                                  tCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTc                                  >  1:4565706/1‑71 (MQ=255)
                                   cAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTcc                                 >  1:3145928/1‑71 (MQ=255)
                                   cAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTc                                  >  1:1039764/1‑70 (MQ=255)
                                    aaTTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTg                                                 <  1:1954897/54‑1 (MQ=255)
                                    aaTTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCaa                                         >  1:1238011/1‑62 (MQ=255)
                                    aaTTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTcc                                 >  1:2720687/1‑70 (MQ=255)
                                    aaTTTGTTCACGCGTTAGCCCCTCTTCCCCCTTgg                                                                    >  1:4567393/1‑35 (MQ=255)
                                      tttGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAg                                                          <  1:2491188/43‑1 (MQ=255)
                                       ttGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGc                              >  1:3355728/1‑70 (MQ=255)
                                       ttGTTCACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGc                              >  1:2154554/1‑70 (MQ=255)
                                          ttCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCaag                           <  1:360923/70‑1 (MQ=255)
                                           tctcgcgTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCagga                          <  1:95016/67‑1 (MQ=255)
                                            cACGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGaa                                                        >  1:1596760/1‑39 (MQ=255)
                                            cACGCGTTAGCCACTCTTCCCCTTTGGCTTTTTCCAGaa                                                        >  1:3966245/1‑39 (MQ=255)
                                             aCGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAAt                                                       <  1:3509533/39‑1 (MQ=255)
                                             aCGCGTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAAt                                                       <  1:2438034/39‑1 (MQ=255)
                                               gcgTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATg                                                      >  1:3661442/1‑38 (MQ=255)
                                               gcgTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATc                                              >  1:2178903/1‑46 (MQ=255)
                                                 gTTAGCCCCTCTTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAgcgcgc                   >  1:824442/1‑71 (MQ=255)
                                                       cccTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCTCAAGAAGCGCGCCTGGCg             <  1:2551347/71‑1 (MQ=255)
                                                       cccTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCg             <  1:17194/71‑1 (MQ=255)
                                                         ctctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGAc           >  1:9170/1‑71 (MQ=255)
                                                         ctctTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCt                                    <  1:1621011/46‑1 (MQ=255)
                                                         ctctTCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGAc           >  1:3963412/1‑71 (MQ=255)
                                                          tctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCt                                             <  1:1032232/36‑1 (MQ=255)
                                                          tctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTc                                           <  1:847534/38‑1 (MQ=255)
                                                          tctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGAc           >  1:2257440/1‑70 (MQ=255)
                                                           ctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAgcgcgc                   >  1:566958/1‑61 (MQ=255)
                                                           ctTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAgcgcgc                   >  1:871690/1‑61 (MQ=255)
                                                            ttCCCCCTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCt                                    >  1:2327503/1‑43 (MQ=255)
                                                             tCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCaag                           <  1:2201357/51‑1 (MQ=255)
                                                              ccccTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCAt       <  1:144220/70‑1 (MQ=255)
                                                              ccccTTTGGCTTTTTCCAGAATGAGTTAATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCAt       <  1:3725476/70‑1 (MQ=255)
                                                                  tttGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCggg  >  1:1043158/1‑71 (MQ=255)
                                                                  |                                                                      
GAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGTTAGCCCCTCTTCCCCTTTGGCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGG  >  minE/1136472‑1136608

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: