Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,763 G→A 65.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,7630GA65.4% 82.2 / 73.4 81intergenic (‑134/+162)yqeF/lysSpredicted acyltransferase/lysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (26/27);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTT  >  minE/1932697‑1932831
                                                                  |                                                                    
aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGt                                                                   <  1:2020001/70‑1 (MQ=255)
aaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGt                                                                   <  1:2929275/70‑1 (MQ=255)
 aaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGttt                                                                 >  1:3686116/1‑71 (MQ=255)
 aaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGttt                                                                 >  1:3302701/1‑71 (MQ=255)
  aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTg                                                                >  1:227862/1‑71 (MQ=255)
   aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTg                                                                >  1:1820570/1‑70 (MQ=255)
   aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                               <  1:6872/71‑1 (MQ=255)
   aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGc                                                               <  1:3332355/71‑1 (MQ=255)
     aGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCaa                                                             <  1:3697110/71‑1 (MQ=255)
     aGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCaa                                                             >  1:2671831/1‑71 (MQ=255)
      gTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAAc                                                            <  1:238748/71‑1 (MQ=255)
       tGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAAcc                                                           <  1:396964/71‑1 (MQ=255)
          tttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACt                                                        <  1:1238681/71‑1 (MQ=255)
          tttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                        <  1:3594304/71‑1 (MQ=255)
             ttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTaa                                                      <  1:2443861/70‑1 (MQ=255)
             ttATCACAAAAGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAAc                                                     >  1:703995/1‑71 (MQ=255)
               aTCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAAcc                                                    <  1:3298916/70‑1 (MQ=255)
                   caAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACt                                                        <  1:444331/62‑1 (MQ=255)
                   caAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACccac                                                 <  1:711887/69‑1 (MQ=255)
                    aaaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGttt                                                                 >  1:1901056/1‑52 (MQ=255)
                     aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAAc                                                            >  1:3091731/1‑56 (MQ=255)
                     aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTaa                                                      >  1:3836110/1‑62 (MQ=255)
                      aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGt                                                                   <  1:831622/48‑1 (MQ=255)
                      aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCaa                                              >  1:1328981/1‑69 (MQ=255)
                      aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCaa                                              >  1:1069696/1‑69 (MQ=255)
                      aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACcca                                                  >  1:2115585/1‑65 (MQ=255)
                      aGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACcca                                                  >  1:914075/1‑65 (MQ=255)
                       ggAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATa                                            <  1:3314606/70‑1 (MQ=255)
                        gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCa                                                          >  1:271425/1‑55 (MQ=255)
                        gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATaga                                          <  1:1270837/71‑1 (MQ=255)
                         aaaTATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATagag                                         >  1:199428/1‑71 (MQ=255)
                            tatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCa                                                          <  1:1247822/51‑1 (MQ=255)
                                 cTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                  >  1:548645/1‑70 (MQ=255)
                                 cTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                 <  1:1121292/71‑1 (MQ=255)
                                  tGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                  >  1:1055590/1‑69 (MQ=255)
                                  tGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGa                                  >  1:2432192/1‑69 (MQ=255)
                                   gAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAg                               >  1:2123775/1‑71 (MQ=255)
                                   gAGCAGCAGTCAGAGACATAACTGGCACGTAAGGt                                                                   >  1:1135732/1‑35 (MQ=255)
                                   gAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTaa                                                      >  1:3210257/1‑48 (MQ=255)
                                    agcagcAGTCAGAGACATAACTGGCACGTAGGGTTTGCAAc                                                            >  1:3044318/1‑41 (MQ=255)
                                    agcagcAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACcca                                                  >  1:2342821/1‑51 (MQ=255)
                                    agcagcAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAgg                              >  1:2943157/1‑71 (MQ=255)
                                    agcagcAGTCAGAGACATAAATGGCACGTAAGGTTTGCAACCACTAACcca                                                  >  1:1907028/1‑51 (MQ=255)
                                     gcagcaGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAg                               >  1:2211150/1‑69 (MQ=255)
                                      cagcagTCAGAGACATAACTGGCACGTAGGGTTTGc                                                               <  1:205775/36‑1 (MQ=255)
                                        gcagTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGt                          <  1:3346316/71‑1 (MQ=255)
                                         cagTCAGAGACATAACTGGCACGTAAGGTTTGCaa                                                             <  1:3869854/35‑1 (MQ=255)
                                          agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAgggg                                      <  1:870562/57‑1 (MQ=255)
                                          agTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTa                        <  1:2413032/71‑1 (MQ=255)
                                            tCAGAGACATAACTGGCACGTAGGGTTTGCAACCACt                                                        <  1:2938119/37‑1 (MQ=255)
                                            tCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATag                                           <  1:1995093/50‑1 (MQ=255)
                                            tCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGtt                         <  1:981639/68‑1 (MQ=255)
                                            tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTa                                    <  1:140714/57‑1 (MQ=255)
                                            tCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTa                                    <  1:447114/57‑1 (MQ=255)
                                             cAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                 >  1:2834728/1‑59 (MQ=255)
                                             cAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGAt                                 >  1:3798679/1‑59 (MQ=255)
                                             cAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtc                     >  1:2829020/1‑71 (MQ=255)
                                              agagACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtct                    >  1:2451562/1‑71 (MQ=255)
                                              agagACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtc                     <  1:803121/70‑1 (MQ=255)
                                              agagACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                       <  1:490953/52‑1 (MQ=255)
                                              agagACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtct                    <  1:207220/71‑1 (MQ=255)
                                                   cATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCAt               <  1:3105313/71‑1 (MQ=255)
                                                     tAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAAtctc                   <  1:2601437/65‑1 (MQ=255)
                                                     tAACTGGCACGTAATGTTTGCAACCACTAACCCACCAATagag                                         <  1:2181103/43‑1 (MQ=255)
                                                     tAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAggg                                       <  1:719858/45‑1 (MQ=255)
                                                      aaCTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCg                           >  1:2804534/1‑56 (MQ=255)
                                                      aaCTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATAc             >  1:217817/1‑70 (MQ=255)
                                                       aCTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACt            >  1:290831/1‑70 (MQ=255)
                                                        cTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACtt           >  1:2966368/1‑70 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       >  1:1045653/1‑71 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       <  1:251644/71‑1 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       <  1:3245856/71‑1 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       >  1:1536392/1‑71 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       <  1:2216762/71‑1 (MQ=255)
                                                           gCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAAcc       >  1:3906007/1‑71 (MQ=255)
                                                            cACGTAAGGTTTGCAACCACTAACCCACCAATaga                                          >  1:737593/1‑35 (MQ=255)
                                                            cACGTAAGGTTTGCAACCACTAACCCACCAATAGAgggg                                      <  1:1068254/39‑1 (MQ=255)
                                                               gTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGt   <  1:2863169/71‑1 (MQ=255)
                                                               gTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGt   >  1:3380421/1‑71 (MQ=255)
                                                                tAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATa                                >  1:2694740/1‑41 (MQ=255)
                                                                tAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCCCCCATACTTAACCTGGtt  >  1:2417842/1‑71 (MQ=255)
                                                                  |                                                                    
AAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTT  >  minE/1932697‑1932831

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: