Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,763 G→A 79.4% intergenic (‑134/+162) yqeF ← / ← lysS predicted acyltransferase/lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,7630GA79.4% 196.8 / 46.4 102intergenic (‑134/+162)yqeF/lysSpredicted acyltransferase/lysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (9/12);  new base A (30/51);  total (39/63)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 106 total aligned reads displayed.

TAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGT  >  minE/1932695‑1932830
                                                                    |                                                                   
tAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAagg                                                                   >  1:2833520/1‑71 (MQ=255)
 aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAagg                                                                   <  1:2542123/70‑1 (MQ=255)
 aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGt                                                                  >  1:5480343/1‑71 (MQ=255)
  aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAagg                                                                   <  1:1475786/69‑1 (MQ=255)
  aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGt                                                                  <  1:466989/70‑1 (MQ=255)
  aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGtt                                                                 <  1:484226/71‑1 (MQ=255)
  aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGt                                                                  <  1:2724405/70‑1 (MQ=255)
   aaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGttt                                                                <  1:2087578/71‑1 (MQ=255)
    aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTg                                                               <  1:4568146/71‑1 (MQ=255)
    aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTg                                                               <  1:2322765/71‑1 (MQ=255)
    aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTg                                                               >  1:4992381/1‑71 (MQ=255)
     aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                              <  1:4526543/71‑1 (MQ=255)
      aaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                              <  1:1993085/70‑1 (MQ=255)
        gTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAAc                                                           <  1:1081794/71‑1 (MQ=255)
        gTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAAc                                                           <  1:274317/71‑1 (MQ=255)
          gATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAagg                                                                   >  1:205636/1‑61 (MQ=255)
            tttttATCACATAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                       <  1:3988933/71‑1 (MQ=255)
            tttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                       <  1:2596593/71‑1 (MQ=255)
            tttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                       <  1:189235/71‑1 (MQ=255)
             ttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGttt                                                                >  1:945466/1‑61 (MQ=255)
             ttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCaa                                                            >  1:3210075/1‑65 (MQ=255)
              tttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTaa                                                     >  1:4467836/1‑71 (MQ=255)
              tttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                              >  1:4635771/1‑62 (MQ=255)
               ttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAAc                                                    >  1:4384860/1‑71 (MQ=255)
                  tCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAAcc                                                   <  1:224332/69‑1 (MQ=255)
                    acaAAGGAAATATGCCTGAGCACCAGTCAGAGACATAACTGGCACGTAAGGtt                                                                 >  1:5559763/1‑53 (MQ=255)
                      aaaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCa                                                             <  1:774399/55‑1 (MQ=255)
                      aaaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCa                                                             <  1:1151202/55‑1 (MQ=255)
                       aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAAccc                                                  <  1:1527812/65‑1 (MQ=255)
                       aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAAt                                            <  1:7641/71‑1 (MQ=255)
                         ggAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATag                                          <  1:1744577/71‑1 (MQ=255)
                          gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCa                                                         >  1:2459952/1‑55 (MQ=255)
                          gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCa                                                         >  1:2816374/1‑55 (MQ=255)
                          gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGAAAcc                                                          <  1:5080674/54‑1 (MQ=255)
                           aaaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATagag                                        >  1:206936/1‑71 (MQ=255)
                             atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAggg                                      <  1:5524019/71‑1 (MQ=255)
                              tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                              <  1:2963592/46‑1 (MQ=255)
                              tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                       <  1:5629948/53‑1 (MQ=255)
                              tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                       <  1:5674813/53‑1 (MQ=255)
                              tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAgggg                                     <  1:4119660/71‑1 (MQ=255)
                              tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAgggg                                     <  1:5541179/71‑1 (MQ=255)
                                tGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTa                                   <  1:3654125/71‑1 (MQ=255)
                                tGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTa                                   <  1:2967586/71‑1 (MQ=255)
                                tGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTa                                   <  1:4917049/71‑1 (MQ=255)
                                  ccTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                 <  1:4189176/71‑1 (MQ=255)
                                   cTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                 >  1:3184108/1‑70 (MQ=255)
                                    tGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                 >  1:2712129/1‑69 (MQ=255)
                                      agcagcAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                <  1:5269958/68‑1 (MQ=255)
                                      agcagcAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAgg                             >  1:1890582/1‑71 (MQ=255)
                                      agcagcAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACAACTAACCCACCAATAGAGGGGTAGATAg                              >  1:2326976/1‑70 (MQ=255)
                                        cagcagTCAGAGACATAACTGGCACGTAAGGTTTGCAACACTAACcca                                                 >  1:1053511/1‑48 (MQ=38)
                                           cagTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGtt                        <  1:559915/71‑1 (MQ=255)
                                            agTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:845870/71‑1 (MQ=255)
                                            agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:3469119/71‑1 (MQ=255)
                                            agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:3856124/71‑1 (MQ=255)
                                            agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:2676576/71‑1 (MQ=255)
                                            agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:3033175/71‑1 (MQ=255)
                                             gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTaa                      <  1:2553984/71‑1 (MQ=255)
                                             gTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                 <  1:2054094/60‑1 (MQ=255)
                                             gTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:5497369/70‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACTACTAACCCACCAATa                                           <  1:807888/49‑1 (MQ=39)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACccac                                                <  1:1758029/44‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACccac                                                <  1:287773/44‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:4457034/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:2067098/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:1621634/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:5649948/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:1338540/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           <  1:2680932/49‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      >  1:4861620/1‑54 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      >  1:243982/1‑54 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      >  1:1958093/1‑54 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      >  1:4777539/1‑54 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      >  1:813669/1‑54 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                       <  1:667527/69‑1 (MQ=255)
                                              tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAt                     >  1:5142764/1‑71 (MQ=255)
                                                agagACATAACTGGCACGTAGGGTTTGCAACCACTAACccacca                                              >  1:2075483/1‑44 (MQ=255)
                                                agagACATAACTGGCACGTAGGGTTTGCAACCACTAACccacca                                              >  1:507672/1‑44 (MQ=255)
                                                agagACATAACTGGCACGTAGGGTTTGCAACCACTAACccacca                                              >  1:4467516/1‑44 (MQ=255)
                                                agagACATAACTGGCACGTAAGGTTTGCAACCACTAACccacca                                              >  1:3871246/1‑44 (MQ=255)
                                                agagACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTaa                      <  1:1465022/68‑1 (MQ=255)
                                                 gagaCATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAgggg                                     <  1:3092639/52‑1 (MQ=255)
                                                  agaCATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTcc                 <  1:1611959/71‑1 (MQ=255)
                                                     cATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCa               <  1:4975227/70‑1 (MQ=255)
                                                     cATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCAt              <  1:5055109/71‑1 (MQ=255)
                                                      aTAACTGGCACGTAGGGTTTGCAACCACTAACcca                                                 >  1:817344/1‑35 (MQ=255)
                                                      aTAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTcc                 <  1:3133098/67‑1 (MQ=255)
                                                      aTAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTcc                 <  1:1765872/67‑1 (MQ=255)
                                                       tAACTGGCACGTAGGGTTTGCAACCATTAACCCACCAATAGAGGGGTAGATAGGGCGt                         <  1:4131597/58‑1 (MQ=255)
                                                       tAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGt                         >  1:5208252/1‑58 (MQ=255)
                                                        aaCTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACt           >  1:568684/1‑71 (MQ=255)
                                                         aCTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                           >  1:2993213/1‑38 (MQ=255)
                                                           tGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                      <  1:1114710/41‑1 (MQ=255)
                                                            ggCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAggg                            >  1:2186694/1‑50 (MQ=255)
                                                            ggCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAggg                            >  1:1548937/1‑50 (MQ=255)
                                                             gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc      <  1:5537407/71‑1 (MQ=255)
                                                               aCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTg    <  1:4927887/71‑1 (MQ=255)
                                                                cGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTgg   >  1:3652166/1‑71 (MQ=255)
                                                                 gTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGt  >  1:1681499/1‑71 (MQ=255)
                                                                 gTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGt  >  1:4150523/1‑71 (MQ=255)
                                                                    |                                                                   
TAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGT  >  minE/1932695‑1932830

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: